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ABSTRACT 


This  thesis  presents  the  results  of  an  evaluation 
and  extension  of  an  I.B.M.  digital  computer  program  designed 
for  control  systems  analysis. 

The  original  program  proved  useful  in  the  determina¬ 
tion  of  root  locus  diagrams,  and  the  modified  program  in  the 
determination  of  Bode  plots,  Nyquist  plots,  and  transient  res¬ 
ponse  of  control  systems. 

The  transient  response  calculation  uses  the  Z-trans- 
form  ability  of  the  original  program  plus  a  power  series  in¬ 
version.  This  method  is  a  modified  form  of  the  Boxer  and 
Thaler  method (1) ,  and  is  used  in  the  belief  that  it  provides 
more  overall  flexibility. 

The  transient  response  of  continuous  systems  can  be 
determined  using  this  program,  if  a  fictitious  sample  and  hold 
device  is  added  to  the  system  input.  A  lead  is  also  added  to 
counteract  the  lag  introduced  by  the  sample  and  hold  device. 
Thus,  for  the  continuous  system,  the  Z-transform  of  the  pro¬ 
duct  of  the  closed  loop  transfer  function,  the  sample  and  hold 
transfer  function,  and  the  fictitious  lead  is  calculated.  This 
Z-transfer  function  is  then  inverted  to  obtain  points  on  the 
transient  response  curve.  Sampled-data  systems  are  analyzed 
directly  without  the  addition  of  any  fictitious  lead  or  hold 


device . 
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Seven  examples  are  presented  to  illustrate  the  man¬ 
ner  in  which  the  program  can  be  employed  to  handle  different 
control  problems  in  the  design  or  analysis  such  as: 

1.  the  effect  of  pure  time  delay; 

2.  the  effect  of  sampling  an  analog  control  system; 

3.  the  design  of  digital  controllers;  and 

4.  the  design  of  analog  controllers. 
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1.  INTRODUCTION 

This  thesis  consists  primarily  of  the  results  of 
an  evaluation  and  extension  of  an  I.B.M.  digital  computer 
program  for  use  in  Control  Systems  Analysis.  The  investiga¬ 
tion  was  to  determine  the  program's  capabilities  and  limita¬ 
tions  as  applied  to  chemical  engineering  process  control  prob¬ 
lems  as  well  as  examine  the  program  for  possible  extension. 

It  was  thought  that  the  program  should  satisfy  the 
following  conditions: 

1.  It  should  provide  a  significant  contribution  to  the  over¬ 
all  systems  analysis. 

2.  It  should  be  user  orientated. 

3.  The  outcome  of  the  analysis  should  be  presented  in  the 
most  easily  understood  form. 

4.  Program  options  should  be  user  specified  so  that  the 
analysis  can  be  tailored  to  the  user's  specific  problem. 

5.  Computation  time  should  be  kept  to  a  minimum. 

A  controlled  system  can  be  represented  very  simply 
by  the  use  of  block  diagrams.  A  simple  feedback  control 
system,  in  block  diagram  form,  is  shown  in  Figure  1,  primarily 
as  an  introduction  to  the  terminology  used  throughout  this 
thesis . 

Referring  to  Figure  1,  the  part  of  the  control  system 
containing  G^(s),  that  is  the  system  between  1  and  2  in  Figure 
1,  will  subsequently  be  designated  as  the  ''forward  loop''  of  the 
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Figure  1 

\ 

General  Block  Diagram 
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Nomenclature 
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Math 

Model 
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G-^  ( s)  G2  ( s) 

0Q  ( s)  —  G^ ( 

L  L 


Characteristic  Equation 


1  +  G1(s)G2(s)  = 


0 
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system.  A  path  such  as  the  one  containing  G^Cs)  and  extend¬ 
ing  from  2  to  1  will  be  designated  as  a  "feedback  loop". 

G^(s) ,  ^2(3) ,  and  G^Cs)  will  be  referred  to  as  transfer  func¬ 
tions  and  are  defined  as: 


ei  (s) 

where  0q(s)  is  the  Laplace  transform  of  the  output  signal  of 
the  block  and  0^(s)  is  the  Laplace  transform  of  the  input  sig¬ 
nal  for  all  I.C.  of  0Q(t)  and  0^(t)  equal  to  zero. 

A  main  object  in  control  work  is  to  be  able  to 
analyze  a  system  such  as  that  shown  in  Figure  1,  for  its  dy¬ 
namic  characteristics.  From  this  study,  information  neces¬ 
sary  to  determine: 

1.  The  type  of  controller 

2.  The  control  parameter  settings 
should  be  derived  so  that 

1.  The  stability  of  the  system  is  satisfactory 

2.  The  control  criteria  are  met. 

A  method  commonly  used  for  analyzing  a  system  for 
stability  and  the  one  upon  which  this  control  system  analy¬ 
sis  program  is  based  is  the  Root  Locus  Technique.  A  Root 
Locus  is  defined  to  be  the  trajectory  in  the  "s"-plane  fol¬ 
lowed  by  a  root  of  the  characteristic  equation  as  some  para¬ 
meter  of  the  corresponding  system  is  varied  continuously.  The 
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parameter  which  is  varied  in  this  analysis  is  the  overall  loop 
gain  K.  In  Figure  1  the  characteristic  equation  for  the  system 
is  1  +  G^(s)G2(s)  =  0  where  K  is  included  in  one  of  the 

transfer  functions,  but  is  usually  separated  for  the  purpose 
of  a  root  locus  analysis. 

This  characteristic  equation  determines  the  transient 
response  of  a  system  to  any  input  and  thus  the  location  of  its 
roots  can  be  used  to  indicate  stability  or  instability  for 
given  values  of  K. 

Analysis  of  control  systems  using  frequency  response 
or  root  locus  techniques  has  generally  been  a  laborious  task 
and  because  of  this  a  really  complete  investigation  is  seldom 
done.  With  the  digital  computer  it  is  now  possible  to  greatly 
reduce  the  time  expended  in  performing  these  analyses  and  thus, 
greatly  increase  the  extent  of  their  application.  The  ulti¬ 
mate,  of  course,  is  to  have  the  analysis  and  design  calcula¬ 
tions  performed  completely  by  the  machine  with  only  the  sub¬ 
jective  decisions  left  to  the  user. 

The  investigation  was  conducted  with  this  ultimate 


goal  in  mind. 
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2 .  LITERATURE  SURVEY 

This  thesis  deals  with  the  role  of  the  digital  com¬ 
puter  in  the  formation  and  use  of  the  various  control  system 
stability  analysis  aids,  such  as  the  Bode,  Nyquist,  Root  Locus 
and  Transient  Response  curves. 

Texts  found  to  be  of  particular  use  in  outlining 
the  basic  theory  behind  the  Root  Locus,  Nyquist  and  Bode  Plot 
methods  of  analysis  were  those  by  Evans(4)  and  Schilling ( 18 ) . 

A  summary  of  papers  dealing  with  frequency  response 
methods  is  presented  by  Oldenburger ( 16 ) .  This  is  a  collection 
of  basic  papers  in  the  field  and  provides  an  excellent  list  of 
references  for  further  search. 

The  data  derived  from  frequency  response  studies  con¬ 
sists  of  different  values  of  amplitude  ratio  and  phase  angle 
between  the  input  and  output  of  a  sinusoidally  forced  system 
for  different  frequencies.  This  type  of  data  lends  itself  to 
analysis  by  the  Bode  and  Nyquist  methods  but  it  cannot  be  used, 
as  is,  in  a  Root  Locus  analysis,  since  this  type  of  analysis 
depends  on  having  the  transfer  function  of  the  system.  Some 
methods  have  been  worked  out  to  obtain  the  transfer  functions 
for  certain  systems  from  this  type  of  data.  A  graphical  method 
is  presented  by  Ganapathy  and  Krishna (6)  and  a  numerical  method 
by  Staffin  and  Staff in (20).  An  analytical  method  for  obtaining 
a  root  locus  of  a  given  characteristic  equation  is  outlined  by 
Chang  (2),  and  the  analogue  computer  is  used  in  conjunction  with 


' 
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the  generalized  Mitrovic  method  in  a  procedure  outlined  by 
Kokotovic  and  Siljak(ll). 

The  above  methods  are  means  of  analyzing  a  system 
for  its  transient  response  characteristics  without  actually 
arriving  at  the  true  transient  response  solution.  Wolfang, 
Wagner  and  Zoss(22)  have  presented  a  method  whereby  the  tran¬ 
sient  response  of  a  system  is  obtained  from  frequency  response 
plots  using  discrete  point  data.  This  method  appears  suitable 
for  digital  computer  application.  A  means  of  finding  the 
transient  response  of  a  system  given  the  input  using  infinite 
matrices  in  time  and  sampled-data  theory  is  outlined  by  Dorf(3). 
This  method  is  particularly  suited  to  digital  computer  appli¬ 
cations  as  it  involves  matrix  inversion,  addition,  and  multi¬ 
plication. 

Digital  computers  are  beginning  to  make  themselves 
felt  as  control  instruments,  for  their  versatility  as  well  as 
their  capability  of  handling  a  large  number  of  control  loops. 

A  discussion  of  the  advantages  and  disadvantages  of  Direct 
Digital  Control,  the  economic  considerations,  and  minimum 
sampling  periods  found  to  be  effective  for  certain  types  of 
systems  is  presented  by  Klock  and  Schoef f ler (10) . 

With  the  digital  computer  arrives  the  so-called 
sampled-data  system,  and  the  extension  of  the  methods  applied 
to  the  continuous  system,  specifically  the  root  locus  method, 
involves  the  theory  of  the  Z-transform.  An  outgrowth  of  this 


I 
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is  the  use  of  the  Z-transform  in  predicting  the  transient  res¬ 
ponse  of  a  system  to  a  disturbance.  Work  has  been  done  along 
this  line  by  Fowler (5)  and  by  Boxer  and  Thaler (1). 

To  use  the  Z-transform,  a  basic  knowledge  of  it  is 
necessary  and  this  can  be  gained  from  texts  by  Jury (9),  Lin- 
dorff (13) ,  and  Kuo (12) . 

Monroe (14)  has  published  a  text  which  is  primarily 
concerned  with  the  synthesis  of  digitally  controlled  systems. 

The  various  factors  which  aid  in  the  selection  and  design  of 
a  digital  controller  are  covered.  Tou(21)  in  his  paper  ap¬ 
proaches  the  design  of  a  "dead-beat"  controller  for  a  sampled- 
data  system  using  the  state-variable  technique.  A  complete 
analysis  of  a  sampled-data  system  using  the  Z-transform  is 
given  in  a  paper  by  Slaughter ( 19) .  The  analysis  is  very  com¬ 
plete  and  shows  the  use  of  the  Z-plane  root  locus  in  design. 

Some  actual  control  problems,  test  equipment  and 
the  methods  used  to  come  to  a  solution  are  presented  by  Rock (17). 
This  is  generally  a  review  of  testing  equipment  and  the  infor¬ 
mation  which  can  be  gained  from  each.  Murrill  and  Smith (15) 
have  written  a  paper  on  the  importance  of  correct  controller 
settings  and  some  present  rules  of  thumb  used  to  arrive  at 
these  settings. 
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3.  Z -TRANSFORM  THEORY 

3 . 1  The  Z-Transform 

The  Control  Systems  Analysis  Program (7)  was  origin¬ 
ally  designed  for  use  as  a  convenient  method  for  obtaining  a 
root  locus  in  either  the  s  or  Z-plane.  This  capability  allows 
its  use  in  analyzing  various  types  of  control  systems  such  as 
some  sampled-data  systems,  systems  with  pure  time  delay,  linear 
systems,  or  a  combination  of  pure  time  delay  and  sampled-data, 
the  latter  system  containing  a  pure  time  delay  along  with  a 
sampling  device.  One  of  the  objectives  of  the  design  of  this 
program  was  to  have  it  analyze  sampled-data  systems  and  pure 
time  delay  systems.  To  do  this,  a  routine  to  convert  an  s- 
transfer  function  to  its  corresponding  Z-transfer  function 
was  worked  out  and  included  in  the  program. 

There  are  a  number  of  classical  methods  by  which  a 
Z-transform  can  be  calculated  from  the  Laplace  Transform,  such 
as  the  method  of  complex  convolution,  by  partial  fractions 
where  the  function  is  broken  down  into  partial  fraction  identi¬ 
ties  which  can  then  be  transformed  directly  to  the  Z-form,  or 
by  the  direct  numerical  evaluation  of  the  series, 

G(Z)  = 

where  wq  is  the  sampling  frequency. 

The  computation  required  for  a  Z-transform  of  a  com¬ 
plex  system  by  any  of  the  above  becomes  very  tedious,  even  for 


1  n=°° 

—  y  G(s+inw  ) 
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machine  computation  where  applicable.  Because  of  this,  a  dif¬ 
ferent  method (8)  of  computing  the  Z-transform  from  an  s-trans- 
form  was  incorporated  into  this  program  by  which  is 

particularly  suited  to  machine  calculation. 

If  the  degree  of  s  in  the  numerator  of  the  transfer 
function  G(s)  is  two  or  more  less  than  that  of  the  denominator, 
equation  (1)  can  be  written  as 

G(Z)  = 

or  in  another  form 

N(Z) 

D  ( Z ) 

The  order  of  the  denominator  of  G(Z)  must  be  the 

same  as  that  of  the  denominator  of  G(s).  This  correspondence 

is  realized  from  the  conformal  mapping  of  the  left  hand  side 

of  the  s-plane  into  a  unit  circle  centered  at  the  origin  in 

sT 

the  Z-plane  when  the  transformation  Z  =  e  is  used.  This 
transformation  only  maps  the  poles  onto  a  different  plane 
so  that  the  number  of  roots  remains  unchanged.  Thus,  D(Z), 
which  is  the  denominator  of  the  Z-transfer  function  can  easily 
be  computed  by  a  direct  transformation  of  the  poles  of  G(s) 
in  the  s-plane  onto  the  Z-plane.  With  D(Z)  known,  equation 
(3)  can  be  written  in  the  form 

(4) 


D(Z) 


n=+°° 


n=-< 


G(s  +  inw  ) 
o 


1  n=+°° 

—  I  G(s  +  inw  ) 
_  L  o 

T  n=-°° 


(2) 


Z=e 


sT 


1  n=+°° 

—  y  G(s  +  inw  ) 
u  o 

T  n=-°° 


(3) 


Z=e 
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N(Z) 
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where  N(Z)  is  an  unknown  Z-polynomial  and  is  assumed  to  be  of 
the  same  degree  as  D(Z) .  The  order  of  N(Z)  cannot  be  greater 
than  that  of  D(Z)  for  a  physically  realizable  system.  However, 
N(Z)  may  be  of  a  lesser  degree  than  D(Z)  and  for  this  case,  the 
coefficients  for  those  assumed  powers  of  Z  which  do  not  exist 
will  turn  out  to  be  equal  to  zero  in  the  computer  program  solu¬ 
tion. 

For  any  given  complex  value  of  s,  the  complex  func¬ 
tion  represented  by  the  right  hand  side  of  equation  (4)  can 
be  evaluated  and  put  in  the  form  of  a  complex  number,  a  +  16. 
Then,  using  n+1  even  different  values  of  s  where  n  refers  to 
the  degree  of  the  Z-transform  and  n  is  an  odd  integer,  a  set 
of  n+1  equations,  such  as  the  following  can  be  formed: 


c 

o 

+ 

ciYi 

+ 

C2y12 

+ 

o.o  + 

II 

£ 

1 — 1 
>- 
G 

U 

al 

0 

+ 

ci6i 

+ 

C26!2 

+ 

.  o  .  + 

C  6  n  = 
n  1 

61 

c 

o 

+ 

C1y2 

+ 

C2y22 

+ 

...  + 

C  y"  = 
n '  2 

a2 

6  Ll  6n 

0  +  C,  +  ...  +  C  — 2-i- 

1  n 


n+1 


2  2  2 
where  Z  =  y  +  i<$  r  and  yn  and  6n  represent  the  real  and  ima¬ 
ginary  parts  of  Zn.  The  subscripts  of  the  y's,  6's,  a's,  and 
B's  correspond  directly  to  the  different  values  of  s  for  which 
they  were  calculated.  If  n  were  an  even  integer,  the  last 
equation  in  (5)  would  be: 


. 


, 
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c  +  c_ 

o  1 


Y 


n+2 


+  ...  +  C 


Y 


n 


n+2 


a 


n+2 


n 


(6) 


Solving  this  set  of  equations  for  the  constants , , 
will  then  define  the  numerator  N(Z)  as: 

N(Z)  =  Cq  +  C1Z  +  C2Z2  +  C3Z3  +  ...  +  CnZn  (7) 


Thus,  the  Z-transform  of  the  system  G(s)  is 


G  ( Z )  = 


C  +  C  z  +  C  z2  +  ...  +  c  zn 
o  1_ 2_ n 

D  ( Z ) 


(8) 


The  root  locus  can  then  be  obtained  for  the  sampled 
system  either  in  the  s  or  Z-plane  from  equation  (8). 


3.1.1  The  Degree  of  the  Numerator  One  Less  Than  That 

of  the  Denominator 


This  is  a  special  case  where  the  term  e(0+)  in 
equation  (1)  is  not  equal  to  zero.  For  this  eventuality  it 
is  necessary  to  make  a  change  in  equation  (2).  If  one  had  a 
transform  such  as: 


G  ( Z )  = 


I  (z) 


C  +  C-|  z  +  . . .  +  c  z 

o_ 1_ n 

D  +  D-.Z+  ...  +  DZ 
o  1  n 


n 


n 


(9) 


its  equivalent  difference  equation  (22)  would  be 


D  0  (0)  +  D  -.0  (-1)  +  ...  +  D  0(-n)  =  C  1(0) 

no  n-1  o  on 


+  C  I  (-1)  +  .  .  .  +  C  l  ( -n) 
n-l  o 


(10) 


/ 


12 


where  0  ( — i)  represent  the  values  of  the  output  function  at 

past  sampling  instants  and  the  -i  in  parentheses  determines 
the  sampling  instant.  The  term  e(0+)  in  equation  (1)  represents 
a  sampled  value  of  the  input  function  at  time  zero  and  the 
only  term  in  equation  (10)  which  could  represent  this  is 
C  1(0),  therefore  must  equal  e(0+).  A  corrected  form  for 
equation  (4)  for  the  case  of  the  numerator  degree  being  one 
less  than  that  of  the  denominator  is 

D  ( Z )  n=+°°  D  ( Z ) 

N  ( Z )  =  -  T  G(s  +  inw  )  +  C  -  (11) 

m  L  o  n  0 

T  n=-°°  2 


Thus,  a  modification  to  the  set  of  equations  (5)  is  neces¬ 
sary  in  order  to  include  this  extra  term.  When  this  term  is 
included  and  equation  (1)  is  expanded  and  evaluated  in  the 
same  manner  as  equation  (4)  we  get  a  system  of  equations  in 
the  form: 


C  +  C-,  y +  ...  +C  y-i  =  a + 
o  1 ' 1  n '  1  1 


CnD(V 


cnD ( 6 i ) 

0  +  C.6,  +  ...  +  C  «1  =  B1  +  — - — 


(12) 


which,  on  collecting  terms  can  be  put  into  the  form 

n  D(y1) 

C  +  C-iYt  +  •••  +  C  (y-|  ~ 

oil  n  1  2  i 

D  (  5 .  ) 


(13) 


0  +  C,6,  +  ...  +  C(6n 
11  n  1 


n 


-)  =  6, 


2 


-  • 
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where  D(y^)  and  D(6^)  represent  the  real  and  imaginary  parts 
of  D  ( Z )  with  D  (Z )  calculated  for  s  =  s^. 

3 . 2  Z-Transform  Inversion 

There  are  several  methods  for  inverting  Z-transforms 
to  arrive  at  the  time  solution.  The  Partial  Fraction  method 
and  the  Inversion  Formulae  method  are  particularly  suited  for 
hand  calculation,  whereas,  the  Power  Series  method  is  better 
siited  for  machine  calculation. 

The  Partial  Fraction  method  involves  breaking  a  Z- 
transform  into  its  partial  fractions  and  using  a  Z-transform 
table  to  find  the  corresponding  time  functions.  This  method 
is  probably  the  best  to  use  if  the  transfer  function  has  a 
pole  with  a  non-zero  imaginary  part.  The  Inversion  Formulae 
is  basically  another  hand  method  which  is  most  easily  applied 
to  transfer  functions  with  simple  poles  on  the  real  axis. 

The  Power  Series  method  involves  a  continuous  division  of  the 
numerator  of  the  Z-transform  by  its  denominator  so  that  a 
power  series  in  Z  ^  results.  In  this  method,  if  the  transformed 
function,  E(Z),  includes  the  input  signal 

N(Z) 

E  ( Z )  =  GI(Z)  =  -  (14) 

D  ( Z ) 

the  coefficients  of  the  power  series  represent  the  value  of 

t  h 

the  time  function  e(nT)  at  the  n  sampling  instant.  To  show 


. 
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this  let: 


E(Z)  = 


U  Z  4*3.  ,Z  +  ...+U..ZH"e 
m _ m-1 _ 1 _ o 

b  Zn  +  b  nZn~1  +  ...  +  b,  Z  +  b 
n  n-1  1  o 


(15) 


then,  dividing  the  numerator  by  the  denominator  by  long  divi¬ 
sion  will  yield: 


E  (Z)  =  AZm  n  +  BZm  n  1  + 


(16) 


where  n  >_  m  for  the  system  to  be  physically  realizable.  Now, 
since  the  Laplace  Transform  of  a  sampled  time  function  can  be 
shown (12)  to  be: 


(t) 


=  E*  ( s)  =  l  e(nT)e 


-nTs 


(17) 


n=0 


and 


,n 


=  e 


nTs 


(18) 


it  can  clearly  be  seen  that  A,  B,  ...  correspond  to  the  res¬ 
pective  values  of  e (nT) .  If, however,  E(Z)  does  not  include 
the  input  but  can  be  represented  as: 

9  (Z) 

E  (  Z  )  =  G  (  Z  )  =  — - 

I  (Z) 

then  the  coefficients  cannot  represent  the  value  of  the  time 
function  at  the  sampling  instants  but  only  some  ratio  of 
Output  to  Input.  In  order  to  arrive  at  values  for  the  output 
time  function,  one  must  first  obtain  the  output  as  some  func¬ 
tion  of  the  input  and  the  transfer  function. 


. 
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Given  the  system 


I(s)  >(  I*(s) 

T 


l - * 

I  T 


6*(s) 


6(s) 


Then 


0(s)  =  G  ( s)  I*(s) 

In  terms  of  the  variable  e*(s) 


(19) 


0*(s)  =  G*(s)  I*(s) 

Now  following  equation  (17) 


Similarly 


and 


G*(s) 


l  G(nT) 
n=0 


-nTs 

e 


I*  (s) 


l  I(nT) 
n=0 


-nTs 

e 


9*  (s) 


l  0 (nT) 
n=0 


-nTs 

e 


(20) 


(21) 


(22) 


(23) 


Now  for  a  system  composed  of  linear  time-invariant  elements 
(25)  ,  the  response  to  a  unit  impulse  applied  to  the  system 
input  at  t  =  0  is 


G*  ( t)  =  l  G  (nT)  6  (t-nT)  (24) 

n=0 

For  an  impulse  modulated  signal 


n=+°° 

=  l  I  (nT)  (5  (t-nT) 

n=-oo 


I*  (t) 


(25) 
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If  the  system  is  linear,  then,  superposition  applies  and  the 
sampled  response  to  an  impulse  train  is  given  by 


0* 


o 


(t) 


l 


n=-«> 


l  G(nT)  I(nT-mT) 
m=0 


6 (t-nT) 


(26) 


The  output  sample  sequence  which  is  of  interest  is  represented 
by  the  bracketed  term  of  equation  (26)  so  that 


0 (nT)  =  l  G(mT)  I(nT-mT)  (27) 

m=0 

or 

oo 

0  (nT)  =  l  G(mT)  6  (t-mT)  I  (nT)  (28) 

m=0 

Now  transforming  equation  (20)  yields 


0(Z)  =  G  ( Z )  I  ( Z ) 


(29) 


Now  G(Z)  can  be  represented  as  an  infinite  series 


of  the  form 


G  ( Z )  =  A±Z  1  +  A2Z  2  +  A3Z  3  +  ... 


(30) 


From  equation  (30)  it  can  be  seen  that  G(mT)  in 
equation  (28)  can  be  replaced  by  A^  for  0  <  n  <_  °°  thus: 


0  (nT)  =  V  A  6  (t-mT)  I  (nT) 
L  r,  m 

m=0 


where  A  are  the  coefficients  of  G(Z). 
m 


(31) 
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This  shows  then,  that  when  using  the  transfer  func¬ 
tion  to  arrive  at  a  value  of  the  time  function  0 (nT) ,  it  is 
necessary  to  sum  the  product  of  Am  and  I (nT)  for  m  going  from 
zero  to  n  instead  of  simply  the  coefficients  of  the  Z-trans- 
form  being  the  values  of  the  time  function  at  that  sampling 
instant . 


. 
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4.  THE  ROOT  LOCUS 


The  control  systems  analysis  program  does  the  cal¬ 
culations  necessary  to  establish  the  root  locus.  The  control 
system  is  specified  by  its  feed-forward  and  feedback  loop 
transfer  functions.  Points  on  the  root  locus,  defined  by 
G(s)/s=a+ib  =  -1/K,  are  found  by  a  systematic  grid  search 
and  interpolation  procedure  over  a  user  specified  area  in 
the  s-plane. 

The  Root  Locus  Technique  is  based  on  satisfying 
the  characteristic  equation  of  the  system  and  any  point  satis¬ 
fying  this  equation  is  a  point  on  the  root  locus.  The  char¬ 
acteristic  equation  can  be  expressed  as: 

1  +  K  G  ( s)  =  0  (32) 

where  G(s)  is  the  open  loop  transfer  function  of  the  system, 
and  K  is  the  loop  gain.  Rearranging  we  arrive  at 

G (s)  =  -  |  (33) 

Now  for  various  values  of  s,  s  =  6  +  iw,  values  for 
G(s) ,  G(s)  =  q  =  u  +  iv,  can  be  calculated.  Since  -1/K  has 
to  be  a  real  number,  v  must  equal  zero  for  equation  (33)  to 
be  satisfied.  The  procedure,  then,  is  to  calculate  values  of 
G(s)  for  various  values  of  s  along  a  line  of  constant  6  or 
w  in  the  complex  plane  and  note  the  sign  of  v.  When  a  change 
in  the  sign  of  v  occurs,  it  signifies  that  the  root  locus  has 
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been  crossed.  The  point  in  the  s-plane  which  satisfies  equa¬ 
tion  (33)  can  then  be  located  more  accurately  by  a  search  be¬ 
tween  grid  points.  The  loop  gain  corresponding  to  this  value 
of  s  is  then  calculated  from 

K  =  -  —  (34) 

u 

This  type  of  procedure  is  carried  out  over  the  whole  specified 
scan  area. 

The  method  is  similar  to  the  principle  incorporated 
in  Evan's  Spirule,  but  is  suited  to  use  by  a  high-speed  digi¬ 
tal  computer.  Once  the  root  locus  has  been  determined  as  a 
list  of  points,  from  which  a  graph  may  be  obtained,  the  re¬ 
sults  can  be  examined  to  determine  system  stability,  (A. 9. 2). 

The  criterion  which  is  used  in  this  type  of  analysis 
is  that  for  a  given  value  of  loop  gain  K,  all  the  roots  of 
the  characteristic  equation  must  exist  in  the  left  half  of  the 
s-plane.  This  guarantees  that  all  the  roots  have  negative 
real  parts  and  thus  the  time  solution  in  response  to  a  bounded 
disturbance  will  not  become  unbounded  for  increasing  time. 

This  is  defined  as  a  stable  system. 

The  values  of  K,  for  which  the  real  part  of  each 
root  of  the  characteristic  equation  is  less  than  zero,  can 
be  determined  from  either  the  complex  Z  or  s-plane  for  a 
sampled-data  system.  Generally,  the  Z-plane  root  locus  is 
almost  exclusively  used  for  analysis  for  sampled-data  systems 
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and  the  s-plane  root  locus  for  continuous  systems.  The  Control 
Systems  Analysis  program,  however,  has  made  the  s-plane  avail¬ 
able  for  the  stability  analysis  of  both  types  of  system,  thus, 
the  user  does  not  have  to  be  familiar  with  both.  To  arrive  at 

the  s-plane  root  locus  for  a  sampled-data  system  the  follow¬ 
ing  operation  sequence  is  performed. 

a)  The  Z-transform  is  calculated. 

b)  The  first  scan  point  in  "s"  is  determined  according  to 
user  specification. 

c)  This  scan  point  is  converted  to  its  equivalent  in  Z. 

d)  A  system  number  in  the  Z-plane  is  calculated. 

e)  This  system  number  is  examined  for  its  proximity  to  a 
root  locus . 

f)  If  it  is  on  the  root  locus,  the  point  is  printed  out 
along  with  its  value  in  the  s-plane.  If  it  is  not,  the 
next  scan  point  is  selected  in  s  and  the  process  re¬ 
peated  from  b) . 

There  are  some  undesirable  characteristics  to  this 
procedure;  the  scan  area  in  Z  is  determined  by  that  specified 
in  the  s-plane.  This,  in  some  cases,  can  by  very  undesirable. 
Points  in  the  Z-plane  are  determined from  the  s-plane  by  the 
mapping  function  Z  =  e  .  By  using  scan  points  determined 
in  the  s-plane,  it  can  be  seen  that  the  area  in  the  Z-plane 
over  which  the  scan  takes  place  is  severly  limited.  For 
example,  using  this  program  Z  can  be  evaluated  only  in  the 
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first  quadrant.  In  order  to  reach  zero  on  the  Z-plane,  s 
would  have  to  be  given  the  value  of  minus  infinity.  This 
severely  limits  analysis  of  a  system  in  the  Z-plane. 

It  has  been  mentioned  that  the  root  locus  of  a 
sampled-data  system  can  be  converted  to  the  s-plane.  It  might 
be  thought  that  this  should  be  enough  for  a  good  stability 
analysis,  however,  analysis  of  a  sampled  system  in  the  s-plane 
does  hide  certain  influences,  particularly  the  full  influence 
of  the  sampling  interval  T. 

It  was  noted  in  the  transient  response  calculations 

for  continuous  systems,  that  very  small  sampling  intervals 

caused  multiple  poles  at  Z  =  1.  The  occurrence  of  which  led 

sT 

to  system  instability.  This  is  obvious  since  the  Lim  (e  )  =  1. 

T-*0 

This  instability  problem  does  not  usually  occur  in  physical 
systems  because  of  the  large  sampling  intervals  used.  The 
s-plane  shows  the  instability  of  the  system  in  this  case,  but 
it  is  hard  to  find  the  reason  until  the  root  locus  of  the  sys¬ 
tem  is  also  examined  in  the  Z-plane.  The  reason  for  the  in¬ 
stability  problem  then  becomes  obvious. 

In  order  to  avoid  the  problems  given  above,  a  pro¬ 
gram  change  was  made  which  now  provides  a  Z-plane  root  locus 
determination  over  a  scan  area  set  in  the  Z-plane.  This  is 
an  optional  calculation  and  requires  a  flag  set  by  the  user. 

With  this  flag  set  only  the  Z-plane  root  locus  is  determined. 
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5.  THE  CONTROL  SYSTEMS  ANALYSIS  PROGRAM 
5 . 1  The  Original  Program 

This  program  was  originally  capable  of  calculating 
a  root  locus  for  a  linear  feedback  system,  with  or  without 
time  delay,  and  also  capable  of  calculating  a  Z-transform  of 
a  transfer  function  which  was  specified  in  Laplace  transform 
notation.  Using  the  Z-transform  conversion  ability  it  chould 
then  provide  the  root  locus  in  the  s-plane  for  sampled-data 
systems . 

In  order  for  the  program  to  function,  it  was  neces¬ 
sary  to  arrange  the  block  diagram  of  the  system  to  be  analyzed 
into  forward  and  feedback  loops.  When  this  was  accomplished, 
the  transfer  functions  were  entered  with  those  in  the  feed¬ 
back  loops  entered  first  and  those  of  the  forward  loops  last. 
The  program  would  then  evaluate  and  combine  these  transfer 
functions  to  arrive  at  an  overall  system  number  for  that 
value  of  s  considered. 

This  aspect  of  the  program  has  a  limitation  in 
that  the  components  of  the  block  diagram  must  be  manipulated 
by  the  user  in  order  to  achieve  the  necessary  input  form. 

As  this  limitation  was  not  insurmountable,  only  an  enquiry 
into  the  programming  necessary  to  provide  this  function  for 
the  user  was  made.  Various  programming  languages  could  be 
used  in  this  capacity,  however,  no  attempt  was  made  to 
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implement  this,  as  the  main  concern  was  the  calculating  abil¬ 
ity  of  the  program,  and  the  types  of  systems  and  problems  for 
which  the  program  could  be  used. 

The  original  program  consisted  of  15  subroutines 
and  the  mainline  program  DK1.  These  subroutines  and  their 
relation  to  each  other  are  shown  in  Figure  (2) . 

Each  subroutine  has  one  specific  purpose  or  step 
in  the  overall  calculation  to  perform  and  can  be  called  a 
number  of  times  by  other  subroutines  to  perform  its  calcula¬ 
tion.  A  list  of  these  subroutines  with  their  purposes  is 
provided  in  Table  (1) . 

5 . 2  The  Present  State  of  the  Program 

During  the  investigation  of  the  C.S.A.  program  it 
was  seen  that  the  Bode  and  Nyquist  stability  analysis  plots 
could  be  easily  incorporated.  This  was  done  and  the  neces¬ 
sary  subroutines  tested  using  the  example  problem  from 
Schilling (18) . 

A  method  was  also  established  whereby  the  transient 
response  of  certain  systems  could  be  obtained  using  the  Z- 
transform  capability  of  the  program  as  a  basis.  The  method 
by  which  this  is  done  was  outlined  previously  in  Section  3.2. 
This  method  depends  on  the  program's  ability  to  calculate  Z- 
transforms  for  transfer  functions  in  s.  The  Z-transform  pro¬ 
gram  requires  that  the  denominator  of  the  s-transfer  function 
be  factored  into  first  and  second  order  terms.  A  Share  Library 
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TABLE  1 

Routine 

Name  Purpose 

DK1 

To  read  and  organize  the  flow  of  data  accord¬ 
ing  to  the  control  card. 

SCAN 

Computes  the  s-plane  scan  points  and  contains 
the  logic  required  to  identify  the  crossing 
of  a  root  locus. 

BITE 

Computes  the  open  loop  system  number  for  a 
set  of  forward  and  feedback  loops. 

SAMP  2 

Computes  independent  complex  system  numbers 
for  two  sets  of  transfer  functions  which  are 
then  multiplied  together.  For  instance,  it 
does  the  computation  necessary  for  the 
evaluation  of 

D  ( Z )  n=a  2tt 

-  £  G(s  +  in  — )  . 

T  n=-a  T 

ZCO 

Computes  the  Z-transform  of  a  first  or 
second  order  polynomial  in  s. 

ARITH 

Used  with  BITE  to  compute  system  numbers. 

The  required  arithmetic  is  carried  out. 

SAMP 

Used  to  compute  G*(s).  System  number  equals 
a  +  i3.  ERROR  is  also  checked  in  this  sub¬ 
routine  . 

POWER 

Computes  the  powers  of  s  and  Z. 

CDIV 

Divides  complex  numbers. 

CMUL 

Multiplies  complex  numbers. 

TRANS 

Transfers  desired  data  changes  from  input 
locations  to  the  proper  locations  when  mul¬ 
tiple  runs  are  specified  on  the  control  card. 

OUTPUT 

Contains  all  the  output  statements  except 
those  error  statements  such  as  SINGULAR  and 
ERROR. 

POLY 

Multiplies  polynomials  with  numerical  coef¬ 
ficients.  Resultant  degree  must  be  <.35 . 
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Routine 

SCANZ 

SOLVE 


TABLE  1  (continued) 


Name 


Purpose 


Determines  the  set  of  equations  to  be 
solved  to  get  the  numerator  coefficients. 

Solves  the  algebraic  equations  set  up  in 
SCANZ  by  the  Gauss  Elimination  Method. 
Will  print  ,!  singular"  if  the  determinant 
is  singular. 
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Program(23)  was  used  to  calculate  the  roots  for  a  polynomial 
greater  than  second  order  and  these  roots  were  then  used  as 
input  data  for  the  C.S.A.  program.  The  Share  program  is  not 
incorporated  into  the  original  Control  Systems  Analysis  pro¬ 
gram,  however,  this  could  be  done. 

A  modification  was  made  to  the  Control  Systems 
Analysis  program  so  that  a  Root  Locus  for  a  sampled-data  sys¬ 
tem  could  be  calculated  and  plotted  in  the  Z-plane.  The  modi¬ 
fication  consists  of  another  calculation  path  in  which  the 
Z-plane  is  scanned  in  the  same  manner  as  the  s-plane  when 
finding  a  root  locus  for  a  continuous  system.  In  order  to 
implement  this,  a  flag  is  set  by  the  user  and  the  scan  area 
in  the  Z-plane,  over  which  the  user  is  interested,  is  set  in 
the  same  manner  as  that  in  the  s-plane  for  a  continuous 
system. 

Provision  was  made  for  the  plotting  of  these  calcu¬ 
lated  Root  Locus  points  using  the  University  of  Alberta  Com¬ 
puting  System's  "Autoplot"  routine. 

The  more  complete  Z-plane  root  locus  has  been  found 
particularly  useful  in  the  design  of  sampled-data  or  digital 
controllers  for  given  systems.  This  aspect  will  be  treated 
later . 

The  new  calculation  path  designed  to  give  the  Z- 
plane  root  locus  was  checked  for  proper  functioning  by  the 
comparison  of  the  digital  solution  and  that  presented  in  an 
article  by  Slaughter (19)  for  the  same  problem. 
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Keeping  the  user  of  the  program  in  mind,  it  was 
thought  that  a  more  explanatory  form  of  Input  Data  printout 
should  be  available.  Consequently  another  output  routine, 
"DPRINT",  was  written  which  prints  out  the  control  card 
specifications  along  with  their  explanation.  This  allows 
the  user  to  make  an  immediate  check  of  the  input  data  without 
having  to  refer  back  to  the  literature. 

A  main  step  in  the  calculation  of  the  transient  res¬ 
ponse  for,  or  the  study  of,  a  system  controlled  with  a  sample- 
data  controller  was  the  formation  of  the  closed  loop  transfer 
function  (Section  7.4).  The  formation  of  this  transfer  func¬ 
tion  required  two  of  the  three  runs  usually  necessary  for  the 
solution,  in  addition  to  the  hand  multiplication  and  addition 
of  polynomials  with  the  consequent  loss  in  accuracy.  The  pro¬ 
gram  has  been  modified  so  that  only  two  runs  of  the  program 
are  necessary  and  all  the  multiplication  and  addition  of 
polynomials  is  done  by  the  machine. 

Figure  (3)  and  (4)  constitute  a  flow  diagram  show¬ 
ing  the  calculation  paths  of  the  present  modified  program. 

The  added  calculation  paths  and  blocks  are  indicated  by  dashed 


lines . 
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Figure  3  -  Calculation 
Flow  Diagram  of  Modified 
C.S0A0  Program 
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6.  PROBLEMS  FOR  WHICH  THE  PROGRAM 
IS  NOT  APPLICABLE 


All  problems  submitted  to  this  program  must  have 
all  their  transfer  functions  specified  in  either  the  Laplace 
or  Z-transform  notation  and  therefore  it  will  only  handle 
linearized  problems „ 

Some  other  types  of  problems  cannot  be  treated  by 
this  program.  For  the  case  where  there  is  a  pure  time  delay 
either  in  a  continuous  or  sampled-data  system,  the  root  locus 
can  be  determined,  the  transient  response  cannot.  The  dif¬ 
ference  lies  in  the  fact  that  the  root  locus  is  determined  by 
the  open  loop  and  the  transient  response  by  the  closed  loop. 

A  pure  time  delay  in  a  loop  is  treated  by  converting 
it  to  Z-notation,  (Appendix  A. 9.2).  For  a  continuous  system 

with  one  time  delay  of  length  T,  there  would  be  only  one 
sT 

Z  =  e  .  For  a  sampled  system,  the  Z  representing  the  delay 

sT 

would  be  different  than  the  Z  =  e  defined  by  the  sampling 

period  T.  In  both  cases,  the  difficulty  is  the  same,  the 

different  terms  cannot  be  separated  in  the  closed  loop. 

In  Figure  (5)  such  a  system  is  presented,  where 

e  represents  the  time  delay.  It  can  be  seen  that  the 

-sT  -1 

open  loop  transfer  function  is  simply  e  G(s)  or  Z  G(s). 

The  terms  in  Z  and  s  can  be  separated  and  dealt  with  by  the 
program.  Now  G(s)  can  be  represented  as  N(s)/D(s).  Letting 
G(s)  be  represented  in  this  fashion,  the  closed  loop  transfer 


slbneri  yXne  Lllv  ai  ot->leiscia  baa,  nolifiion  no-®  10 


Y6l9b  omii  «uq  6  ei  mitt  9^"  9EB0  3ci *  riJf1J 


32 


Figure  5 

Control  Loop  With  Time  Delay 


Forward  Loop  Transfer  Function  =  G(s) 

—  c 

Feedback  Loop  Transfer  Function  =  e 
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function  of  the  system  becomes: 


N  ( s ) 
D  ( s ) 


1 


.  "ST  N(s) 

+  e  d7sT 


(35) 


or 


(s) 


N(s) 

D  ( s )  +  Z_1N(s) 


(36) 


For  equation  (36)  the  terms  in  Z  and  s  cannot  be 
separated,  and  because  of  this  the  transient  response  for  this 
type  of  problem  is  out  of  range  of  the  program's  capabilities. 
A  similar  analysis  holds  for  the  case  of  pure  time  delay  in 
the  feed  forward  loop. 
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7,  APPLICATIONS 

The  main  object  of  this  research  was  to  find  the 
types  of  chemical  engineering  control  problems  the  program 
could  be  applied  to  and  to  investigate  a  possible  widening 
of  its  applications. 

There  are  two  main  classes  of  control  systems,  linear 
and  non-linear;  non-linear  being  the  majority,  and  also  the 
more  difficult  to  treat.  Each  main  class  may  be  divided  into 
two  subclasses,  one  is  the  continuous  and  the  second  is  the 
discrete  or  sampled  system.  Except  for  one  case,  the  program 
is  unable  to  cope  with  the  non-linear  type  systems,  however, 
it  deals  very  well  with  most  of  the  linear  types.  The  follow¬ 
ing  will  be  a  summary  of  the  problems  encountered  in  treating 
each  type  of  problem  using  the  Control  Systems  Analysis  pro¬ 
gram  and  an  outline  of  the  steps  used  to  carry  the  problem 
to  solution. 

7.1  Root  Locus  for  Continuous  Linear  Systems  with  or 

Without  Pure  Time  Delay 

The  simplest  control  problem  is  the  single  feedback 
loop,  linear,  continuous  system.  No  problems  were  experienced 
when  analyzing  this  type,  however,  an  outline  of  the  data  which 
must  be  provided  about  the  system  is  presented  so  that  the  ex¬ 
ceptions  can  be  noted  in  following  types  of  problems.  To  enter 
this  type  of  problem  for  stability  analysis,  a  minimum  number 
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of  specifications  must  be  made.  There  must  be  specified: 

a)  the  number  of  feedback  and  forward  loops; 

b)  the  order  of  each  polynomial  in  the  transfer  function 
and  whether  it  is  in  Z  or  s  notation; 

c)  the  scan  area  and  increment  on  both  axis;  and 

d)  the  number  of  different  solutions  wanted  for  the  problem. 
Parameter  changes  may  be  made  for  the  same  problem. 

With  the  above  information  a  root  locus  calculation 
can  be  performed. 

The  next  type  of  system  is  one  in  which  there  is  a 
pure  time  delay  in  either  the  forward  loop,  or  in  the  feed¬ 
back  loops.  In  treating  this  problem,  the  Z-transform  cap- 

•  •  —  s  T 

ability  of  the  program  is  used.  The  factor  e  representing 

the  pure  time  delay  is  entered  into  the  transfer  function  as 
-1  sT 

Z  ,  since  Z  =  e  .  No  Z-transform  calculation  of  the  s-transfer 
function  is  specified  by  the  user,  however,  the  time  delay  is 
entered  as  data  along  with  the  highest  power  of  Z.  This,  of 
course,  is  equal  to  one.  An  example  problem.  Appendix  (A. 9. 2), 
showing  the  effect  of  various  pure  time  delays  on  the  stability 
of  a  control  system  is  shown  along  with  the  input  information 
necessary  for  the  correct  designation  of  the  problem. 

For  the  above  two  types  of  problems  only  the  root 
locus  was  mentioned  as  the  stability  criterion,  however,  be¬ 
cause  of  additions  made  to  the  program,  Nyquist  plots  and  Bode 
plots  are  also  available.  These  are  produced  when  and  only 
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when  specifically  asked  for  by  the  user.  This  involves  set¬ 
ting  the  span  of,  and  the  increment  of,  the  frequency  over 
which  the  user  wishes  to  examine  the  system.  The  user  must 
also  set  a  flag  in  order  to  direct  the  program  to  do  either 
or  both  of  the  calculations  necessary  to  produce  the  two  plots. 
An  example  problem  has  been  included.  Appendix  (A. 9. 3),  show¬ 
ing  the  computed  graphs  for  these  three  types  of  analyses. 

7 . 2  The  Study  of  Sampling  Interval 

The  insertion  of  a  sample  and  hold  device  into  a 
loop  has  a  destabilizing  effect  on  a  control  system.  The  de¬ 
gree  to  which  this  effect  is  felt  by  the  system  depends  a 
great  deal  on  the  sampling  frequency  used.  Generally,  the 
smaller  the  sampling  interval,  the  less  the  destabilizing 
effect  the  sampler  has  on  the  system.  It  would  be  advantage¬ 
ous  to  a  systems  engineer  to  be  able  to  determine  beforehand 
what  effect  sampling  intervals  would  have  on  a  given  system, 
and  to  be  able  to  determine  where  the  limit  of  stability  for 
the  system  occurs.  It  would  then  be  possible  to  specify  the 
sampling  rate,  and  design  on  this  basis  with  some  assurance 
of  overall  system  stability. 

The  program  is  well  suited  to  this  type  of  study. 
Using  this  program,  a  system  can  usually  be  studied  through  a 
root  locus  analysis  to  determine  for  what  gains  instability 
occurs  for  various  sampling  periods.  Study  using  the  transient 
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response  depends  on  the  placing  of  the  sample  and  hold.  A 
different  method  than  the  one  outlined  in  Section  3.2  may 
have  to  be  used  if  the  sampler  is  not  in  the  forward  loop. 

In  this  way  the  program  could  be  used  as  a  valuable  tool  in 
the  conversion  of  a  continuous  control  system  to  a  sampled 
system  in  preparation  for  digital  control. 

7 . 3  Transient  Response  Calculation 

7.3.1  Continuous  Systems:  Change  in  Set-Point 

With  the  addition  of  another  subroutine  "PDIV" , 
Appendix  (A. 5),  it  was  found  that  the  Control  Systems  Analysis 
program  could  be  used  in  the  calculation  of  the  transient  res¬ 
ponse  of  some  systems.  The  calculation  for  continuous  systems 
involves  obtaining  the  Z-transform  of  the  closed  loop  s-domain 
transfer  function  and  inverting  this  Z-transform.  The  inver¬ 
sion  results  in  values  for  the  weighting  sequence  of  the  out¬ 
put.  The  mathematical  formulation  is  outlined  in  Section  3.2. 

A  continuous  system  may  be  represented  as  shown 
in  Figure  6-a,  for  which  the  open  loop  transfer  function  is: 

GoL (s)  =  G1(S)  G2(s) 

and  the  closed  loop  transfer  function  is: 


(37) 
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Using  the  closed  loop  transfer  function  the  block  diagram 
shown  in  Figure  (6-a)  can  be  converted  to  that  shown  in 
Figure  (6-b) .  By  inserting  a  sample  and  hold  before  the  block 
representing  the  system,  the  continuous  system  can  now  be 
transformed  to  a  sampled  system  to  which  Z-transform  theory 
can  be  applied,  Figure  (6-c) . 

As  would  be  expected,  the  sampling  period , T , chosen 
for  the  sampler  which  converts  the  continuous  signal  to  the 
closed  loop  system  to  a  sampled  signal  has  a  great  effect  on 
the  overall  response  of  the  system.  Now,  for  the  sampled  sys¬ 
tem  to  exhibit  the  same  type  of  transient  response  as  the 
continuous  system  the  sampler  should  have  little  effect  on 
the  system.  Since  it  is  the  roots  of  the  system  which  deter¬ 
mine  the  type  and  manner  of  response,  a  measure  of  the  effect 
of  the  sampler  on  the  system  is  the  difference  between  the 
root  loci  of  the  sampled  closed  loop  system  and  the  root  loci 
of  the  continuous  closed  loop  system.  When  the  root  loci  of 
the  sampled  and  continuous  systems  agree  as  to  the  position 
of  their  poles  and  zeros  on  the  s-plane,  the  trajectories  of 
the  roots  in  the  complex  area  of  the  s-plane,  as  well  as  in 
the  values  of  the  relative  gains  along  the  locus,  the  inverse 
of  the  Z-transform  multiplied  by  the  value  of  the  input  signal 
at  that  time  will  be  the  value  of  the  output  function  in  the 
time  domain.  The  inverse  of  the  Z-transform  is  obtained  by 
the  power  series  method. 
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Figure  6 

Block  Diagram  for  Continuous  System 
Transient  Response  Calculation 
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When  using  a  sample  and  hold  device  there  is  a  cer¬ 
tain  amount  of  lag  associated  with  it.  For  instance,  in  most 
cases  a  zero-order  hold  will  give  a  lag  of  half  a  sampling 
period.  To  correct  this  lag  it  is  desirable  to  add  a  ficti¬ 
tious  lead.  Similarly,  the  approximation  of  continuous  func¬ 
tions  using  a  first  order  sample  and  hold  device  is,  in  most 
cases,  improved  by  the  addition  of  a  lead  equal  to  one  samp¬ 
ling  period. 

If  the  transfer  function  of  the  closed  loop  turns 
out  to  have  the  degree  of  the  denominator  only  one  greater 
than  that  of  the  numerator,  a  first  order  hold  is  necessary. 
This  eliminates  the  need  for  the  user  to  define  a  new  Z  to 
account  for  the  transport  lag  produced  by  the  hold.  It  would 
be  necessary  to  define  a  second  Z  if  a  zero-order  hold  were 
used  as  the  sample  and  hold  device.  The  use  of  the  first 
order  hold  is  made  necessary  because  the  program  lacks  the 
capability  to  account  for  a  different  Z,  other  than  the  one 
being  used  for  the  transform,  and  at  the  same  time  calculating 
a  special  Z-transform.  This  is  the  case  when  a  function  is 
treated  which  has  its  numerator  one  less  in  degree  than  its 
denominator.  When  the  order  of  the  denominator  is  two  greater 
than  that  of  the  numerator  either  hold  can  be  used. 

An  example  of  each  type  of  system  is  presented  in 
the  Appendix,  (A. 9. 4.1,  A. 9. 4. 2). 
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7.3.2  Continuous  Systems:  Change  of  Load 


All  the  systems  treated  previously  were  disturbed 
by  a  step  change  in  the  set-point  of  the  controller.  In 
chemical  systems,  disturbances  also  occur  because  of  changes 
in  the  load  variable.  Although  the  root  locus  analysis  is 
the  same  for  both  types  of  disturbance,  since  it  deals  with 
only  the  characteristic  equation,  the  overall  response  is 
different.  This  can  be  attributed  to  the  change  in  the  numer¬ 
ator  of  the  transfer  function  of  the  system. 

A  block  diagram  with  the  load  variable  included  is 
shown  in  Figure  (7-a) .  The  equation  representing  this  system 
is 

G  (  s )  G  ( s ) 

e  (s)  =  - 2 -  9.  (s)  +  - - -  0  (s)  (39) 

1+G  (s)G  (s)  1+G  ( s )  G  (  s ) 

m  p  m  p 


and  when  studying  the  effect  of  load  change,  0^(s)  can  be  set 

to  zero  and  for  G  (s)G  (s)  =  G  T  (s) : 

m  p  oL 


6o(s) 


gl(s) 

1  +  GoL(s) 


eT  (s) 


(40) 


A  block  diagram  representing  equation  (40)  is  shown 
in  Figure  (7-b)  which  can  be  condensed  to  that  of  Figure  (7-c) . 

Thus,  to  arrive  at  the  response  to  a  step  input  or 
any  other  disturbance  in  the  load  variable,  it  is  only  neces¬ 
sary  to  follow  the  steps  previously  outlined  for  the  response 
of  a  continuous  system  to  a  change  in  the  set-point 
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Figure  7 


Block  Diagram  Sequence  for  the  Transient  Response 
for  a  Continuous  System  (Load  Change) 


(a) 


(b) 
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point.  It  should  be  noticed  that  the  denominator  of  the  trans¬ 
fer  function  remains  the  same  for  both  a  load  and  set-point 
disturbance.  This  is  advantageous,  as  it  means  the  denomina¬ 
tor  of  the  closed  loop  transfer  function  must  only  be  factored 
once  to  make  the  transient  response  both  to  load  and  set- 
point  change  available. 


7.3.3  Sampled  Systems:  Set-Point  Change 

The  method  for  obtaining  the  transient  response  of 
these  systems  is  somewhat  simpler  than  that  outlined  previ¬ 
ously  for  the  continuous  system.  For  this  case  it  is  not 
necessary  to  match  root  loci  in  order  to  determine  the  effect 
of  the  sampler,  as  the  sample  and  hold  device  are  integral 
parts  of  the  system.  Because  of  this,  the  Z-transform  of  the 
system  transfer  function  can  be  used  directly  to  arrive  at 
the  transient  response  for  various  types  of  inputs.  The  in¬ 
version  routine  used  in  the  continuous  case  can  also  be  used 
here.  The  standard  rules  for  combining  Z-transforms  apply 
so  that  attention  must  be  paid  to  these  rules  when  arriving 
at  a  transfer  function  to  represent  the  closed  loop  system. 

It  can  be  seen  that  for  some  cases  it  will  be  necessary  to  do 
some  hand  manipulation  of  blocks  in  the  block  diagram  as  the 
program  does  not  have  block  manipulation  capability. 


In  the  procedure  for  obtaining  the  transient  response 
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or  root  locus  for  an  error-sampled  system,  it  was  necessary 
to  multiply,  and/or  add  different  polynomials  together  by 
hand.  This  was  both  tedious  and  inaccurate,  consequently  an 
addition  (Ao3,  A. 4,  A. 6,  A„7)  was  made  to  the  program  to  re¬ 
lieve  the  user  of  this  calculation.  The  addition  can  be  used 
only  for  unity  feedback  error-sampled  systems.  The  main  func¬ 
tion  of  this  addition  is  to  relieve  the  user  of  the  necessity 
of  combining  the  Z-transform  of  the  controller  and  that  of 
the  process  and  to  reduce  the  number  of  runs.  An  example 
problem  is  presented  (A. 9 .6)  to  illustrate  the  method. 

7.3.4  Sampled  Systems:  Load  Change 

The  Calculation  of  the  transient  response  of  a 
system  to  a  load  change  differs  from  that  for  a  set-point 
change.  Figure  (8)  is  a  block  diagram  representing  an  error- 
sampled  control  system  with  unity  feedback.  Given  a  load 
change  and  taking  the  set  point  I(s)  as  zero,  the  equation  for 
the  output  can  be  obtained  from  the  block  diagram  as: 


eTG*(s)  -  D*  (s)  G  *  (s)  0  *  ( s) 

L  L  p  O 


(41) 


Taking  the  Z-transform 


0lGl(Z)  -  D(Z)G  (Z)  eQ(Z) 


(42) 


from  which  the  closed  loop  transfer  function 


' 


45 


0TGT  (z) 

eo  (Z)  =  - — - -  (43) 

CL  1  +  D(Z)G  (Z) 

P 

can  be  formed,  eLGL(Z)  represents  the  Z  [0  (s)  GL  (s)  ]  . 

It  can  be  seen  that  the  input  function  is  not  a 
separate  identity  in  the  numerator,  thus,  it  must  be  Laplace 
transformable.  Because  of  this  the  calculation  of  points  on 
the  transient  response  curve  is  also  different.  These  points 
will  now  be  the  coefficients  of  the  power  series  in  Z  \ 
Section  (3.2),  when  the  power  series  inversion  method  is  ap¬ 
plied,  this  was  not  the  case  for  a  set-point  change.  Thus, 
the  transient  response  calculation  for  a  sampled-data  system 
with  a  load  change  requires  a  variation  in  the  procedure  nor¬ 
mally  used  for  the  transient  response  to  a  set-point  change. 
The  transfer  function  6TGT  (Z) ,  the  numerator  of  equation  (43) , 
must  be  calculated  from  eT  (s)GT  (s)  before  the  overall  transfer 

La  1j 

function,  equation  (43)  can  be  obtained. 

7 . 4  The  Specification  and  Design  of  Digital  Controllers 
7.4.1  Designing  a  Digital  Controller  for  a  System 

An  investigation  into  the  application  of  the  Control 
Systems  Analysis  program  to  the  specification  and  design  of 
digital  controllers  was  made.  The  particular  application  is 
in  the  determination  of  a  digital  controller  to  control  a 
system  previously  controlled  by  a  continuous  controller. 
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Figure  8 

Block  Diagram  for  Load  Change 
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A  continuous  system  is  represented  by  the  block  dia¬ 
gram  in  Figure  (9-a) ,  where  G  (s)  is  the  transfer  function  of 

c 

the  controller  and  G^fs)  is  the  transfer  function  of  the  pro¬ 
cess  . 

Now,  in  converting  to  a  digital  controller,  which 

could  be  a  digital  computer,  modifications  are  made  to  Figure 

(9-a)  to  represent  the  digitally  controlled  system.  This  is 

shown  in  Figure  (9-b) .  D  (Z)  represents  the  digital  controller 

c 

and  H(s)  the  hold  device  which  will  give  a  constant  analogue 
input  to  the  process  between  sampling  instants. 

In  the  design  of  the  digital  controller,  the  root 
locus  in  the  Z-plane  was  found  to  be  the  most  useful  tool. 

The  procedure  followed  was  to  first  obtain  the  root  locus  in 
Z  of  the  system  without  the  controller,  but  the  hold  included. 
It  can  be  determined  where  the  poles  and  zeros  in  Z  of  the  con¬ 
troller  should  be  placed  from  an  examination  of  this  locus. 

This  step  constitutes  one  run  of  the  program.  With  the  con¬ 
troller  poles  and  zeros  decided  upon  a  new  root  locus  is  cal¬ 
culated  either  in  s  or  Z ,  depending  on  the  user's  preference, 
(Appendix  A. 9. 6).  The  advantage  of  one  type  of  root  locus 
versus  the  other  varies  with  the  problem.  The  object  of  this 
step  is  to  obtain  a  knowledge  of  the  overall  loop  gain  for 
which  the  system  will  be  stable.  In  some  cases  one  root  locus 
is  better  than  the  other  for  this.  The  user  now  picks  the 
overall  loop  gain  for  which  he  wants  to  examine  the  transient 
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Figure  9 

Block  Diagram  Sequence  for  Digital 
Controller  Implementation 
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response.  This  is  the  end  of  the  second  run.  The  third  run 
is  to  obtain  the  transient  response  of  the  system  for  the  loop 
gain  decided  upon  previously  in  the  second  step.  The  program 
always  calculates  the  transient  response  of  the  system  for  the 
loop  gain  K  equal  to  one.  This  means  that  when  an  overall  loop 
gain  is  decided  upon  in  step  2,  it  should  then  be  put  in  the 
numerator  of  either  the  process  transfer  function  or  the  con¬ 
troller  transfer  function.  K  would  then  be  included  in  the 
transfer  function. 

Before  the  modification  mentioned  in  Section  7.3.3, 
the  closed  loop  Z-transfer  function  with  the  value  of  gain 
inserted  in  the  numerator  was  calculated  by  the  user.  This 
involved  both  multiplication  and  addition  of  polynomials  and 
was  necessarily  done  with  the  Z-transfer  function,  and  was 
repeated  for  each  different  value  of  K  chosen.  In  order  to 
free  the  user  of  this  time  consuming  operation,  a  modification 
was  added  which  first  takes  the  Z-transform  of  the  process 
transfer  function,  combines  it  with  the  controller  Z-transfer 
function  provided  by  the  user,  computes  the  closed  loop 
transfer  function  and  from  this  calculates  the  transient  res¬ 
ponse.  To  use  this  modification  the  flag  LOOP (Appendix  A. 5) 
must  be  set  equal  to  one.  The  sample  period , total  time,  units 
of  time  and  the  size  of  the  step  in  the  setpoint  must  also  be 
set  (Appendix  A. 5).  The  order  of  each  term,  the  number  of 
terms,  and  the  coefficients  of  these  terms  making  up  the  con¬ 
troller  transfer  function  in  Z,  must  be  read  in  (Appendix  A. 4). 
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An  example  calculation  is  presented  in  Appendix  (A. 9. 6),  with 
the  accompanying  input  data. 

When  using  a  process  control  computer  for  control, 
the  controller  for  the  process  is  in  the  form  of  a  computer 
program,  so  that  once  the  transfer  function  in  Z  has  been 
determined,  it  must  be  programmed  for  computer  use.  An  out¬ 
line  of  the  method  to  do  this  is  presented  in  Kuo (12).  It 
is  sufficient  to  say  that  once  a  Z-transfer  function  has  been 
determined  for  the  controller,  a  digital  program  can  be  built 
directly  to  realize  this  transfer  function. 

It  should  be  noticed  that  once  the  root  locus  for 
the  open  loop  is  determined,  a  loop  gain  is  picked  from  it 
for  which  the  transient  response  of  the  system  is  calculated. 
Various  gains  can  be  picked  from  this  root  locus  of  the  sys¬ 
tem  to  give  different  transient  responses.  Thus,  for  each 
transient  response  calculation  it  is  not  necessary  to  deter¬ 
mine  a  root  locus.  This  results  in  a  considerable  saving  in 
calculation  time. 

7.4.2  Direct  Replacement  of  a  Continuous  Controller 

by  a  Digital  Controller 

The  program  can  be  used  not  only  in  specifying  a 
digital  controller  for  a  system,  but  also  specifying  a  digi¬ 
tal  controller  to  provide  the  same  response  as  when  the  system 
was  under  continuous  control. 
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For  a  continuous  system.  Figure  (9-a) ,  the  open 
loop  transfer  function  can  be  represented  as: 

GoL(s)  =  GCI(s)Gp<S)  (44) 

The  Z-transform  of  which  is 

GoL(Z)  =  GcIGp<Z>  <45> 

For  an  error-sampled  system,  Figure  (9-b) ,  the 
open  loop  transfer  function  is 

goL  (z)  =  Dc(Z)H  G  (Z)  (46) 

Now  for  the  response  of  the  two  systems  to  be  the 
same,  the  root  loci  must  be  the  same,  or. 


1  +  KGC£Gp(Z)  =  1  +  KDc(Z)HGp(Z) 

from  which 


(47) 


Dc  (Z) 


G  tG  (Z) 

C  I  p 


H  G  (Z) 
P 


(48) 


A  transfer  function  for  D  (Z)  is  to  be  determined. 

c 

Dividing  the  terms  in  the  right  hand  side  of  equation  (48)  in¬ 
to  their  respective  numerator  and  denominator  terms,  equation 
(48)  can  then  be  written  as: 


-  ...  -  k  1 
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D  (Z) 
c 


(49) 


Since  there  is  a  one  to  one  correspondence  between 


the  poles  on  the  s  and  Z-planes,  they  can  be  separated  as 
shown  in  the  numerator  of  equation  (49) .  The  zero  order  hold 
used  does  not  in  itself  contribute  any  poles  to  the  denomina¬ 
tor  of  equation  (49) ,  therefore,  equation  (49)  can  be  written 
as  above. 


Now  cancelling  the  identical  poles  in  equation  (49) 


the  final  Z-transform  necessary  to  duplicate  the  continuous 
action  is  then 


(50) 


DC(Z) 


D  t(Z)N„N  (Z) 
cl  H  p 


Thus,  in  order  to  duplicate  a  continuously  controlled 


system  using  a  digital  controller,  the  controller  transfer 
function  must  be  calculated  in  this  manner,  for  it  is  not 
equivalent  to  the  corresponding  Z-transform  of  the  s-transfer 
function  for  the  continuous  controller. 


It  should  be  noticed,  that  in  the  above  derivation 


sampling  period  is  an  independent  factor  and  can  be  set  by  the 
designer.  For  this  sampling  period  the  Z-transform  equation 
representing  the  digital  controller  can  be  derived  which  will 
give  the  same  transient  response  to  a  disturbance  as  the 
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original  continuous  controller.  Thus,  by  this  method,  digital 
controllers  may  be  designed  to  give  a  specified  transient  res¬ 
ponse  for  a  maximum  sampling  interval,  the  response  being 
specified  as  that  of  the  continuous  system  for  the  same  dis¬ 
turbance  . 

For  small  sampling  periods  the  Z-transform  of  the 
continuous  controller,  used  to  represent  the  digital  con¬ 
troller,  gives  approximately  the  same  response  when  it  is 
substituted  directly .  As  the  sampling  periods  become  larger, 
this  direct  substitution  approximation  becomes  less  accurate, 
and  from  the  derivation  shown  previously  in  this  section, 
they  cannot  be  expected  to  be  equivalent.  Since  large  samp¬ 
ling  periods  are  desirable,  the  above  method  seems  particularly 
suited  for  the  design  of  digital  controllers  to  replace  the 
analogue  controllers  in  existing  control  loops.  An  example 
illustrating  this  method  is  presented  in  Appendix  (A. 9. 7). 
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8 .  RESULTS 

1.  It  has  been  shown  by  example,  (Appendix  A. 9. 2),  that  this 
program  can  be  used  for  the  root  locus  stability  analysis 
of  a  continuous  system  with  a  dead-time  in  either,  or 
both  the  forward  or  feedback  loops.  Bode  and  Nyquist 
plots  can  also  be  derived,  (Appendix  A. 9. 3). 

2.  The  root  locus  calculation  was  found  to  be  dependable, 
and  gave  accurate  results  for  all  systems  treated.  These 
systems  were  linear  and  of  three  types: 

a.  The  conventional  algebraic  transfer  function  in  s. 

b.  Pure  time  delay  added  to  type  "a". 

c.  Sampled-data  systems. 

3.  Bode  and  Nyquist  subroutines,  (Appendix  A.l,  A. 2)  were 
added  to  the  main  body  of  the  program  and  the  operation 
of  these  subroutines  was  checked  using  an  example  from 
(18) .  The  plots  corresponded  exactly  with  those  given 
in  the  reference.  The  resulting  plots  are  shown  in 
Appendix  ( A . 9 . 3 ) . 

4o  No  use  can  be  made  of  this  program  when  frequency  res¬ 
ponse  data  alone  is  available.  The  program  can  only 
treat  a  problem  when  its  components  are  represented  as 
transfer  functions  in  s  or  Z . 

5.  Using  the  Z-transform  capability  of  the  program,  the 
transient  response  of  a  number  of  different  types  of 
systems  can  be  obtained. 
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a.  Linear,  algebraic  systems  in  s. 

The  transient  response  for  this  type  of  system  can 
be  calculated  quite  accurately.  The  method  used  is 
outlined  in  sections  (7.3.1)  and  (7.3.2).  Accuracy 
of  the  transient  solution  can  be  very  sensitive  to 
the  agreement  between  root  loci.  To  ensure  the  valid¬ 
ity  of  the  transient  solution  the  root  loci  agreement 
should  be  within  about  one  percent  for  the  gain  at 
each  point  on  the  root  locus,  the  poles  and  zeros 
should  agree  exactly,  and  the  trajectories  of  the 
roots  in  the  complex  plane  should  fall  one  on  the 
other.  Generally,  the  smaller  the  sampling  period, 
the  closer  is  the  agreement  between  the  two  root  loci. 

b.  Linear,  algebraic  systems  in  s  with  dead-time. 

The  method  which  is  used  for  the  calculation  of  the 
transient  response  requires  that  the  denominator  of 
the  closed  loop  transfer  function  be  factorable  to 
first  and  second  order  terms  in  s.  For  dead-time  in 
this  system  this  is  not  possible,  therefore,  the 
transient  response  for  this  type  of  system  cannot  be 
obtained  by  the  method  outlined  in  Section  (7.3.1). 

c.  Sampled-data  system. 

No  matching  of  root  loci  is  required  in  order  to  find 
the  transient  response  of  this  system.  It  is  found  by 
calculating  the  closed  loop  transfer  function  in  Z  and 
inverting  this  transfer  function,  Sections  (7.3.3,  7.3.4). 
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cL  Systems  with  a  digital  controller 

The  transient  response  of  a  system  such  as  this  can 
be  obtainedo  Care  must  be  taken  when  combining  the 
Z-transform  polynomial  of  the  controller  and  that  of 
the  rest  of  the  loop.  If  one  sampling  period  is  some 
multiple  of  the  other,  the  submultiple  sampling 
method (12)  will  most  likely  have  to  be  used.  The 
closed  loop  transfer  function  of  this  would  then  be 
inverted  to  give  the  transient  response. 

6.  Sampling  effect  on  a  system  (7.2) 

Studies  can  be  made  of  this  effect  either  through 
the  use  of  the  root  locus  or  the  calculation  of  the 
transient  response  of  the  system  to  a  disturbance. 

If  the  sampler  is  in  the  forward  loop,  the  transfer 
function 


C(Z) 


G1(Z) 


(51) 


1  +  G1(Z)G2(Z) 


I(Z) 


is  obtained  and  the  transient  response  is  obtained  using 
the  method  outlined  in  Section  (7.3.3). 

If  the  sampler  is  in  the  feedback  loop,  the  Z-trans¬ 
form  is 


C(Z) 


(52) 


1  +  g1(z)g2(z) 


where  the  input  is  not  separable  from  the  transfer 
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function.  For  this  case,  the  method  outlined  in  Section 
(7c 3.4)  must  be  used. 

7.  The  effect  of  a  digital  controller. 

The  program  may  be  used  to  study  the  effect  different 
digital  controllers  would  have  on  a  control  loop.  Design 
of  a  digital  controller  for  a  certain  loop  can  be  done  us¬ 
ing  information  obtained  from  the  root  locus  plot  for  the 
system.  An  outline  of  the  method  is  presented  in  Section 
(7.4.1)  an  example  is  included  in  the  Appendix,  (A. 9. 6). 

8.  Other  uses  of  the  program. 

a.  Demonstrative  for  students 

The  program  is  easy  to  use  and  a  student,  after 
approximately  a  half  hour  explanation  of  the  data  in¬ 
put  method,  should  be  able  to  use  the  program  for  Root 
Locus,  Nyquist  and  Bode  calculations.  At  present,  the 
use  of  the  program  in  calculating  transient  response 
of  continuous  or  sampled-data  systems  requires  more 
than  cursory  knowledge  of  the  program.  The  primary 
requirement  for  this,  and  for  the  design  of  digital 
controllers  is  knowledge  of  the  Z-transform  which  will 
allow  the  user  to  make  subjective  decisions  regarding 
the  calculation  procedure.  These  procedures  are  out¬ 
lined  in  Sections  (7.3)  and  (7.4). 

b.  As  a  research  tool 

It  could  be  used  for  studying  the  effect  that  vari- 
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ations  in  integral  time,  loop  gain,  and  in  the  sampled 
data  case,  sampling  periods  have  on  system  stability. 

The  effect  of  different  control  configurations  could 
also  be  studied, 
c.  "On-Line"  Simulation 

Most  chemical  systems,  being  fairly  slow  in  response 
could  be  simulated  by  this  program  as  long  as  they 
met  the  following  specifications 

(1)  All  equations  in  the  blocks  between  the  output 
and  input  are  Laplace  transformable. 

(2)  There  is  no  dead-time  present  in  the  loop. 

See  Section  (7.2). 

9.  Grid  Size  Specification 

The  increments  for  the  grid  search  used  in  the  root 
locus  calculation  determine  the  information  which  can  be 
obtained.  If  the  grid  is  coarse,  fewer  scan  lines  cross 
a  locus  and  less  points  on  it  are  calculated.  With  a 
coarse  grid  it  is  also  possible  to  miss  a  pole  or  zero. 

The  user  should  always  check  this  possibility.  A  missed 
pole  or  zero  may  be  spotted  easily  by  noting  the  magni¬ 
tude  and  sign  trend  in  the  value  of  the  loop  gain  given 
for  each  point  on  the  locus. 

Grid  size  specification  depends  largely  on  the  purpose 
of  the  calculation.  A  fine  grid  may  give  more  points  on 
the  root  locus  but  it  also  means  more  computer  time.  If 
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it  is  necessary  to  have  all  the  poles  and  zeros  of  the 
system  picked  out,  a  good  criteria  to  follow  is  to  set 
the  grid  boundaries  large  enough  to  include  the  largest 
significant  pole  or  zero.  The  real  axis  grid  increment 
should  then  be  specified  as  approximately  one-third  the 
distance  between  the  two  closest  roots  so  as  to  place  at 
least  two  grid  points  between  them.  The  lower  boundary 
of  the  scan  area  is  set  to  be  the  real  axis  as  the  root 
locus  is  symmetrical  about  the  real  axis. 

10.  The  Accuracy  of  the  Transient  Response  Determination 

For  a  continuous  system,  the  accuracy  is  strongly 
dependent  on  the  agreement  obtained  between  the  two  root 
loci.  Section  (7.3).  Example  problems  are  included  in 
the  Appendix,  (A. 9. 4,  A. 9. 5),  to  illustrate  this  depen¬ 
dence  . 

For  a  digitally  controlled  system,  the  transient  res¬ 
ponse  obtained  may  be  assumed  correct.  The  calculation 
requires  no  matching  of  root  loci, and  hence,  involves 
only  the  inversion  of  the  closed  loop  transfer  function 
in  Z.  Caution  must  be  exercised  to  prevent  the  use  of  a 
sampling  period  greater  than  half  the  bandwidth  of  the 
signal.  For  example,  if  a  sine  input  were  used,  the 
sampling  frequency  should  be  at  least  twice  that  of  the 
input  signal. 


*00-1  9fUt  SB  8XXB  IBS'S  srfi  ©d  o*  d©8  ai  £>B1&  rtBOa 

ai  bebi/Ioni  sis  emaldcysq  slgffiBxa  .(£.  )  n ocios-:  »xr< 
-a*qsb  ailtt  •aM*e«.Ul  <W  ,(8.e.A  ,fr.e.A)  (»iV:»qA 


o  51/  9»  if  •/  c  :o  e  *  - nocJuBD 

'  *  6^:- 


60 


11.  Accuracy  of  the  Root  Locus  Calculation 

A  half  interval  search  is  used  between  scan  points 
for  a  minimum  of  4  iterations.  The  next  interval  is 
then  required  to  be  less  than  0.005  before  the  search 
is  discontinued. 

12.  Program  Requirements 

a.  The  transfer  function  for  each  block  in  the  block 
diagram  be  either  in  s  or  Z  notation. 

b.  If  a  Z-transform  calculation  is  to  be  made,  the  form 
of  the  s  terms  in  the  denominator  must  be  either 

s  +  a 
or 

s2  +  bs  +  c 

c.  The  block  diagram  should  have  no  inner  loops. 
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13.  Pole-Zero  Cancellation 

This  presents  no  special  computational  difficulty. 
For  this  case  the  word  pole  or  zero  or  nothing  at  all 
might  be  printed  out.  It  is  up  to  the  user  to  check 
the  poles  and  zeros  found  by  the  program  and  determine 
the  cause  if  the  number  printed  out  does  not  agree  with 
that  expected  for  the  specific  problem. 
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9.  CONCLUSIONS 

1.  This  program  can  be  used  effectively  for  the: 

a.  Stability  analysis  of  a  continuous  system. 

b.  Stability  analysis  of  a  sampled-data  system. 

c.  Stability  analysis  of  either  system  with  pure  time 
delay. 

d.  Transient  response  determination  for  a  continuous 
system. 

e.  Transient  response  determination  for  a  sampled  or 
digitally  controlled  system. 

f.  Study  of  the  effect  of  sampling  period  in  a  system 
containing  a  sample-hold  device. 

g.  Demonstration  and  comparison  of  Root  Locus,  Nyquist 
and  Bode  Plots  as  stability  analysis  tools. 

h.  "On  Line"  simulation  of  a  process  provided  the 
specification  mentioned  in  the  results,  Section 
(9 . 8-c) ,  are  met. 

2.  The  accuracy  of  the  simulation  of  a  continuous  system  de¬ 
pends  on  the  match  of  the  root  loci  for  the  sampled  and 
continuous  cases c 

3.  The  number  of  points  on  the  root  locus,  but  not  the  ac¬ 
curacy  of  these  points,  is  determined  by  the  grid  in¬ 
crement  set  by  the  user.  The  accuracy  to  which  these 
points  are  determined  is  independent.  For  a  sampled 
system  the  accuracy  of  the  points  on  the  root  locus  may 
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be  affected  by  the  number  of  terms  used  in  the  series  to 
calculate  the  Z-transform. 

4.  The  program  can  be  used  for  the  design  and  specification 
of  a  digital  controller  for  a  system. 

5.  The  program  can  be  used  to  design  a  digital  controller  to 
replace  a  continuous  controller  which  will  give  the  same 
response  as  the  continuous  case. 

10.  RECOMMENDATIONS 

1.  Examine  the  possibility  of  adding  the  Share  Library  Pro¬ 
gram  (2  3)  which  determines  polynomial  roots,  to  the  Control 
System  Analysis  program. 

2.  Further  expansion  of  the  program  so  that  the  closed  loop 
transfer  function  can  be  calculated  by  the  machine  for  non 
unity,  multiple  feedback  systems. 

3.  Subroutine  PDIV  be  expanded  to  include  disturbances  other 
than  the  step  form  currently  used,  such  as  sine  wave, 
point  input,  and  square  wave. 

4.  Subroutine  ADD  be  expanded  to  handle  systems  other  than 
the  unity  feedback  systems  presently  handled. 

Examine  the  possibility  of  approximating  pure  time  delays 
using  polynomial  approximations  for  the  calculation  of 
the  transient  response  of  a  system  containing  a  time  delay 
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11.  NOMENCLATURE 


K 

q 


s 


T 

Z 

D(Z) 


Dc(Z) 
e(0  +  ) 


G  (s) 
G(Z) 


GcL(s) 


OL 

I*(s) 
9  (Z) 


0i  (s) 

0L(S) 

90(s) 

0* . (s) 
1 


closed  loop  gain  of  the  system 

system  number  for  G(s)/s=a+ib  where  a  and  b  are 
numeric 

independent  variable  of  the  Laplace  transform 

sampling  interval/period 

independent  variable  of  the  Z-transform 

denominator  of  a  Z-transform 

Digital  controller  transfer  function  in  Z 

the  value  of  the  sampled  input  signal  at  t  =  0+ 

transfer  function  in  s 

transfer  function  in  Z 

closed  loop  transfer  function 

open  loop  transfer  function 

sampled  input  signal 

Z-transform  of  a  sampled  signal 

general  input  variable 

load  variable 

output  variable 

sampled  general  input  signal 

sampled  output  signal 
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SUBROUTINE  BODE 

Section 

Title 

Page 

A. 1.1 

Subroutine  Listing 

69 

A.  1. 2 

Subroutine  Flow  Diagram 

72 
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SUBRCUT  INE  BODE  I NBOD , OMEG A , DOMEG , OMEGF ) 

c  ^  ^ 

C** _ THIS  SUBROUTINE  CALCULATES  THE  MAGNITUDE  AND  PHASE 

C**  ANGLE  GIVEN  A  TRANSFER  FUNCTION  WHEN  I  IOMEGA  IS  SU6- 
C**  ST  I  TL  TED  FOR  S. 

C»* _ _ _ _ _ ________ _ _ _ _ _ __ 

C**  the  lgg  magnitude  IS  printed  out  under  the  heading 

C**  MAGNITUDE.  IT  IS  ALSO  PLOTTED  ON  THE  Y-AXIS  VS.  LOG 

C** FREQUENCY. 

c** 

C**  THE  ANGLE  THETA  IS  PRINTED  OUT  IN  DEGREES.  IT  IS  ALSO 

C»» _ PLOTTED  IN  DEGREES.  THETA . IS . DEFINED  AS  THE  PHASE 

C**  LAG  CF  THE  SYSTEM. 

C** 

C»*  .  DAIA _ READ  IN  BY  DILL. _ FORMAT  .  .  .  3  E I Q .  5  .  I  5  .  .  I  S  T  CARD 

C**  AFTER  M-VECTOR,  STARTING  IN  COLUMN  6,  ENDING 

C**  WITH  C.C.  40. 

C*  * _ __ _ ! _ J _ j _ „ _ ________ 

C**  IF  THIS  SUBROUTINE  IS  TO  BE  USED,  THERE  MUST  BE 
C**  SPECIFIED 

JO** Q.ME G A  =  S  T  A  R  T I  MG— F-RELULiLCY 

C**  COMEG  -  INCREMENTAL  FREQUENCY 

C**  OMECF  =  FINAL  FREQUENCY 

■£■♦♦_. _ _ NBCJC  =..  1  *  EL  AG  FOR  THIS  PROGRAM.. . 1  SIGNIFIES  ON 

C** 

_ DIMENSION  F  (  7, ICO) _ 

DIMENSION  JO(6C),A(500),M(500),S(2,3b),Z(2,20,6),Q(2), 
1C<  2) «D( 2 ) «B( 

135,3M,CA(36),SLM{35),w(2,10) 

COMMON  JO  ,  A  ,  M  ,  S  ,  Z  , 

1W 

COMMON  G ,C_ t  D_ t B 

REAL  OMEGA, DOMEG , OMEGF 
INTEGER  N BCD 
SAVE  =  CMEGA 

X  I=-C.  IE-1C 
X 2= I  .OE-  10 
J=C 

WR ITEI6, 1  ) 

C** 

£** _ MAJLJJL-  TITLE , _ ___ _ _ 

C** 

NYQ=NBOC 

4  IF  (OMEGA  .GF.  OMEGF )  GO  TO . 5 

N  =  2 
N  H=  1 

. ...QME.GA =OM.EG A.+ D  O  ME  G _ _ _ _ _ _ i 

S  IN,  1  )=OMEGA 
S ( NH  ,  1  )  =0  •  0 
GAIL  SCAN  (NYQ) 


C** 
C ** 


CALCULATE  THE  MAGNITUDE 
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/^^AG  =  SG«T(0;(  1)^C(  I)+Q(2)^Q(2)  ) 

C** 

C**  CALCULATE  THE  PHASE  ANGLE 

C** WHAT  QUADRANT  ARE  WE  DEALING  WITH 

C** 

I  F ( C (  L )  .LT.  XI )  GO  TO  9 

. 1F(Q(1) ,GT.  X2 ) GO  TO  9  _ ______ _ __ _ 

C** 

C**  IS  THE  ANGLE  =  90  OR  270  DEGREES 

C  *  * _ 

I F  (  Q  (  2  )  .  L  T  .  0,0  GO  TO  8 
TFET  A  =90 . 0 

_ G  O  TC  14 _ _ _ _ _ _ 

8  THETA=270.C 
GO  TC  14 

9  IF  (CM2)  .1  T .  0,0  GO  TO  12 

IF(C(  1)  .LT.  0.0)  GO  TO  11 

C** 

C»* _ 1ST  QUADRANT  _ _ 

C*  * 

THETA=A  TAN ( Q ( 2 ) /Q (  1) ) *5  7.295  78 
C-Q  TC  14 

C*  * 

C**  2ND  QUADRANT 

C** _ _ _ _ _ _ _ 

11  THETA=18C.0  +  ATAMQ(2)/Q(  l  )  ) ♦57.295  78 
GO  TC  14 

L 2  IFIU1I  .LT.  O.C)  GO  TO  13 

r** 

C**  4 TH  QUADRANT 

C ** „ .  _ _ _ _ _ _ 

THETA=360.0+ATAN(Q(2)/Q(i) >*57.29578 
GO  TC  14 

C»* _ 

C*  *  3RD  QUADRANT 

C*  * 

_ 13 _ THETA  =  180.0  +  ATAMQ(  2)/Q(  1  )  )  *57. 295  78  _ _ _ 

C** 

C**  PHASE  LAG 

£** _ 

14  THETA=THETA-360.0 
EMEGA=ALOG10 (OMEGA ) 

BMAG=ALQG1Q(AMAG) _ _ _ _ _ _ _ 

C** 

C*  *  WRITE  OUTPUT 

C** _ _ _ 

WRITE(6,2)  OMEGA, BMAG, THETA 

C  *  ^ 

C*»  THE  RESULTS  ARE  WRITTEN  UN  TAPE  3  FOR  THE  AUTOPLOTTER. 

C**  THE  NAG.  PLOT  IS  FIRST.  X-AXIS - LCG  FREQUENCY 

C**  Y-AXIS - LOG  MAGNITUDE  RATIO 

C**  PHASE  ANGLE  PLOT  X-AXIS - LOG  FREQUENCY _ 

C**  Y-AXIS - PHASE  ANGLE 

C**  PLOT  LOG  AMPLITUDE  VS.  LOG  FREQUENCY. 

C**  PLOT  PHASE  ANGLE  VS.  LOG  FREQUENCY. 
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WRITE  (3,3)  BMEGA , BMAG 

J  =  J+  1 

E  (  1 1  J  )  =BME.GA 

E ( 2 »  J )  =  THETA 
GC  TC  A 

25. . „.CJi£.GA=.SAV.£_ . . . . . . . . . . 

W  R I T  E  (  3 » 7 ) 

N  =  J 

_ EO  6  J  -  1  _ 

6  WRITE (3,3)  E(1,J),E(2,J) 

RETURN 

1  . . FORMA T  (  1HL,58X, 1 ItiBGQE  POINTS, //4QX ,9H FREQUENCY, 12X  ,  9 

IhN  AGN ITUDE , 6 
IX,  11FPHASE  ANGLE) 

2  FORMAT  <3?X,3F?C.8) 

3  FORMAT  (2F13.A) 

7  FORMAT  ( 6 H  END  ) 

. . END. _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ 
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Figure  10 

Flow  Diagram  for  Subroutine  Bode 
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SUBRCUTINE  AN Y G  ( NYQ, OMEGA , DOMEG , OMEGF ) 

c** 

£** _ THE  ACTUAL  FREQUENCY  IS  GIVEN  IN  THE  PRINTOUT  UNDER 

C**  FREQUENCY. 

C**  THE  LOG  TO  THE  EASE  10  OF  THE  FREQUENCY  IS  PLOTTED  ON 

£** _ J  HE X-A.Xi.JS..*. . . . . . . . . 

C** 

C**  THIS  PROGRAM  CAN  BE  USED  FOR  A  NYQUIST  PLOT  ANALYSIS. 
C** _ _ _ 

C**  VARIABLES  TQ  BE  SPECIFIED 

C**  ONEGA  =  STARTING  FREQUENCY 

£** _ _ _ _ DQKEG  .«  INCREMENTAL  FREQUENCY . . _ . . . . 

C**  OMEGF  =  FINAL  FREQUENCY 

C**  C.C.  6 — C.C.  35  INCL. 

C*» _ FERMAT _ 3E1C.5  FOR  THE  FIRST  THREE  VARIABLES. _ 

C*  *  NYQ  =  I*  FLAG  FOR  THIS  SUBROUTINE 

C**  FOR  NYQ,  INSERT  A  I  IN  C.C.  5 

C*«  THE  LATA  CARD  CONTAINING  THE  ABOVE  INFORMATION  MUST  BE 
C**  THE  1ST  CARD  AFTER  THE  M-VECTOR  CARD  OR  CARDS. 

C*  $ 

C  I  HENS  ION  JO (60) , A( 500) ,M ( 500) , S(2,35) , Z( 2,20,6) ,Q( 2) , 
IC (2) ,C(2  )  » B( 

135,36  )  t  CA  (36  )  ,  SUM  (35)  ,W(2,10) _ _ 

COMMON  JC  ,  A  ,  M  «  S  ,  Z  ♦ 

1W 

_ CQMMCN _ e _ «_£ _ t _ D _ _ i _ B _ 

REAL  OMEGA, DOMEG, OMEGF 
INTEGER  NYQ 

_ W  RITE . (6,1) _ _ 

A  IF  (CMECA  .GE.  OMEGF)  GO  TO  5 
N  =  2 

_ !S_H=JL _ _ _ 

C** 

C**  INCREMENT  THE  FREQUENCY 
C»* .  . _ 

CMEGA=CMEGA+DOMEG 

C** 

C»*  INITIALIZE  THE  FREQUENCIES. 

C** 

SIN,  I  )  =0M EGA 

SINH,  I)=0.0  . _ . . . . . . . . . . . 

C** 

C**  C ALCLL ATE  COMPLEX  NUMBER  OBTAINED  FOR  I  IOMEGA  SUBSTITU 

_ 1TED  FOR  S. _ 

C*« 

CALL  SCAN  (NYO) 

. g=-g  (JLI _ _ 

c** 

C**  WRITE  OUTPUT 

C** _ _ 

WR  I T  E ( 6  *  2 )  OMEGA, Q( I)  ,Q(2) 

C** 

£**  WRITE  C  N  TAPE  3  FOR  A  U  TCP  LOT  TER.  +  V  E  AND  -VE  IMAGINARY 
1  PARTS  OF  TH 
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C**  COMPLEX  NUMBER  AND  THE  REAL  PART.  AXIS - Y= I M. ,  X=RE.. 

C  ^  ^ 

V*  R  I  T  E  (  3  »  3  )  Q  (  I  )  ,  Q  (  2  )  ,  G 

GO  TO  A 

5  NYG-NYQ+1 
RETURN 

_ 1 _ FORMA I  .  ( It  1 1  4 HMY QU  I  S I  POINT S / / 4 2 X ,  5 H QM E G A .  14X,9HR 

1EAL  PART,?/, 

19HIMAG  PART) 

2  F  G  R  M  A  T  ( a2JUJLEL2II«  8  ) 

3  FORMAT  (3E13.4) 

ENC 
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Figure  11 

Flow  Diagram  for  Subroutine  ANYQ 
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SUBRCUTINE  MULT  ( A , ADEN  1 , NPRO , M ) 

c** 

C»* _ IH.1S  .SUBROUTINE  GUIDES  IME . CALCULATIONS  NECESSARY  TO 

O*  FORM  THE  FORWARD  LOOP  TRANSFER  FUNCTION,  GIVEN  A  CON- 

C**  T ROLLER  IN  Z  AND  THE  PROCESS  TRANSFER  FUNCTION  IN  S 

■C.**. . „ANJQ . LS...C.ALJLE.D. . WHEN  LQQP=1. _ SEE  PDIV. _ _ _ _ _ _ 

C** 

C **  SUBROUTINE  A  Mi  A  L  G  FORMS  THE  NUMERATOR  AND  DENOMINATOR 

_C** _ TERMS  OF  THE  DIGITAL  CONTROLLER  TRANSFER  FUNCTION  INTO 

C**  ONE  POLYNOMIAL  EACH. 

C** 

C  *  * _ A  D  E  N  =  S  T  GRACE . V  E  C  T  O  R  EUR  Q  E  N  0  M INATO  R  C  O  E  F  F  I C  T  E  N  T  S ... . . . 

C**  ANUM=STGRAGE  VECTOR  FOR  NUMERATOR  COEFFICIENTS 

C**  M1=ST GRACE  VECTOR  FOR  THE  NUMBER  OF  TERMS  IN  EACH 

C*» _ 

C**  SUBROUTINE  PMPY  IS  AN  I.B.M.  LIBRARY  SUBROUTINE 
C**  IT  SERVES  TO  MULTIPLY  THE  TRANSFER  FUNCTIONS  OF 

C»* PR Q C ESS  AN C  C ON T ROLLER. THE  NUMERATORS  FIRST  ,  DEN- 

C**  OMINATCRS  SECONC.  SEE  PMPY  FOR  VARIABLE  LIST 
C ** 

_ VARIABLE  LIST _ 

C**  - STORAGE  AREA  FOR  PROCESS  NUMERATOR 

C**  A  DEN  1 - STORAGE  AREA  FOR  PROCESS  DENOMINATOR 

X** _ JYERC.---NUMBER . OF . TERMS .  IN . THE  DENOMINATOR  . 

C**  (v - M-VECTCR 

C** 

f ;  £  *  »  ❖  $  ❖  $  ^  ^  »  *  $  V  £  *  £  Z  V  4  $  £  *  »  £  #  »  »  »  »  £  3  ❖  ❖  »  »  »  »  4  »  »  »  $  $  £  £  £  £  #  $  »  »  »  »  »  # 

C  l MENS  ION  A ( 500 )  , ADEN  1 ( 36 )  , ADEN (20) , ANUM( 20) , Ml ( 10) ,M( 
I5C0)  ,Z ( 36  )  ,  A 

_ 1NUMU20) _ _ _ _ _ _ „ _ _ _ _ _ 

C*  * 

C**  C ALCLL ATE  THE  CONTROLLER  NUMERATOR  AND  DENOMINATOR 

C*» POLYNOMIALS 

C*  * 

CALL  AMALG ( ADEN, ANUM, Ml ) 

. . NP.RQ.=  IA„BS1_NPRQJ _ _ _ _ _ _ __ _ 

CO  1  1  =  1, NPRO 
ANUM1  (  I  )  = A (  I+4CC ) 

1  At  1  +  4CQ)  =  C«Q 

C** 

C**  MULTIPLY  THE  NUMERATORS 

C*» . _ . _ .  . .  . . .  ...  .  . . 

CALL  PMPY  (  Z  ,  IDIMZ,MH  2)  ,  A  NUM  ,  NPRO  ,  ANLM  l  ) 

CO  2  1=1 , IDIMZ 

_ 2 _ A  (  I +4C0 )  =  Z ( I ) _ 

M ( 2)  =  ICIMZ-l 
M  (  l  )  =  i 

C** . . . . 

C**  MULTIPLY  THE  DENOMINATORS 
C** 

_ CALL  P  M  P  Y  (  Z  t  IDIMZ. Ml(4)  .ADEN  ,  NPRO ,  ADEN  1  ) _ _ 

CO  3  1=1 , ICIMZ 
3  ACENI  (  U  =  Z(  I  ) 

£** . . . 
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RESET  R - V  EC  TGR  VALUES 

M  3  )  -  1 

R ( A ) ^  ID IM7-  1 

NPRC=  ID IMZ 

RETURN 

_ JEJSi.C... _ _ _ _ _ _ _ _ _ _ _ _ _ L _ _ _ . . . . . 
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Figure  12 


Flow  Diagram  for  Subroutine  MULT 
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SUBRCUT  INE  AMALG ( ADEN , ANUM , M 1 ) 

DIMENSION  ACEN ( 20 ) , ANUM (20) , M 1 (  L 0 )  , A  1  ( 40 )  , DG 1 ( 40 ) , B 1 ( 3 

1 ,4  ) 


C** 

C** 

£** 

THIS  ROUTINE  IS  CALLED  FROM  SUBROUTINE  MULT.  ITS  PUR¬ 
POSE  IS  TO  MULTIPLY  THE  TERMS  IN  THE  NUMERATOR  AND 

c** 
c** 
c  *  $ 

DENOMINATOR  OF  THE  DIGITAL  CONTROLLER  TRANSFER 

FUNCT ION. 

ONE  TERM  IN  THE  NUMERATOR  AND  ONE  IN  THE  DENOMINATOR 

c  ^  ^ 
o* 
. c  *  * 

IS  FORMED  FROM  THESE. 

IN  ESSENCE,  GIVEN  THE  NUMBER  OF,  AND  THE  DEGREE  OF 

EACH  OF  THE,  POLYNOMIALS,  THIS  ROUTINE  WILL  SUCCtS- 

c** 

c** 

f  *  if  5;  $  * 

SIVELY  MULTIPLY  THEM  TOGETHER  TO  ARRIVE  AT  THE 

POLYNOMIAL  PRODUCT. 

INPUT  DESCRIP T  I.L1N 

c** 

c** 

£** 

THE  INFORMATION  CONCERNING  THE  DIGITAL  CONTROLLER  TO 

BE  INSERTED  IN  THE  LOOP  IS  READ  IN  FROM  THIS  SUBROU¬ 
TINE.  IF  LOOP  IS  NOT  SET  =  1  (SEE  SUBROUTINE  PDIV), 

c*  * 

£**$$$ 

c»* 

THIS  SUBROUTINE  IS  NOT  REACHED  AND  NO  DATA  IS  READ. 

DATA  INPUT 

FOR  LOOP  SEE  PCIV 

c** 

c** 

C** 

THE  NUMBER  OF  POLYNOMIALS  AND  THE  DEGREE  OF  EACH  MUST 

BE  READ  IN  FOR  BOTH  NUMERATOR  AND  DENOMINATOR  USING 

FORMAT  1 C  1 3 .  THE  NUMERATOR  DtSCRIPTION  PRECEDES  THE 

C** 

C*  * 

c*  * 

DENOMINATOR  DESCRIPTION.  IF  THERE  WERE  THREE  POLYNOM¬ 
IALS  IN  THE  NUMERATOR  OF  1ST, 2ND,  AND,  1ST  ORDER  RES¬ 
PECTIVELY  THE  INPUT  INFORMATION  WOULC  APPEAR  AS 

c** 

c** 

c** 

3121.  THE  DENOMINATOR  DESCRIPTION  IS  THE  SAME 

AND  COMES  IMMEDIATLEY  AFTER  THE  NUMERATOR.  ORDERED 

FROM  RIGHT  TO  LEFT. 

c** 

c$$ 

c** 

THE  COEFFICIENTS  ARE  ENCODED  IN  THE  SAME  MANNER  AS 

THOSE  OF  THE  A-VECTOR  FOR  AN  ORDINARY  S-PLANE  ROOT 

LCCUS  DATA  INPLT.  THEY  ARE  ENTERED  IN  ORDER  FROM  RIGHT 

C** 

C  *  *  *  ** 
C** 

TO  LEFT. 

VARIABLE  DEFINITIONS 

ANUM (  I  ) =S TOR AGE  AREA  FOR  POLYNOMIAL  DENOMINATOR 

c** 

c** 

C** 

ANUM ( I )  =  STORAGE  AREA  FOR  POLYNOMIAL  NUMERATOR 

Ml  (  J  )  =  f\  U  M  .  AND  DEN.  DESCRIPTION  VECTOR 

PRCGRAN  L I M I  TAT  IONS 

C*  * 

C*  * 

c  *  * 

PRODUCT  NUMERATOR  OR  DENOMINATOR  DEGREE  MUST  BE  LESS 

THAN  19 

TEN  Ml  VECTOR  ENTRIES  IS  THE  MAXIMUM 

c  ^  ^ 

REAL' (5,  ICC)  (MUI), 1  =  1, 10) 

N  A  =G 

M  H  =  M  1  (  1  ) 

CO  1  1=1, MH  _ _ . 

1 

NA=NA+M1 (  1  +  1 ) 

NTN  =  NA  +  M1  (  1  ) 

N  A  =0 

NB=M  1  (  1  )  +  2 
NB1=NB+Ml  ( N 8 ) - 1 
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H  M  f  I  =  I  1  J 
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0  =  A  2 
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2  NA=NA+M  1  (  I  +  1  ) 

NT=NA+M1 ( NB ) +NTN 

READ (  5, 1C l >  <  A1  (  I  >  ,  I  =  1  ,  NT ) 

NCQU I =0 

JR  =  C 

C*  *  CCUNTER  FOR  COEFFICIENT  DATA 

. - . J.M=.Q 

20  JM=JM+1 

NCCU  T  =  NCOU  T  +  1 

_ NP=M1( JM)-1 _ _ _ 

C**  NUMBER  OF  POLYNOMIALS  TO  BE  MULTIPLIED 

12  J  M  =  J  M  +  1 

_ _ NTl-IAi  SIMM  JM)  ) 

C**  DEGREE  GF  FIRST  POLYNOMIAL 

NT=NT  1  +  1 

£** _ NUMBER  OF  TERMS  IN  FIRST  POLYNOMIAL _ 

CO  15  1=1, NT 
J  R  =  J  R  +  i 

. JL5 _ DjG  1  (  I  )  =  A 1  (  J  R  ) _ 

C**  MCVL  TERMS  CF  FIRST  POLYNOMIAL  FRCM  A1  VECTOR  TO 

C**  TEMPORARY  STORAGE 

_ IF ( NP  .LT.  1)  GO  TO  24 _ 

CO  19  I  B=  I  , NP 

C**  SET  UP  TEST  TO  DETERMINE  IF  MORE  MULTIPLICATION 

_ LS.  .NECESSARY _ ....  . . . 

J  M  =  J  M  1 

NT2=IABS(M1 ( JM) )+l 

£** _ NLMPER  CF  TERMS  IN  SECOND  POLYNOMIAL _ _ _ 

CO  16  1=1, NT2 
JR= JR+ i 

__ _  CG  16  J  =  1 ,  N  T 

16  8  l  (  I, J)=DGl( J)*A1( JR) 

C*  *  MULTIPLY  FIRST  POLYNOMIAL  COEFFICIENTS  BY  NEXT 

£** POLYNOMIAL  COEFFICIENTS 

L  =  NT  1+NT2 

C**  NUMBER  OF  TERMS  IN  COMBINED  PGLVNCMIAL 

_ . . . .C.Q . 17  K  =  1 

17  CG1(K)=0.0 

CO  18  1  =  1, NT2 
CG  L8  J=  1 ,  N  T 

K=  I  +  J-l 

C**  COEFFICIENT  LOCATION 

18  CG  1  (  K  1  =DG  1(K  ItMlI-iJ) 

(!*♦  REPLACE  POLYNOMIAL  TERMS  IN  TEMPORARY  STORAGE  BY 

C*  *  LAST  COMB  I  NEC  POLYNOMIAL  PRODUCT 

_ NT  1=NT1+NT2— 1 _ 

C**  RESET  CCUNTER  TG  DEGREE  CF  COMBINED  POLYNOMIAL 

19  NT=NT 1+ 1 

24 _ I FCNCCU I . •  GEj>__Z  1  GO . I.Q . 22 _ _ — _ _ 

M  1  (  1  )  =  1 
CO  21  J  =  1  ♦  N T 

_ ANUM ( J ) =  D  G 1 ( J  1 _ 

21  CG  1  ( J  )  =0  •  0 
M  1  (2  )=NT 

_ GO  TG  20 _ _ _ _ _ _ _ _ _ _ _ _ _ 
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0=«L 

AT  AO  T  l  HUH  0  9  •  1  MBTAU)3 

J)=*l 
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YPA22 J 33  21  **3 
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ST  +iT/*J 

JAIMO/Y JOS  13/1 « i m • ) 0  v  I  2X-  -•  T  1  m 3  j  0  J Z1  *0 

J, I =  *  S 1  33 
0.3-01)103  W 
ST / , I = I  Bi  33 
T/,J=L  Bi  J3 
I-L+I =* 

VJ  I  TA30J  T/.j  I  31  3  3  330  **0 

1 1  ,  1  )  1  !*■(/)  I  03- (  >» )  1  03  b  I 
1  13  V  [  :  !  IOS  3 

I  .  )0  )>  S  JAI  13/yj  13  03M  / 0 3  I2AJ  **0 

f  —  s  T  t-  ♦  I  I  /  =  I  T  / 

jA  I  -i  ;/J  S  j  i  \'\  3  i  .<  i  'SOJO  F  ‘1  j  T  lOJ  T323S  0*3 

IMTtf»T/  PI 
SS  T  l  S  .30.  T  <  OH)  3T  AS 

I  =  (  I  )  I  vi 
T  %  I =L  IS  33 
dll  )  i=  (  U  3U»V 

) • 0= ( L ) i 03  IS 
1  =  (  S  )  J  m 
OS  3T  00 


-85- 

22  CO  23  J  =  1  t N T 
ADEN ( J ) =CG 1 ( J ) 

23  C G  1  ( J  )  =0 • C 

_ NJJ.3  )  =  1 _ 

V  1  (A  )=NT 

ICO  FORMAT ( 10 13 ) 

LQ.1 . --E.Q-8.MAXi.9m,.5J _ _ _ 

RETURN 

END 
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T  ,1=1  00 

(l  )  !  --  (D/J  IA 
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i  -  (  ,  u 

U  =  (  tM  !  '•» 

i  t  1)1  - ooi 
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APPENDIX  A. 5 
SUBROUTINE  PDIV 


Section 

Title 

Page 

A. 5.1 

Subroutine  Listing 

88 

A. 5. 2 

Subroutine  Flow  Diagram 

91 

Figure  14 
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SUBRCUTINE  PD  I  V ( A , ADEN , N , OT , T , FMT , ST EP , LOOP ) 

CR 

C**  THIS  SUBROUTINE  IS  USED  FOR  THE  INVERSION  OF  A  Z-TRAN- 

C**  S  FORM  .  IT  DIVIDES  THE  NUMERATOR  OF  ThE  Z— TRANSFORM  BY 

CR  ITS  D  ENOM  INA  TOR  . 

C*«  THE  Z-TRANSFORM  TREAT  ED  HERE  IS  THAT . OF  THE  RATIO  OF 

C**  CUTPIT  TO  INPUT v  Thus  THE  EVALUATION  OF  THE  TRAN  S I  EM  T 

C**  RESPONSE  CONSISTS  OF  SUMMING  THE  PRODUCT  OF  THE  Z-TRAN 

CR  -SFQRM  POLYNOMIAL  COEFFICIENTS  AND  ThE  TIME  VALUE  OF 

CR  THE  INPUT  FUNCTION  AT  THE  VARIOUS  SAMPLING  INSTANTS. 
CRR*  INPUT  VARIABLES 

CR  DT  =  INCREMENT  IN  TIME  EQUAL . T L  T HE  S AM PLING 

CR  TIME. 

CR  T  =  MAXIMUM  TIME  FOR  WH I Ch  SOLUTION  IS  REQUIRED 

CR  FMT (  1  )  =  UNITS  OF  TIME  ( M I N . « HRS. » SEC S .  ) _ 

CR  STEP  =  SIZE  OF  STEP  USED  TO  PERTURB  THE  SYSTEM 

CRR*  PROGRAM  LIMITATIONS 

C**  PROGRAM  DIMENSION ING  ALLOWS  FOR . NO  MORE  THAN  450  POINT 

C**  S  ON  THE  TRANSIENT  CURVE. 

C.  *  $  *  *  *  DATA  INPUT 

CR  CT.T.FMT(l)  ARE  ENTERED  UN  THE  FIRST  CARD  AFTER  THE  M- 

C**  VECTOR  CARD.  FORMAT  IS  2E10.5,A6,  START  AT  C.C.  46 

CR  END  AT  C.C.  71. 

C** _  M-VECTOR  CARD.  START  WITH  STEP  IN  C.C.  1  WITH  FORMAT 

CR  E10.5,I5 

CR  LOOP  AND  STEP  ARE  READ  IN  DK1. 

CRR* _ CLOSED  LOOP  ROUTINE _ 

C**  ThE  FLAG  LCCP  IS  USED  TO  ENTER  THIS  ROUTINE.  WHEN 

C**  LCOP  =  1  AND  T  IS  NOT  =  0,  THE  ALGORITHM  FOR  CALCULAT¬ 

OR  ING  TEL  CLOSED  LOOP  TRANSFER  FUNCTION  QF  THE  SYS  I  EM, 
C**  SUBROUTINE  ADD,  IS  ENTERED  AND  USING  THIS  TRANSFER 
C**  FUNCTION  ThE  TRANSIENT  RESPONSE  IS  CALCULATED. 

CR _ 

CR  IF  LCOP  IS  NOT  SET  =  1,  NO  MULTIPLICATION  OF  THE  PRO- 
CR  CESS  TRANSFER  FUNCTION  BY  THE  CONTROLLER  FUNCTION 

.  C  R OCCURS  .  . _ . __ 

C**  IF  T  =  0,  PCIV  IS  NOT  ENTERED. 

C  *  * 

»»»»»**  *»*»»*»*  aRRRRRRRRR^RRRR 
DIMENSION  ANUMI 500) ,ADEN( 36) ,0(450 ) , FMT ( 12),A(500),SUM 
1  ( 45C ) 

IJM  T  E  G  E  .R  N . .  .  _  _ _ _ _ _ _ _ _ _ 

WR  if E(6f 7) 

WRITE(6,8)  C  T  ,  F  M  T  (  1  ) 

_ WR  ITE  (  6,  16  ) _ _ 

N=  I  ABS ( N  ) 

C=T/D  T 

_ IMT  =  C . . . . . . . . __  _ _ _ _ _ _ _ _ _ 

IN  =  N  +  1 
CO  I  1  =  1, N 

_ J  = IN—  I _ 

1  ANUM ( J  )  = A (  1+40C ) 

C** 

C*»  DIVISION  OF  C N  E  POL  Y NO MI  A L  BY  ANOTHER 
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i!>.  U,  ifl  21  .  117UJ32  HOI  Ha  HO  3  3M  IT  MUM  IX  AM  =  T  **3 

(  .A  t  .2  ,  .it  I  H )  J  h  M  7'  2  7  1  '  U  =  If  )  TM  1  **3 

.-T2V2  17  a  ...f  7  M  7  0  1 2U  M3  T  A  30  311  d  -  H  J  7  2  **3 

2H01TATIMIJ  MAHOUHM  *♦*♦*3 
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.1*0  /l  GAO:  j?  A  <7372  3MA  MOO  J  **3 

mnaa  hooj  oj^ojo  ****>3 

j-n  .3/lTUOfl  2M1  TH3  0T  GD2U  2 1  <7033  OAJ1  3HT  **3 

IOJA  JHT  ,0  »  I 

•  ■  A  I 

i..  A  AM  21  T  ,  ,  I2U  G  /  A  03*37/3  21  ,C11A  3'llU3HaU2  **3 

.  j  1  I2MOM23  J  H  ...  11 JOO  i  ■-  y 

*  *3 
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lOIlJitH  H3  J  j  "iH 1  i/!00  3HT  VO  HI  I  731  JO  JJ32HAHT  22J0  **3 

.2HJ000  **3 

•  G  3  H  3  T  H  3  7  OH  cl  VICIH  ,  0  =  T  31  **3 
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1  «J JjTHI 
IN  ,•'  )  3TI 

(in  1  »  T  3  (  8  ,  d  )  iTI  fl  fj 
(01,0)  IT1MW 
(  M  )  2  3  A I  *  H 
1  □  \  T  *  3 
3  =  T  M  1 
] f H=HI 
H  .  1  ” I  I  03 
I  -/lift 

(  )  ,)A+  1)  A=  (  t)  MJHA  I 

♦  f  0 

I  )hA  Vi  JAIMOHYJUM  3 H 0  10  /1I2IVI3  **0 
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c** 

I F ( LCOP  .NE.  1 )  GC  TO  23 

NT=N 

r.** 

C** 

STORE  OPEN  LOOP  VECTOR 

c** 

„.EG  24 . J=  1  »NT . .  . 

24 

ANUM (  J+4CC )  = ADEN ( J ) 

c** 

c** 

C A LCLLATE  CLQSEC  LOOP  TRANSFER  FUNCTION 

C*  * 

CALL  ADD { A  » ADEN  » N ) 

23 

T  =  0  .  C  . .  . . .  ... 

I  M=  l 

C** 

c*» 

CALCULATE  TEE  QUOTIENT 

c** 

2 

1  =  IN 

C(  IM  )  =  ANUM ( 1) AMHiMM _ __ _ _ _ 

c** 

C**  CALCLLATE  POINTS  ON  TRANSIENT  CURVE. 
C  ^  ^ _ 

IF  UM  .LE.  1 )  GC  TO  20 
I CM= IM-  l 

_ SUN  (  LH I y SUM t  LQJN  )  +C(  IM  ).*S.TJE£ 

GC  TC  18 

20  SUM (  1  )  =  C (  1  )  *  STEP 

_ I  DM= I M-  1 _ 

18  WRITEOtlO)  T »  SUM  (  I  M ) 

C*  * 

C**  INCREMENT  TEE  TIME. 

C** 

T  =  T  +  C  T 

C** _ 

C**  CALCULATE  THE  NUMERATOR. 

C*  * 

_ CJQ....J . JgJLrM _ _ _ _ _ 

1  =  1-1 


3  ANUM(J)=ANUM(J)-ADEN( I)*C( IM) 

_ 


C*$ 

c  ** 
c** 

REPLACE  OLD  NUMERATOR  COEFFICIENTS  WITH  THOSE  JUST 

CALCULATED. 

13 

ANUM (  IN  )=C.C 

CO  13  J  =  2  »  IN 

ANUMI  J-l  )  =  A  N  U  Mi  (  J  ) 

C*  * 

c** 

CHECK  FOR  SOLUTION 

LIMIT. 

c** 

I  F ( I M  .GE.  IMT)  GC 

TO  11 

I M= I N  +  1 

GC  TC  2 

1 1 

I F ( LCOP  .NE.  1)  GO 

TO  22 

c** 

LCCP=LCQP  +2 

*o 

£S  )T  J,  l  I  .  H.  H03J  )  J  1 

H=  T  H 

♦  *3 

. 

•  *  3 

i  ,  i  =L  AS  J  I 

{  t  )  /  ,  '  -  l  OOA+L  )  MU  A  AS 

**3 

IK;  I  h.  I  *  HI  J  ::j  jo  h  /  jjojad 

*  *  3 

(  t  ,  .  ,  A  )  A  J  J  A  3 

3.0=1  tS 

i  =  M ! 

T  Ill’ll  J  IT  J  !  A  J  3  3  J  A  3  **•-> 

$*3 

Hl  =  i  S 

(  )  .All)  •  -  (  l  ) 

*  J 

.  3VHU3  TH3I2HA«T  ¥  '  04  3TA JU3JA3  **0 

**3 


S  i  >T  0  3  (I  .  3J.  1)31 

I  -VII  =  M3  1 

M  3  T  Z  *(  i  I  )  3  +  (  v  I  )  U  l  =  I  v  I  )  M  J  l 

til,  3T  .33 

s-lT  <>  *(  I  )  3=  (  I  )  MJ2  os 
i  -  v  I  *  1 0  1 
c  iT'vja ,i  loMUiifl^ 

$  *3 

.  -i  n  3  iT  T  JM3H3HI  **0 

**3 

T3+7=T 

**3 

.flUTA#3MUI'  30T  3TAJU3JO  *.*0 

*  *3 

/  »  I  =  l  £  0  3 
1-1  =  1 

(  I  ) 3* (  1  )  : .  -  (  t  )  i /  = ( l )  A  £ 

*  *3 

,  u  HI  ft  W  ’  TO  131  11.03  HO 1  A  >  i  Ml  K  a  JO  33AJMJH  **D 

.337AJJ3JA3  *  *3 

*  #  3 


3.3=  {  I)  MIKA 
m  » s = t  ei  03 

(  l  )  1  [ .1  /•  A  =  {  J -L  >  MU' A  £  1 

♦  >>3 

.  T  1  UJ  HO  1  TO  JO,'  HO  I  #3J-0  **0 

*  O 

11  OT  3 J  I TM1  .JO.  MI  HI 

i  f Ml =  1 

S  JT  03 

SS  OT  0,  (I  .3/.  4 03  J  )  T  1  J  I 

u  ho:jj  =  9j:)j 

**3 


-90- 

C**  RESTORE  OPEN  LCOP  VECTOR  IN  ORIGINAL  POSITION  IN 
C**  PREPARATION  FOR  ROOT  LOCUS  CALCULATION. 

C  *  * 

011-21  J  =  UNT 

ADEN ( J )  =  A  N  U  M (  J+400 ) 

21  ANUM ( J+4C0 ) =0. C 

2  2  W  R  I  T  E  l  6 , 6  )  (  C  l  J  )  ,  J  s  L  ,  I  M  )  _ _ 

WRITE(6i  15)  STEP 

WRITEI6,  L  7  I  (SUM  I  >,I=1,IM) 

Vv  R  I  T  F  (3 , 1 4  ) 

14  FORM  A  T { 6H  END  ) 

12  FORMAT  ( 1 H 1 ) 

5  FORMAT  (  5E  15-  51 _ _ _ _ _ 

6  FORMAT ( 6E2C. 6/ ) 

7  FORMA T ( IHl ,48X , 34H INSTANT ANEGUS  TIME  DOMAIN  RESPONSE/5 

_ 1AJU2.1L- FOR...  SU _ 

1CCESIV6  MULTIPLES/6  3X  *6H0F  TAU) 

8  FORMAT! 1H3,5X, 34HTIMES  IN  MULTIPLES  CF  TAU.  TAU  =  ,Fb 

_ _  1  .  3  ,  A  6  )  _ 

9  FORMAT! 2E1C. 5, A6  ) 

10  FORMAT { 2E  1  3*  5 ) 

_ L5 _ FORMAT  (  1HJ»5X>5  2HRESPGNS£  CUiUE  PQLNJS  FOR  .  A.  STEP  INPU _ 

IT.  STEP  SIZE 
1  =  »  F  6 . 3  ) 

16  FORMAT! 1H1U  5X«  3 5HZ— TRANSFORM  POLYNOMIAL  COEFFICIENTS/ ) 

17  FORMAT  ( 1HJ»6E20.6) 

RETURN 

E1L0 


Kl 


-  .v  - 

.  1  I  I  <?  H  J  A  I  3  !  >1  1  ^  1  WO  I  3V 

.  !  (  I  1  !  ■ 


S  .1  j  J  /  J  <M  J  ^  )  1  £  J  H 


1  *  I =t  IS  I 

(  j  j.n  K,  )  '  \  =  (  L  ) 

j  .  =  (  -W,_  )  AU  ' 

(  M  I  ,  1  *  L  t  a  )  J  M  t  d  )  3  T  I  «  w 
M  \^  [  ;M  ,  •  )  :l  M  W 
(MI  ,  1=  I  ♦  l  >  (  VI  id)  dll  WW 

{  Ml  ,  M  I  Oi 
(  1  }  )  1  A  MJ  ■  < 

(  I  H  I  )  T  AM*  H 
(  M  lie  )  !  A>.W I  i 
(\d.OS Jd) FAhVO^ 

■  j  3  </»  A  1 ,  A  I  '  1  H  *  6  t  X  8  *  t  l  '1  )  1  A  M  WO  3 

J2  jF»C.l  H£S,X»J 

(  J  MHO,  ’  \<  J  -i  A  I  I  JJV  iVl^^JI 

i  .  .  i  :  3J«U  TJUM  V I  23M1THA6  .<•  . 

1 6n  ,<Mi  as )  i  am*  H 


*o 

♦  O 

*  *  0 


s 

ss 


\  1AM  ']  3  <11 


A I 
SI 

e 

6 


,31. >  A  >  <  21 -ins  3V!  u.1  32i(392MasaUe.uHiiiA«»av. 
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APPENDIX  A. 6 
SUBROUTINE  PMPY 


Title 


Page 


Subroutine  Listing 
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-94— 

SUBROUTINE  PMPY(Z,  IDIMZ,IDIMX,X,IDIMY,Y) 

c*$ 

C**  THIS  SUBROUTINE  MULTIPLIES  TUG  POLYNOMIALS  TOGETHER* 

C***«*  VARIABLE  DEFINITIONS 

C** 

C**  z - VECTOR  OF  RESULTANT  COEFFICIENTS*  ORDERED  FROM 

C**  SMALLEST  TO  LARGEST  POWER. 

C**  IDIMZ - DIMENSION  OF  Z  (CALCULATED) 

C**  x - VECTOR  OF  COEFFICIENTS  FOR  FIRST  POLYNOMIAL* _ 

C**  ORDERED  FROM  SMALLEST  TO  LARGEST  POWER 

C**  I  C  I  MX - DIMENSION  OF  X  (DEGREE  IS  IDIMX-l) 

C* * . . Yr —  VECTOR  OF  C Q E F FJ C  LENT S  FjQ.R_.SE CON D„P OLYNGMIAL, 

C**  ORDER  EC  FROM  SMALLEST  TO  LARGEST  POWER 

C**  ICIMY - DIMENSION  GF  Y  (DEGREE  IS  IDIMY-l) 

C*«  

DIMENSION  Z( 36) ,X(20) , Y(20) 

_ I F ( ICIMX*ICIMY)10,10.2Q _ _ _ 

10  I  D  I M  Z  =0 
GO  TC  50 

__20 _ IDIMZ=IDIMX+IDIMY-1 _ 

DC  30  1  =  l  »  ID  IMZ 
30  Z (  I ) =0*0 

_ _ cd. . A.c . liiiiom- _ _ _ _ _ _ _ 

CC  40  J=1  , ID IMY 
K=  I+J-l 

_AL _ Z  1  K  I  =  X  (  I  )*Y( J)+Z(K) _ 

50  RETURN 
END 


'  ' 

-•  «  ■  ■  : ;  r  :  1  j 


1  ;■  .  :  l 

1031  I  ■  )  1 

.  .  Y  I  !  q  f/)  1  |  •  I  c:  T  ril  ?'  1  H  ’  1  1  '  ' '  ^ 

.  |  j  jvM  I  1  I 

83^-iq  TejJWAJ  0T  T<J3J  J>  *2 

(  -  :  »  I 


**d 
**3 
*o 

*  O 

*  *  3 
**3 
$  $  3 
**3 
**3 

*  *3 

,****♦*♦***3 

( o<; )  v f  i  >s  )  /.t  i  I  ^  l  a  ■ y  l  i 

,  1  .  H  (  V  1  I  I  1  1  M  1 

,.=  VM.il 
0?  3T  33 
i  -  r  i  ' ;  +  ■  '  i  I  l  -  v  l  M 
SM I ai *  1*1  3£  )3 

0 . 0=  (  I  )  S 
M  1 1  «  I  =  I  3  ^  J  3 
r  41 Ul  t  I =t  3 A  03 
I  -L  M  =  # 

(  *  )  5  ♦  (  l  )  y  *  ( l  )  x  =  (  a  )  \ 

14  9  U  T  3  H 
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SUBROUTINE  ADD 


Title 


Page 


Subroutine  Listing 
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SUBRCUT INE  ADD ( A , ADEN L ♦ N2 ) 

£ $  £  £  $  $  $  $  £  $  $  A-  £  #  $  #  #  #  $  A-  $  £  ❖  ❖  $  $ 

C** 

C*«  THIS  SUBROUTINE  IS  CALLED  FROM  PDIV _ 

C**  IT  ACDS  TWO  POLYNOMIALS  TO  GET  THEIR  SUM 

O*  IT  IS  USED  AT  PRESENT  TO  CALCULATE  THE  CLOSED  LOOP 

C»»  TRANSFER  FUNCT ION  Of  A  LOOP  WITH  UM1Y  FEEDBACK 

£**$$$  variable  definition 

C**  A - STORAGE  AREA  FOR  THE  NUMERATOR  OF  THE  FORWARD  LOOP 

L** _ TRANSFER  FUNCTION _ 

C**  ACEN1 STORAGE  VECTOR  FOR  THE  DENOMINATOR 

C**  N 2 THE  DEGREE  OF  THE  LARGEST  POLYNOMIAL.  IT  IS 

C** _ .  ALWAYS.  THE  DEJ^JDMXMATJDB  when  UStD  FOR  The  ABOVE 

C * *  PURPOSE 

C ** 

*  $_$  »»  »$  fcjfr  $  »  ❖  £  $  »»»  £#  «  #  »  fr  »  **  * 

DIMENSION  A( 500 ) , ADEN1 ( 36) 

NT=1ABS(N2)+1 

_ „_JCjD„1 _ I^JLjJNI_ . . . . . . . . . . . . . . 

I  AOENK  I  )=ACEN1(  I  )+A(  1+400) 

RETURN 

END 


( 

**********************  ********* 

V  IQS  •*)  3  v3JJ/  '  <?  I  11  -'  '  Vt'  <  l  11 

' .-,  ^  •.  a  i  )  ot  ,jr  i  ..<>•<■  ">  ;*  | 

V,.  I  iiOJO  m  JIAJUOJAJ  )T  T'3a3*4  IA  032u  21  '' 

"  , 

1  I  I .  I  330 

«OOJ  0*A>mH  3HI  30  .UTA^W  :  ,  ,  ,  , ; ‘ 

1  ;  t  .JA1  «U)  (  )  “  1  :  J  v  AJ  JHT  I  1  3  i-'-'i  l  3HT  S 

,.  )8A  3311  «03  •  '  '  U*U 


**3 
**3 
>*3 

*  *3 
•MO 

:--  *  >  v  *  3 

>  *D 

*  *3 

*  *  3 
**0 

V  *  ;. ; 

**3 

**d 


******  '  .  ,  |* 

(  Ot  )  1  3Q A  ,  (  J  Jc*  )  Jl  2  43SI  I 

i  f  (  ^  '  )  2  J  A  1  =  7  / 

T /  ,  1  =  I  i  33 

(  (>3^  +  1  )  +-(1)1  ^  I  »  -  (  I  )  i  J  -1 

/1UJT3H 
Qtf  3 
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APPENDIX  A. 8 

ORIGINAL  PROGRAM,  SUBROUTINE  MODICIATIONS 

1.  Subroutine  DKl 

a.  A  programming  error  was  corrected  so  that  titles 
for  the  various  problems  could  be  read  in  with  the 
problem.  A  change  of  format  was  necessary. 

b.  Statements  were  added  so  that  output  for  the  com¬ 
puting  system's  autoplotter  could  be  obtained.  This 
made  it  possible  to  get  graphical  output  for  the 
solution. 

c.  New  read  statements  were  added  to  accommodate 
subroutines  BODE,  ANYQ,  PDIV,  and  MULT. 

d.  Modified  so  that  subroutine  POLY  is  called,  when¬ 
ever  a  transient  response  calculation  is  called  for. 

e.  Statements  added  to  obtain  the  Z-plane  root  locus. 

2.  Subroutine  POLY 

a.  Changed  to  print  out  a  second  order  polynomial  under 
both  the  factored  polynomial  and  the  polynomial 
headings . 

b.  The  subroutines  PDIV  and  MULT  are  called  from  this 
subroutine . 

3.  Subroutine  OPUT  -  Modified  to: 

a.  Write  out  Z-plane  root  locus  points  when  called  for 
by  the  user. 


' 


I 


1  .?  ’ 
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b.  Write  on  Tape  3  for  the  autoplotter  when  graphs  for 
the  s  or  Z-plane  root  loci  are  requested.  These 
plots  are  not  available  for  the  same  run. 

4.  Subroutine  SCAN  -  Modified  to: 

a.  Contain  the  branch  to  subroutine  NYQ  which  calculates 
points  for  the  Nyquist  Plot. 
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APPENDIX  A. 9.1 

EXAMPLE  PROBLEMS  AND  SOLUTIONS. 
TRANSFER  FUNCTION  INDEX 


Appendix 

Section 


Transfer  Functions 


A. 9. 2 


GOL(S) 


30 (s+2.0) 


-sT 


s  (s-3 . 0)  (s+10.0) 


A. 9. 3 


GcL  ^ s  ^ 


30 (s+2.0) 


s  (s-3 . 0)  (s+10 . 0) 


A. 9. 4 


GoL(s) 


3.6 (0.2s+l) 

0 . 2s  ( 5 . 0s+l) 


GcL(s) 


0 . 72 ( s+5 . 0) 


s  +  0.92s  +  3.6 


esTH1 (s) GeL (s) 


esT  1-e  ST  9  0 . 72 ( s+5 . 0) 

-  ( - )  (-2 - ) 

T  s  s  +0.92s+3.6 


A. 9. 5 


3.6 (10s+l) 


GOL  <s)  _ 


10s  (s+1)  ( 5 . 0s+l) 


0.072 ( 10s+l) 


GcL(s)  ” 


(s+0.0875) (s  +1.1125S+0.8226) 


esT/2Ho(s)GeL(s) 


=  esT/2  , 


1-e 


-sT 


-)  ( 


0.072 (lOs+1) 


(s+0.0875) (s  +1.1125S+0.8226) 


-) 


8 
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A. 9. 6 


A. 9. 7 


H  (s)  G  (s) 
o  p 


0.1 


s  (s+1) (s+0 .2) 


First  Example 


Dc(Z) 


Z-0.834 


Z(HQ(s)G  (s) 


— |  Zi-U  . 
>  =  ( — 
— '  Z-l. 


■)Z 


0 


1-e 


-sT 


0.1 


(s+1)  (s+0. 2) 


Second  Example 


DC(Z) 


[ 


Z  (Ko(S)Gp( 


(Z-0.95)  (12) 

7 

i  -ST 

1-e 

0.1 

(Z-1.0) 

Li 

s 

(s+1)  (s+0. 2) 

S+0 . 1 

D  (s)  G  (S)  =  - - 

p  s  (s+1)  (s+0 . 2) 


0.1 

H  ( s)  G  (s)  =  - 

°  p  (s+1) (s+0. 2) 


Dc(Z)  Z^0(s)Gp(s)]  = 


0.59557Z2-0.539112Z-0. 00269 

7 

l-e“sT  0.144 

(Z-l)  (0.034278Z  +  0. 023014) 

Lx 

s  (s+1) (s+0 .2) 
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APPENDIX  A. 9. 2 

ROOT  LOCUS  PLOT  FOR  A  SYSTEM 
WITH  PURE  TIME  DELAY 


Purpose : 

To  illustrate  the  procedure  used  when  dealing  with  a 
system  containing  a  pure  time  delay. 

To  illustrate  the  effect  of  pure  time  delay  on  a 
system  through  the  comparison  of  root  loci  of  the  system  for 
various  delays. 


Transfer  Functions: 

Without  Delay 


(s) 


30 (s+2.0) 
s  (s-3 . 0)  (s+10 . 0) 


With  Delay 


(s) 


30  ( s+2 . 0) 
s  (s-3 . 0)  (s+10 . 0) 


e 


-sT 


Outline : 

The  input  data  used  for  this  problem  is  shown  in  order 
of  input  under  the  heading  INPUT  DATA  on  page  105 .  The  input 
format  is  shown  for  each  data  group  in  a  data  set.  The  control 
card  data  and  the  M-vector  data  are  printed  out,  in  expanded 
form,  by  subroutine  DPRINT  for  each  run.  This  provides  the  user 
with  a  check  of  the  input  data.  Further  expansion  of  subroutine 


DPRINT 


.  ■  ~ 


•  ■ 


' 
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could  be  made  so  that  the  transfer  function  or  functions  of 
the  problem  is  also  printed  out.  This  expanded  version  of 
the  control  cards  is  in  itself  a  shortened  description  and 
for  further  information  the  user  is  referred  to  (24) . 

The  root  locus  for  this  system  without  the  time  delay 
is  presented  by  Schilling (18) .  Schilling's  diagram  and  the  one 
obtained  through  the  use  of  the  C.S.A.  program  agree  exactly. 
The  root  locus  for  the  system  without  the  added  delay  is  pre¬ 
sented  so  a  comparison  may  be  made  between  its  root  locus 
and  those  of  the  same  system  with  different  delays  added  to 
it. 

Time  delays  of  1,  0.1  and  0.01  minutes  were  added  to 
the  transfer  function  for  the  comparison.  The  results  ob¬ 
tained  in  the  form  of  root  locus  plots  are  shown  in  Graph  (1)  . 
The  trend  of  the  root  locus  towards  that  of  the  system  con¬ 
taining  no  time  delay  as  the  delay  is  decreased  is  apparent. 
This  type  of  change  with  the  variation  of  the  time  delay  is 
expected  and  serves  to  verify  the  applicability  of  the  pro¬ 
gram  to  this  type  of  problem. 

It  should  be  noticed  that  the  transfer  function  con¬ 
tains  a  positive  pole  which  imparts  some  peculiarities  to  the 
root  locus  diagram  for  this  system,  for  there  is  a  lower 
limit  of  loop  gain,  K  =  1.39,  for  which  the  system  is  stable. 
Further  mention  of  this  will  be  made  in  the  Appendix  Section 


A. 9. 3. 


\o  aoi&iav  bsbnsqxs  siriT  .duo  bedai*q  o»I»  ei  mi  Ido-  q 

-  .  si  babbe  add  uioridiw  •<**  c  :  ‘  ' 
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Graph  1 

Root  Loci  for  Pure  Time  Delay 
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If  a  time  delay  term  is  present  in  the  problem,  four 
more  data  entries  must  be  made  than  were  made  for  the  problem 
without  this  dead-time,  (Appendix  A. 9. 3).  To  make  these  en¬ 
tries  clearer,  the  transfer  function  can  be  converted  to: 

30 ( s+2) 

G  =  - 

s  ( s— 3 )  (s+10)  (Z  +  0) 

sT 

where  Z  =  e 

The  transfer  function  is  now  a  function  of  both  s  and 
Z  and  this  dictates  the  extra  four  data  entries.  These  en¬ 
tries  occur  in  the  A  and  M  vectors  and  the  control  card. 

The  value  of  the  time  delay,  T,  is  entered  in  A (7) 
in  the  A-vector.  The  coefficients  of  the  Z-polynomial  are 
also  entered  in  the  A-vector.  The  method  of  entering  these 
differs  in  no  way  to  that  of  entering  those  pertaining  to  s 
polynomials.  When  the  format  of  the  denominator  of  this 
function  is  encoded  in  the  M-vector,  the  degree  of  the  Z- 
polynomial  is  indicated  by  a  negative  integer,  see  INPUT  DATA 
under  M-Vector,  in  this  case  -1.  When  a  Z-term  is  included 
in  a  transfer  function  it  is  also  necessary  to  specify  the 
highest  power  of  Z  in  the  problem.  This  is  done  on  the  control 
card  by  setting  variable  Jo  (7)  which  corresponds  to  F7  in  the 
expanded  input  data  list. 

The  program  will  accommodate  up  to  six  different  time 
delays,  thus,  the  number  of  time  delays  which  it  can  expect 
must  be  set  in  JO (8)  or  F8  in  the  expanded  input  data  list. 

For  this  case  JO (8)  is  equal  to  one. 


.  •  £>•  fib  .  ‘ 


(0  +  3)  (01+8)  (£-8)8 


n  1  1o  ip£b  3ri^  ^icboov-M  ab?  c  baboons  i  noiioni/5 

ii  >  :.i a  o  c  >:t*bon  jt  :  •  jt  x.  tc  oriV 

.  f;  i  ir.fa  buqn .  »l)!i|axL  s;'j  i  83  :to  (c')Cl-  ru  £*i  iEJjm 
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INPUT  DATA 


CCNTRCL  CARD  CATA  (FORMAT  2413) 


2C  8  ICC  l  I  100-20 

TITLE _ 

EXAMPLE  SCHILLING  PG.  214.  DEAC-T I  ME  =  1 


A-VECTCR  ENTRIES  (COEFFICIENT  DATA.  FORMAT  9E8.5) 

-1  .0  C .  5 _ _ CL._0 _ 20.0  4.C _ l  .0  _ 

2 . C  l.C  3C.0  IC.C  i.C  -3.0  L.O  0.0 

1.0  C.C  1.0 

M-VECTCR  ENTRIES  (FORMAT  2413) 

_ 2 _ 10  4  11 _ J  -1 _ _ 

OTHER  FLAGS  AND  INPUT  DATA.  THESE  ARE*  N YQ , OMEGA , DOMEG, 
CMEGF*NBOD,NZRT  »DT»T*FMT(  1) , STEP, LOOP. 

THESE  VARIABLES  ARE  ENTERED  AS  ABOVE  ACCORDING  TO  THE 
FOLLOWING  FORMAT  IF  NOT  APPLICABLE  TO  THE  PROBLEM 

_ NOTHING  NEEC  BE  ENTERED.  THE  CARDS  MUST  STILL  BE 

INCLUDED. 

1ST  CARC  15, 3E  10. 5,21 5,2E10.5, A6 

2ND  CARD  E1C.5,  15 

FOR  THIS  PROBLEM  THE  SUPPLEMENTARY  DATA  ENTRIES  ARE, 

NIL 

DATA  END _ 


MAO  Ijq/I 


{  TAM*H>  AT  A3  3HAJ  J  J « T  Vi  1 D 

0  S-  0  0  I  1  1  '  }  1  y 

3  J  T  II 


'  .  -  n-  ..'■  • 

«  .  3<  i  .AT  AO  1 


,0 


0.  1 

.  J 


J  .  ** 

J.£- 


.  0  s 

0.  J 


0.0 

3.01 


0  .  d 

o.o£ 

o.  i 


2.0 
3.1 
J  .0 


0.  I  - 

o.s 

0.1 
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FI  NC.  CF  A-VECTOR  TERRS 

N  *  20 

F 2  NO.  CF  M -VECTOR  TERRS 

N  =  8 

F  3  NO.  CF  RINS  TC  BE  RADE. 

N  *  1 

FA  ONE  CF  FOUR  OPTIONS  SPECIFIED 

-IF  Z-TRANSFCRR!  TC  BE  CORIPUTED  N  =  0 

-TWO  SAMPLER  SYSTEM.  N  =  N 

-ALL  TRANSFER  FUNCTIONS  IN  ONE  FORWARD 

LCOP  WITH  UNITY  FEEDBACK.  N-0 
-OTHERWISE  N»N 

N  =  0 

F 5  ONE  CF  THREE  OPTIONS  SPECIFIED. 

-FCR  RCOT  L-CCUS  CF  CONTINUOUS  SYSTEM  OR 

FOR  A  SYSTER  IN  Z-FORM.  N  =  (J 
-FCR  RCOT  LOCI  POINTS  OF  A  CNF -SAMPLER 
SYSTER  BUT  NO  Z-FORM  N=l 

-FCR  Z-TRANSFCRM  COMPUTED  OR  ROOT  LOCI 
POINTS  FOR  TWO-SAMPLER  SYSTEM  N--UCH...) 

2 

ii 

I 

o 

F 6  HIGHEST  POWER  CF  S 

N  =  l 

F  7  HIGHEST  POWER  CF  Z 

N  =  l 

F  8  NUMBER  OF  VALUES  ASSIGNED  TO  1 

N  =  l 

F 9  RCOT  LOCI.  -,  +  , OK  EC  TH  FEEDBACK  OPTION  (t,3»0)  N  =  I 


FIG  SCAN  CONTROL  (N  =  0, 1,-1)  V  +  H,H,V 

N  =  V0 

F 1 1  ROC  IF  I  EC  Z-FORR  OPTION  (N=0*l) 

N  =  -0 

M2  REPORT  HEADING  OPTION  (N-+,-2) 

UNUSUAL  Z-FORM  OPTION  (N=i,2) 

N  =  -2 

- 
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F  13 

LOCI  OPT  I  UN  t  USU AL  L Y  N  =  0 

N 

*  0 

F  14 

TERMS  IN  SERIES  FOR  G*(S)#N=0 

GIVES  19  TERMS 

N 

a  -0 

F  16 

IF  Z -TRANSFORM  TO  BE  COMPUTED 

N.NE.O 

N  =  C  EGR  E  E  OF  RESULTING  Z-FORM  DENOMINATOR 

4N  RCCT  LOCUS  POINTS 

N 

=  “0 

F  17 

B-MATRIX  YES  N=i 

N 

*  -0 

F  19 

REAL  PART  N  =  - 1  OK 

N 

=  -0 

M-VECTOR  DATA 


FEEDBACK  LOOP  NO.  1 

SPECIFICATION  OF  THE  FEEDBACK  LOOP  COMPONENTS 
AS  TO  CEGREE  OF  NUMERATOR  AND  DENOMINATOR 


UNITY  FEECBACK 


FORWARD  LOOP  NO.  I 


NO.  CF  TERMS  IN  NUMERATOR  N  =  2 

DEGREES  OF  THESE  TERMS  -VE  INDICATES  Z-FORM 
'  1  C 


NO.  CF  TERMS  IN  DENOMINATOR  N  =  4 

DEGREES  OF  THESE  TERMS  -VE  INDICATES  Z-FORM 
1  I  1  -I 


_1  * 


,,,3Z  M' 


T<US  JA3» 

/ 


i 


m  Vi  ril  I  A  "  3  v  ’  »  w 
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APPENDIX  A. 9. 3 

ROOT  LOCUS ,  BODE ,  NYQUIST  DIAGRAMS 


Purpose : 

To  illustrate  the  modified  Control  Systems  Analysis 
program's  capability  in  producing  these  plots. 

Transfer  Function: 


30 (s+2) 
s (s-3) (s+10) 


Outline : 

The  Root  Locus,  Bode,  and  Nyquist  diagrams  for  this 
system  are  presented  in  Schilling (18) .  The  operation  of  the 
program  was  checked  by  comparing  the  computed  results  and  those 
given  in  the  reference.  Graphs  (2, 3, 4, 5)  agree  exactly  with 
those  given  in  (18) . 

For  instruction  on  how  to  obtain  one,  or  all  three  of 
these  plots,  for  a  system  such  as  the  one  above,  refer  to  the 
listings  of  the  subroutine  or  subroutines  concerned,  (A.1,A.2). 

These  plots  can  be  obtained  for  linear  continuous  sys¬ 
tems,  or  for  systems  with  a  pure  time  delay. 

The  Input  Data  for  this  problem  is  listed  on  pages  115 


to  117. 


A  problem  with  a  pure  time  delay  element  is  not  in¬ 
cluded.  The  problem  format  would  be  the  same  as  that  of  the 
problem’  presented  in  Appendix  (A. 9. 2)  except  for  the  additional 
flags  and  parameters  defined  and  designated  in  the  listings 
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of  the  programs  concerned  such  as  NYQ  and  BODE. 

Graph  2  represents  the  root  locus  of  the  system.  Two 
sets  of  curves  are  apparent,  one  composed  of  black  dots,  the 
other  a  solid  line»  The  solid  curve  is  the  root  locus  dia¬ 
gram  and  the  dotted  curve  is  a  plot  of  the  negative  value  of 
the  loop  gain  for  each  point  satisfying  the  characteristic 
equation.  It  can  be  seen,  how  for  poles,  the  gain  curve  goes 
to  zero  and  for  zeros  it  goes  off  to  infinity,  and  for  each 
branch  of  the  locus  there  is  a  corresponding  gain  branch. 

Since  the  value  of  the  gain  plotted  is  the  negative 
value  of  the  true  gain,  points  on  the  root  locus  are  indicated 
by  negative  values  of  gain.  Points  are  also  calculated  for 
values  of  true  negative  gain,  as  long  as  the  equation 
G(s)  =  -1/K  is  satisfied.  These  are  plotted  and  written  out 
as  positive  gain  values. 

The  root  locus  is  defined  as  the  paths  of  the  roots 
of  the  characteristic  equation  as  the  parameter  K  is  varied 
from  zero  to  infinity.  This  excludes  all  negative  values  of 

K,  therefore,  these  points,  though  satisfying  the  characteris- 

* 

tic  equation,  must  be  excluded  from  the  root  locus.  The  user 
can  do  this  quite  easily  either  through  the  printed  output  or 
by  the  use  of  the  gain  curve  plotted  with  the  root  locus.  When 
the  gain  curve  has  negative  values,  the  points  corresponding 


to  these  values  are  on  the  true  root  locus. 
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Graph  2 

Root  Locus  Diagram 
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I 


0  '£• 


Ill 


The  Nyquist  plot  for  the  system  is  presented  as  Graph 
3.  The  stability  of  the  system  for  which  the  Nyquist  plot  is 
calculated  is  indicated,  for  various  values  of  loop  gain,  by 
the  number  of  encirclements  about  the  point  -1/K.  These 
clockwise  encirclements  must  equal  the  difference  between  the 
number  of  poles  and  zeros  in  the  right  half  of  the  s-plane. 
From  the  Nyquist  diagram,  (Graph  3) ,  this  condition  is  met 
when  1/K  becomes  greater  than  -0.70  or  K  becomes  greater 
than  1.43  which  is  corroborated  by  the  root  locus  plot, 

(Graph  2) ,  of  the  same  system. 

The  Bode  Plots,  Graphs  (4)  and  (5) ,  were  also  calcu¬ 
lated  for  the  system.  The  system  itself  is  open  loop  un¬ 
stable  because  of  the  pole  in  the  right  half  of  the  s-plane. 
This  positive  pole  makes  it  a  non-minimum  phase  system  for 
which  caution  must  be  exercised  in  applying  the  Bode  criteria. 
For  such  a  system  the  Nyquist  plot  should  be  referred  to,  to 
establish  the  regions  of  stability  before  utilizing  the  Bode 


method. 
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Graph  3 
Nyquist  Plot 
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-0.3 


• 

Graph  4 

Bode  Plot  Magnitude 
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Graph  5 

Bode  Plot  Phase  Angle 
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Oi  n  rd  CO  <D  <c  C  D>rH  0 


200.0. 
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INPUT  DATA 


CONTROL  CARD  DATA  (FORMAT  2413) 


- 

03 

1  C  C  1  C 

C  1 

0  0-2 

0  0 

TITLE 

NYCUIST 

» 8CCE  , AND  ROCT 

LOCUS 

FOR  SCHILLING  EXAMPLE  PG. 

214 

A- V  EC  TC  R 

ENTRIES  (COEFFICIENT 

DATA. 

FORMAT  9E8.5) 

-C.5 

l.C  4.0 

O 

♦ 

o 

12.0 

o 

• 

o 

o 

. 

o 

o 

• 

o 

2  •  C 

l.C  3  C  .  0 

1C.C 

1.0 

-3.0  1.0 

c 

# 

o 

1.0 


M-VECTCR  ENTRIES  (FORMAT  2413) 

2 

l  C  3  1  1  l 

OTHER 

FLAGS  AND  INPUT  DATA.  THESE 
CMEGF,NB0C»NZRT  » DT  ,  T»  FMT (  1  )  , 
THESE  VARIABLES  ARE  ENTERED 

ARE,  NVQ, OMEGA, DGMEG, 

STEP, LOOP. 

AS  ABOVE  ACCORDING  TO  THE 

FOLLOWING  FORMAT  IF  NOT  APPLICABLE  TC  THE  PROBLEM 
NOTE  ING  NEED  BE  ENTERED.  THE  CARDS  MUST  STILL  BE 

INCLUDED. _ 

1ST  CARD  I5,3E10.5,2I5,2E10.5,A6 

2NC  CARD  E1C.5,  15 


FOR 

THIS 

PROBLEM 

THE  SUPPLEMENTARY 

CATA  ENTRIES 

ARE, 

NYC 

=  1, 

OMEGA  = 

O.C,  DOMEG  =  i.O, 

OMEGF - 1 00 • 0 , 

NBOD  =  1 

DATA 

ENC 

HMMaHSHMMMMMI 


AT  I J  HU  1 


((  !*>$  !AMn;n)  AH'  3  3HA3  J3AT/33 

00  S-  00103103  *v  81 

'  3JTIT 


A  i  S 
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u.  - 
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U.  I 


0.  Oi 

o.  r- 


3. SI 
o.  i 


o.o 

o  .OJ 


o.a 

.Of 


3.1 

0.1 


c'  .  0- 
D.S 
0.1 
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FI  NO.  CF  A-VECTCR  VERMS  N  =  16 


F 2  NO.  CF  M-VECTCR  TERMS  N  =  7 

F  3  NO.  CF  RINS  TC  BE  FADE.  N  =  l 

FA  CNE  CF  FCUR  OPTIONS  SPECIFIED 

-IF  Z-TRANSFORM  TC  BE  COMPUTED  N  =  0 
—  TWO  SAMPLER  SYSTEM.  N=N 

-ALL  TRANSFER  FUNCTIONS  IN  ONE  FORWARD  N  -  0 

LOOP  WITH  UNITY  FEEDBACK.  N=0 
-OTHERWISE  N=h  _ 

F  5  ONE  CF  THREE  OPTIONS  SPECIFIED. 

-FCR  RCCT  LOCUS  CF  CONTINUOUS  SYSTEM  OR 
FCR  A  SYSTEM  IN  Z-FORM.  N=0 
-FCR  RCCT  LOCI  POINTS  OF  A  ONE-SAMPLER  N  =  -0 

SYSTEM  BUT  NO  Z-FORM  N=1 
-FCR  Z-TRANSFORM  COMPUTED  OR  ROOT  LOCI 
PCINTS  FOR  TWO-SAMPLER  SYSTEM  N=-(10+...) 

F 6  HIGHEST  POWER  OF  S  N  =  1 

F 7  HIGHEST  POWER  OF  Z  N  =  -0 

F 8  NUMBER  CF  VALUES  ASSIGNED  TO  T  N  =  -0 


F 9  ROOT  LOCI.  -,+,0R  BOTH  FEEDBACK  OPTION  (1,3,0)  N  =  I 
F10  SCAN  CONTROL  (  N  =  0 , 1 ,- l  )  V+H ,  H  ,  V  _  _ _ N  -  0 


Fll  MODIFIED  Z-FORF  OPTION  (N=0,1)  N  =  0 

FI?  REPORT  HEADING  OPTION  (N=+,-2)  N_*  -2 

UNUSUAL  Z-FORM  OPTION  (N=l,2) 


. 

1  11 


' 
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F  1 3  LOCI  OPT  ION t  USUALLY  N«Q _  _ _  N  =  -0 

F14  TERNS  IN  SERIES  FOR  G*(S)tN=0  GIVES  19  TERMS  N  =  -0 

F 16  IF  Z-TRANSFORM  TO  BE  COMPUTED  N.NE.O_ 

N  =  C  EGR  E  E  OF  RESULTING  Z-FORM  DENOMINATOR  N  -  -0 

4N  RCOT  LOCUS  POINTS 


F  1  7 

B-NATRIX 

YES 

N  =  1 

N  =  -0 

F  19 

REAL  PART 

N  =  - 

-10K 

N  =  -0 

M-VECTOR  DATA 


FEEDBACK  LCCP  NO.  1 

SPECIF  ICAT  ICN  OF  THE  FEEDBACK  LOOP  COMPONENTS 
AS  TO  DEGREE  OF  NUMERATOR  AND  DENOMINATOR 


UNITY  FEEDBACK 


FORWARD  L COP  NO.  1 


NO.  CF  TERMS 
DEGREES  CF 
1  0 


IN  NUMERATOR  N  = 
TFESE  TERMS  -V E 


2 

INDICATES  Z-FORM 


NO.  CF  TERMS  IN  DENOMINATOR  N  =  3 

DEGREES  OF  TFESE  TERMS  -VE  INDICATES  Z-FORM 

1  1  l 


p  =  *4  V  J  JAU2U  »  HD  IJf  10  :  /  ' 


0—  *  H 


MM  T  AH  I  '40)430  OtfA  '  ■  t^^JA  -  33  01  2A 
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APPENDIX  A. 9. 4 

TRANSIENT  RESPONSE  CALCULATION  FOR  A 
A  CONTINUOUS  SYSTEM  (1st  order  hold) 


Purpose : 

To  illustrate  the  use  of  the  Z-transform  capabilities 
of  the  program  to  predict  the  transient  response  of  continu¬ 
ous  systems  to  step  disturbances  in  the  set  point. 


Problem  Source: 

The  system  was  designed  to  represent  a  typical  control 
problem.  A  solution  for  the  response  of  this  system  was  first 
obtained  using  an  analogue  computer.  This  was  used  to  check 
the  solution  for  the  same  problem  given  by  the  digital  com¬ 
puter. 


Transfer  Functions: 
Open  Loop 


(s) 


Closed  Loop 


3.6 (0.2s+l) 

0 . 2s (5 . Os+1) 


(s) 


3.6  (0.2s+l) 
s2  +  0.92s  +  3.6 


0.72 (s+5.0) 


s2  +  0.92s  +  3.6 


-uniinoo  So  b8noqeo^  3n&ien&i3  9rfi  Jbibs'iq  c  me  g  '  !o 

,.  ' 


119 


Outline : 

The  analogue  solution  is  calculated  for  a  step  input 
of  1.8.  The  digital  solution  is  for  the  same  step. 

In  order  to  obtain  the  digital  solution  a  certain  pro¬ 
cedure  must  be  followed  and  certain  criteria  must  be  met.  An 
outline  of  the  procedure  follows;  calculations  are  presented 
in  order  to  clarify  each  step. 

1.  Obtain  the  closed  loop  transfer  function  of  the  system. 
This  is  a  user  calculation. 

2.  Using  the  Control  Systems  Analysis  program  determine  the 
root  locus  for  this  transfer  function.  Graph  (6) . 

3.  Decide  on  the  type  of  hold  to  be  used.  In  this  case  a 

s  t 

1  order  hold  is  necessary  because  the  order  of  the 
denominator  is  only  one  greater  than  the  numerator.  Com¬ 
bine  the  hold  transfer  function  and  that  of  the  closed 
loop  process. 

4.  With  this  hold,  pick  a  sampling  period  T  and  calculate 
the  root  locus  of  the  modified  transfer  function,  (Graph 
7) ,  in  the  s-plane  using  the  C.S.A.  program.  If  this 

root  locus  does  not  coincide  with  that  of  the  closed 

function 

loop  transferAwithout  the  hold,  a  new  T  (smaller)  must 
be  picked  and  the  calculation  repeated. 

sT  0.72(s+5.0) 

=  e  H1(s) 


G  ( s ) 


s2  +  0.92s  +  3.6 
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5. 


6. 


The  first  order  hold  is  expressed  as 

-sT 


1  -  e 


H,  (s)  =  -  (• 

T 


■) 


usually 

This  type  of  holdAgives  a  lag  of  one  sampling  period  to 

the  system.  To  offset  this  inherent  lag  it  is  necessary 

sT 

to  add  the  same  amount  of  lead  in  the  form  of  e 
When  the  root  loci  of  the  sampled  system  and  the  continuous 
system  coincide  in  terms  of  root  trajectory,  pole-zero 
position  and  gain,  the  continuous  system  is  simulated 
exactly  by  the  sampled  system.  This  leads  to  the  next 
step  of  evaluating  the  transient  response.  Graph  (8) . 

Set  the  transient  solution  limits,  the  time  increment  and 
step  magnitude.  If  the  root  locus  is  wanted,  keep  JO (16) 
positive,  if  it  is  not,  change  JO (16)  to  a  negative  number. 
JO (16)  is  the  sixteenth  variable  in  the  Control  Card 


Vector. 

7.  The  points  on  the  transient  response  curve  are  calculated 
and  printed  out  on  paper  as  well  as  on  cards  for  the 
Autoplotter,  (Graph  8)  . 

The  analogue  computer  solution  is  shown  along  with  the 
digital  solution,  Graph  (8) ,  for  ease  of  comparison.  It  can 
be  seen  that  for  this  case  where  the  root  loci  have  been 
matched  almost  exactly,  T  =  0.01,  the  agreement  between  the 
two  solutions  is  excellent. 


. 


' 

b-MO  Xo^xioO  3riJ  ni  •Id*i*ev  arid  ei  I 


Analogue  and  Digital 
Transient  Response  Comparison 
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Time  (min 


. 
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The  input  data  for  each  separate  step  is  shown  for 
problem  illustration.  The  data  for  the  sampled  system  and 
the  corresponding  continuous  system,  pages  125  to  130  differ 
for  two  reasons,  a  Z-transform  must  be  obtained,  and  the  sample 
and  hold  plus  the  lead  term  must  be  added  to  the  transfer  func¬ 
tion.  These  additions  convert 


0.72  (s+5) 

s)  =  — 2 

s  +  0.92s  +3.6 

sT  ,  -sT  n 

e  1-e  9  0.72(s+5) 

—  ( - r(-2 - ) 

T  s  s+0 .92s  +  3.6 

(Z-l)  (Z-l)  (0.72)  (s+5) 

.01)  (Z+0)  ( s+0)  (s+0)  ( s2+0 . 92s+3 . 6 ) 

The  denominator  which  is  submitted  for  Z-transformation 
contains  only  the  poles  of  the  system,  no  others.  In  this 
case  the  denominator  submitted  is 

s2  +  0.92s  +  3.6 

and  is  the  first  polynomial  to  be  encoded  in  the  A-vector. 

The  equation  containing  the  hold  and  time  advance  is  then 
entered  in  the  normal  fashion.  For  more  information  regard¬ 
ing  this  see  (24) . 


to 


GcL(s)  “  Z 


or 


GcL(s,Z)  = 


(( 


io~*  tb  OCX  oi  5SI  ,msd3\’8  euoua  inoo  pnibaoqasiiOD  9dd 


-oii:  i  an&?J  ertt  oJ  b9bbe  ed  .-taum  n  1  )  sj  oitS  ax/Iq  blori  baa 


' 


s.t  +  sse.o  +  b 
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INPUT  DATA 


CONTROL  CAR  C  CATA  (FGRMAT  2  A  I 3 ) 

15 

5  1C02C0100-200 

0  0 

TITLE 

R  CO  T 

LOCUS  OF  THE  CLOSED  LOOP  (CONTINUOUS) 

A-VECTCR  ENTRIES  (COEFFICIENT  DATA.  FORMAT 

9E8.5 ) 

-C.  1 

5  •  C 

l.C  1.0  C.O  6  •  C  5 .  C 

1.C  C • 72  3.6  C . 92  l.C 

0.0  0.0 

M-VECTCR  ENTRIES  (FORMAT  2413) 

2 

1  C  1  2 

OTHER 

FLAGS  ANC  INPUT  DATA.  THESE  ARE,  Nr 

Q,  OMEGA, DOMEG, 

CMcGF  »NBOD,NZRT  »DT ,T,FMT(  1  ), STEP, LOOP 
THESE  VARIABLES  ARE  ENTERED  AS  ABOVE 
FOLLOWING  FORMAT  IF  NOT  APPLICABLE  TC 

0 

ACCORDING  TG  THE 
THE  PROBLEM 

NOTHING  NE EC  BE  ENTERED.  THE  CARDS  MUST  STILL  BE 

1ST  CARD  I5,3E10. 5,21 5, 2E10.5,A6 

.  2NC  C ARC  E1C.5, 15 

FOR 

THIS  PROBLEM  THE  SUPPLEMENTARY  CATA  ENTRIES  ARE, 

NIL 

DATA  ENC 


AT  A Q  IJfHM 


(tIA5  I  A  -i  HH)  AT  A3  3*  A  3  J3flT/H.) 


0  0  0  0  S-  0  i 


3  $  0  3  i  2  ?  I 


3  JT  1  1 

t a jjj/ii T/ioD)  quj  naz-uD  3ht  h  >  i33° 

.  I  .  TM310I  ■  >  '  1  '  ' '  J  83T33y-* 


(.0 


. 


j.e 

d .  i 


3.6 
SV  .  ) 


0.3 

d.t 


0.  I 

$T  .3 


3.1 

3.1 


1.3- 
3 .  i.1 


(  U  2  3  I  y  1  /1 3  VOT33V-'! 


S  I  3  1  S 

•  '  •  ?  :  ; :« .1  •  1  ri 

■ 

4«,e.ui3S..e is.e  .out  ,ei  a»*3 

I  «  ?  .  3  I  i  3  W  A  3  3/^ 

A3A3  YflATM3*3Jq<wV',9-'T  M3JH3flS  2HI  "03 

J1  A 


3/13  ttTA3 
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F  1 

NC. 

CF 

A-VECTCR  TERMS 

N  = 

15 

F  2 

NO. 

CF 

M- VECTOR  TERMS 

N  = 

5 

F  3 

NO. 

CF 

RUNS  TC  BE  MADE. 

N  = 

1 

F  <i 

ONE 

CF 

F C UR  OPTIONS  SPECIFIED 

-IF  Z-TRANSFORM  TC  BE  COMPUTED  N=0 
-TWO  SAMPLER  SYSTEM.  N=  N 

-ALL  TRANSFER  FUNCTIONS  IN  ONE  FORWARD  N  -  0 

LCOP  WITH  UNITY  FEEDBACK.  N-0 
-OTHERWISE  N  =  N  _  . . 


F  5  CNE  CF  THREE  OPTIONS  SPECIFIED. 

-FOR  ROOT  LOCUS  OF  CONTINUOUS  SYSTEM  OR 
FOR  A  SYSTEM  IN  Z-FORM.  N=0 
-FOR  ROOT  LOCI  POINTS  OF  A  ONE-SAMPLER  N  *  0 

SYSTEM  BUT  NO  Z-FORM  N=1 
-FOR  Z-TRANSFORM  COMPUTED  OR  ROOT  LOCI 
POINTS  FOR  TWO-SAMPLER  SYSTEM  N=-<10+...) 


F 6  HIGHEST  POWER  OF  S 


N  =  2 


F  7  HIGHEST  FGWER  OF  Z  N  *  0 

F 8  NUMBER  CF  VALUES  ASSIGNED  TO  T  N  =  0 

F9  ROOT  LOCI.  -,  +  ,OR  ROTH  FEEDBACK  OPIION  (  1,3,0)  N  =  I 

F 10  SCAN  CONTROL  (N  =  0,l,-1)  V  +  HtH,V  N  *  0 

Fll  MOCIFIEC  Z-FORM  OPTION  (N=0,1)  N  -  0 

F 12  REPORT  HEADING  OPTION  (N=  +  *-2)  N  =  -2 

UNUSUAL  Z-FORM  OPTION  (N=I»2) 


ei  •  * 


;  t  '  <  :  •  '  >1  •  . 


.aOASOBBB  VT!  U  HTI  *  Ip  1J 

3  7  1 


' 


*oi  r<*o  if  mi* 
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F 1 3  LOCI  OPTION,  USUALLY  N*0  _ _  N  =  0 

F14  TERNS  IN  SERIES  FOR  G*(S),N=0  GIVES  19  TERMS  N  =  -0 

F 1 6  IF  Z-TRANSFORM  TO  BE  COMPUTEO  _N.NE.O  _ 

N«CEGREE  OF  RESULTING  Z-FORM  DENOMINATOR  N  =  0 

♦  N  RCCT  LOCUS  POINTS 

F 1 7  B-MATRI*  YES  N  =  1  N  =  -0 

F 19  REAL  PART  N=- 10K  N  =  -0 


M-VECTOR  DATA 


FEEDBACK  LCOP  NO.  1 

SPECIFICATION^ OF  THE  FEEDBACK  LOOP  COMPONENTS 
AS  TO  CEGREE  OF  NUMERATOR  AND  DENOMINATOR 


UNITY  FEEDBACK 


F CRto ARC  LOOP  NO.  1 


NO.  CF  TERMS  IN  NUMERATOR  N  =  2 

DEGREES  OF  TFESE  TERMS  -VE  INDICATES  Z-FORM 
1  0 


NO.  CF  TERMS  IN  DENOMINATOR  N  =  1 

DEGREES  OF  THESE  TERMS  -VE  INDICATES  Z-FORM 
2 


2HA3T  PI  2:V!t>  »  n  231**32  A  2MP  ; 

. 


2TMI09  2U3QJ  U3« 

f  f.W  23Y  <I*U9-0  TH 


XOJ-*M 


' 


At* 


At  AO  »l  03  V  -« 


21  3^9*00  900  J  XOA60..  if  ■  L  W.3  i.  TAJ  I  *  ’.  039  _ 

1  .,,130  Q*A  *01  1*3*  J*  3*J«r3i  t.  *  2/‘ 


.  14  It'TAMHlW  141  2M*3T  93  -0/ 
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INPIT  DATA 


CONTROL  CARD  DATA  {FORMAT  2413) 


29  14 

1 

C-l  1  2 

5  110 

0  2 

0150 

0  2 

TITLE 

TRANSIENT 

RESPONSE. 

1ST  ORDER 

HOLD. 

SAMPLE 

TIME  =0.01 

MIN. 

A-VECTCR 

ENTRIES  (COEFFICIENT 

DATA. 

FORMA  T 

9E8.5  ) 

-C.5 

C. 

5  2.0 

0.0 

10.0 

LC.O 

0.01 

0.0 

3.6 

C. 

52  1.0 

5.0 

1.0 

C.  722  -  l  .  0 

1.0 

-1.0 

1. 

C  3.6 

0.92 

1.0 

0.  C 

1.0 

0.0 

1.0 

C. 

C  1.0 

0.01 

M-VECTCR  ENTRIES  (FORMAT  2413) 

1  2  C  4  1  C-l-L  5  2  1  1-1  0 

GTFER  FLAGS  ANC  INPUT  DATA.  THESE  ARE,  N VQ , OMEGA , DOMEG , 
CMEGF ,NBOC  »NZRT  »DT , T , FMT ( 1 ) .STEP, LOOP. 

THESE  VARIABLES  ARE  ENTERED  AS  ABOVE  ACCORDING  TO  THE 

_ FOLLOWING  FORMAT  IF  NCT  APPLICAB L  E  T C  THE  P R 0 B LEM 

NOTHING  NEEC  BE  ENTERED.  THE  CARDS  MUST  STILL  BE 
INCLLCEC. 

1ST  CARD  15,3E10.5,2I5,2E10.5,A6 

2ND  CARD  tlC. 5, 15 

FOR  THIS  PROBLEM  THE  SUPPLEMENTARY  DATA  ENTRIES  ARE, 


CT  =  C.Ci,  T  =  4.5,  FMT(l)  =  MIN.,  STEP  =  1.8 
CATA  ENC 


T  A Q  TJ‘1/1 


it  l  I  ■'))  ATA3  JHAD  JOflT  ,03 
S  0 10  S  0  0  I  J  ?  S  J  I -0  i.  A1  PS 

UT 11 

.  i  .  -  3mi  i  j  j-!  •  -  .  j  rai  .3 

l  c  .  3.  1  3  .  A I  /  i  ’  101  H  .  3)  23IHT/TJ  »  iTOjV-A 


.  J 

i  .0 

.  1 

3.01 

. . :) 

o.s 

?.3 

2.3- 

0.  1 

O.i- 

SSV  .3 

0.1 

.( 

0.  I 

s?.o 

d.E 

.  3 

0. 1 

J  .  J 

0.  I 

S  (  .  3 

d.E 

3.1 

O.i- 

I  0 .  J 

0.  I 

3.3 

/ 

0.  I 

(  £  1  AS 

i  A  4  "!  1  ) 

2  !  I  5-7  /ID 

V-N 

0  1- 

1  1 

S  2 

1-  1-  3 

1  A  0 

S  i 

I  ‘«/i  ] ,•! A  0DAJ3  ' 

f 2  ,  ( I) TN3  , T  9 TO «  1  ,  •  *3 

JHT  jl  Ok  I  33A  3 V : J  i A  2:.  3*37/3  3*A  23J8A1PAV  32331 

'I  31  -■  rDl  II  T  A  )3  BMItfCLIXH 

JJ1T  I2J^  2  G  -  A  0  JO  .a3^3TH3'^H  a  J  J  i/1  D/MTJ/t 

.  333  JJ3/1I 

dA, 2.01  as, 2  is, 2  .oj 3e , 2J  opa3  t 2 1 

21  ,2.313  3MA3  3 /IS 

,3flA  <3JflT2!3  A  T  A3  Y« A  1  A  3  1 3 J 4MU2  3  I T  M3J83ftS  2 1 HT  HD3 
b . I  =  4  J  f  2  ,  .  H  I  M  ^  (1)13  ,2. A  =  T  ,13.3  =  TO 


3/13  A  T  A 3 
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F  1 

NO. 

CF 

A-VECTOR  TERMS 

N  * 

29 

F2 

NO. 

CF 

M-VECTOR  TERMS 

N  * 

14 

F  3 

NO. 

GF 

RLNS  TO  BE  MADE. 

N  s 

1 

F4  ONE  CF  FOUR  OPTIONS  SPECIFIED 

-IF  Z-TRANSFORM  TO  BE  COMPUTED  N-0 
-TWO  SAMPLER  SYSTEM.  N=N 

-ALL  TRANSFER  FUNCTIONS  IN  ONE  FORWARD  N  =  0 

LOOP  WITH  UNITY  FEEDBACK.  N=0 
_ -OTHERWISE .  N  =  N 


F 5  UNE  CF  THREE  OPTIONS  SPECIFIED. 

-FOR  RCOT  LOCUS  OF  CONTINUOUS  SYSTEM  OR 
FOR  A  SYSTEM  IN  Z-FORM.  N-0 

_ -FOR  RCOT  LOCI  POINTS  OF  A  ONE-SAMPLER  N  =-ll 

SYS  TEN  BUT  NO  Z-FORM  N=l 
-FOR  Z-TRANSFORM  COMPUTED  OR  RUOT  LOCI 
POINTS  FOR  TWO-SAMPLER  SYSTEM  N=~(10+...) 


F 6  HIGHEST  POWER  OF_  S 


N  =  2 


F 7  HIGHEST  POWER  CF  Z 


N  =  5 


F  8  NUMBER  OF  VALUES  ASSIGNED  TO  T  _  _  N  -  1 

H9  ROOT  LOCI.  -,+,OR  BOTH  FEEDBACK  OPTION  (1,3,0)  N  =  1 

F 10  SCAN  CONTROL  (N  =  0,i,-i)  V  +  H,H,V  _  N  =  0 

F 1 1  MODIFIED  Z-FORM  OPTION  (N  =  0,l)  N  =  0 

F  1 2  REPORT  HEADING  OPTION  (N=+,-2)  N  =  2 

UNUSUAL  Z-FORM  OPTION  (N=l,2) 


>*OTDJV*A  ^0  .0/ 


.3CHH  38  01  enj«  10  .0* 


' 
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H  13  LOCI  OPTION, USUALLY  N=0 _ _ _ _  _ N  0 _ 


r  jl  j 

FL'» 

1  ^  r  i  i 

TERMS  IN 

Ull  *  L  1  n  — v 

SERIES  FOR  G* ( S ) ,N=0 

GIVES  19  TERMS 

N 

s  150 

E  16 

IF  Z-TRANSFORM  TO  BE  COMPUTED 

N.NE.O 

N  =  D  EGR  E  E 

OF  RESULTING  Z-FORM  DENOMINATOR 
♦N  RCOT  LOCUS  POINTS 

N 

*  2 

E  1 7 

B-MATRIX 

YES  N= l 

N 

=  -1 

E  19 

REAL  PART 

N  =- 1 OK 

N 

=  -0 

M-VECTOR  OATA 


FEEDBACK  LOOP  NO.  I 

SPEC  1  F  I  CAT  I  ON  OF  THE  FEEDBACK  LOOP  COMPONEN  T  S 
AS  TO  CEGREE  OF  NUMERATOR  AND  DENOMINATOR 


UNITY  FEEDBACK 


FORWARD  LCCP  NO.  i 


NO.  OF  TERMS  IN  NUMERATOR  N  =  A 

DECREES  OF  THESE  TERMS  -VE  INDICATES  Z-FORM 
1  0  -I  -l 


NO.  CF  TERMS  IN  DENOMINATOR  N  -  5 

DEGREES  OF  THESE  TERMS  -VE  INDICATES  Z-FORM 
2  l  1-1  0 
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APPENDIX  A. 9. 5 


TRANSIENT  RESPONSE  CALCULATION  FOR  A 
CONTINUOUS  SYSTEM  (3  order,  zero  order  hold) 


Purpose : 

To  further  illustrate  the  procedure  to  be  followed  in 
calculating  the  transient  response  of  a  system  to  a  step  input 
using  the  Z-transform. 

Problem  Source: 

This  is  essentially  the  same  problem  as  Problem  A. 9. 4.1, 
however,  a  first  order  capacity  has  been  given  to  the  control 
element,  where  previously  the  control  element  was  represented 
as  a  pure  gain. 

Block  Diagram  Representing  the  Problem: 

This  is  the  same  as  in  A. 9. 4,  the  change  is  that 


G  (s) 
v 


i  s+1 


o 


Parameter  Values 


T 


5.0 


P 


K 


1 


1.2 


K 


2 


1.5 


t .  =  0.2 

l 


T 


P 


1.0 


erid  ***a"~V*  “^i0  ^ 
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Transfer  Function: 

Open  Loop 


(s) 


3.6(0.2s+l) 

0 . 2s  (s+1)  (5.0s+l) 


Outline : 

A  root  locus  analysis  of  the  system  containing  the 
integral  time  =  0.2  showed  that  the  system  would  be  un¬ 
stable  even  for  very  low  loop  gains.  This  root  locus  is  not 
reproduced  herein,  but  by  looking  at  the  transfer  function, 
it  can  be  seen  that  poles  will  occur  at  0.0,  -1,  and  0.2  and 
a  zero  will  occur  at  -5.0.  One  root  started  at  the  -1.0  pole 
and  went  to  the  zero  at  -5.0  and  the  other  two  roots  met  and 
branched  between  -0.2  and  0.0.  Their  paths  then  crossed  the 
imaginary  axis. 

By  moving  the  controller  zero  from  -5.0  to  -0.1  the 
system  can  be  stabilized,  the  roots  starting  at  -0.2  and  -1.0 
now  meet  and  split  into  the  imaginary  plane  and  continue  in 
a  vertical  line  to  plus  and  minus  infinity.  The  system  is 
now  indicated  stable  by  the  root  locus,  (not  reproduced), 
for  all  loop  gains.  Changing  the  controller  zero  from  -5.0 
to  -0.1  is  equivalent  to  changing  the  integral  time  from  0.2 
to  10  minutes  so  that  now  the  system  transfer  function  is 


3.6  (lOs+1) 


10s (s+1) ( 5 . 0s+l) 


uiot^onu'i  !9:ien&:rT 
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The  user  must  now  calculate  the  closed  loop  transfer  function. 
Thus,  for  unity  feedback 


G 


cL 


(s) 


3.6  (lOs+1) 

50s3  +  60s2  +  46.0s  +  3.6 


The  Control  Systems  Analysis  program  cannot  transform  s- 
polynomials  of  order  greater  than  two  to  their  Z-transform 
counterparts,  therefore,  it  is  necessary  to  factor  the  denom¬ 
inator.  This  was  done  using  the  Share  library  program(23) 

2 

and  the  factors  obtained  were  50s  +  4.3752  and  s  +  1.1125s 
+  0.8226.  Thus,  the  closed  loop  transfer  function  becomes 


3.6  (lOs+1) 

gcL  (s)  =  - 5 - 

( 50s+4 .376  7)  (s  +1.1125s+0.8226) 

Following  the  procedure  outlined  in  Appendix  (A. 9. 4)  the  root 
locus  of  the  above  closed  loop  transfer  function  is  calculated 
and  plotted  to  use  as  a  future  comparison,  (Graph  9) . 

Since  the  order  of  the  denominator  of  the  s-transfer 
function  is  now  two  less  than  that  of  the  numerator  it  is 
now  possible  to  use  the  zero  order  hold  as  the  sample.  This 
sampler  is  of  the  form 


Ho(s) 


1  - 


e 


-sT 


s 


and  is  simpler  than  the  first  order  hold.  With  the  zero  order 

usually 

hold  a  built  in  lag  of  half  a  time  period  is/^present.  In 
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Graph  9 

Root  Loci  of  Continuous  and  Sampled  Systems 
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order  to  account  for  this  it  is  necessary  to  add  a  lead  of  half 
a  time  period.  Thus,  the  overall  transfer  function  in  s  is 


GcL<s)  =  esT/2Ho(s) 


3.6  (lOs+1) 


(50S+4.3762) (s  +1.1125s+0.8226) 


Now  the  denominator  polynomials  must  be  in  either  the 
2 

form  s+b  or  s  +as+b  m  order  to  be  correctly  converted  to  the 
Z-plane.  Converting  the  closed  loop  transfer  function  we  get 

sT/2, 


GCL(S)  ” 


'Hq(s)  0.072  ( 10s+l) 


(s+0.0875) (s  +1.1125S+0.8226) 


It  is  now  necessary  to  find  a  sampling  period, T  to  give 
a  root  locus  for  the  sampled  system  in  this  plane  which  will 
correspond  to  the  continuous  case.  In  this  way,  the  effect 
of  the  sampler  is  done  away  with  and  the  sampled  system  will 
correctly  simulate  the  continuous  case.  This  corresponds  to 
step  4  in  the  problem  of  Appendix  (A. 9. 4) . 

From  Graph  9,  which  shows  the  root  loci  of  both  the 
sampled  and  continuous  cases,  it  can  be  seen  that  for  a 
sampling  period  T  =  0.1,  the  root  loci  trajectories  and  poles 
match  almost  exactly.  The  gain  curves  were  also  close  but 
not  exact;  these  are  not  shown.  This  match,  however,  was 
deemed  sufficient  (Section  8. 5. a),  to  give  a  transient  res¬ 
ponse  solution  closely  approximating  the  true  solution.  It 
was  expected,  though,  that  the  imperfect  match  in  trajectory 
and  gain  would  contribute  to  both  a  phase  shift  and  a  dif¬ 
ference  in  amplitude.  Because  the  poles  and  zeros  are  matched. 
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the  form  of  the  response  or  the  oscillatory  character  of  the 
digital  solution,  should  be  the  same  as  that  for  the  con¬ 
tinuous  solution. 

Steps  6  and  7,  (Appendix  A. 9. 4) ,  were  then  carried 
out  so  that  the  transient  curve  no.  (1)  graph  (10)  was  ob¬ 
tained.  The  analogue  solution  is  shown  as  curve  no.  (2) 
graph  (10) .  The  expected  result  was  obtained.  There  is  pre¬ 
sent  a  small  shift  both  in  phase  and  amplitude,  however,  the 
oscillatory  response  is  quite  similar.  For  process  control 
work,  where  the  oscillatory  response  is  the  main  concern  to 
the  designer,  an  error  of  this  magnitude  can  be  tolerated. 
This  error  may  be  eliminated  by  smaller  T,  however,  round-off 
error  and  computation  time  may  become  excessive.  It  should 
be  noted  that  when  a  zero  order  hold  is  used,  a  lead  of  less 
than  a  full  sample  period  is  employed.  This  makes  it  neces¬ 
sary  to  put  a  1  in  position  JO  (11)  of  the  control  card. 

The  input  data  used  for  this  problem  is  listed  on 
pages  137  to  143.  The  disturbance  was  a  step  =  2.0. 
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o  3  4J  ad  u  s  ^  • 


Time  (min 
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INPUT  DATA 


CONTROL  CARD  CATA  (FORMAT  2-413) 

17  6  1C02CG1CG  -2  0  0  0  C 

TITLE _ 

ROOT  LOCUS  OF  TFE  CLOSED  LOOP  (CONTINUOUS) 
A-VECTCR  ENTRIES  (COEFFICIENT  DATA.  FORMAT  9E8.5) 


-0^02 _ CU_2 _ C^Ol _ 0,0  2.0 _ 4_,C _ 

L.O  1C.0  3.6  0.82263  1.1125  l.C  4.3762  50. 0 

M-VECTCR  ENTRIES  (FORMAT  2413) 

2  1  C  2  2  1 

OTHER  FLAKS  AND  IN  PUT  C  ATA.  THE$_E  ARE, N YQ , OMEGA , DCMEGt 

CMEGF, NBObfNZRTtDTt T,FMT( 1 )  , STEP, LOOP. 

THESE  VARIABLES  ARE  ENTERED  AS  ABOVE  ACCORDING  TO  THE 
FOLLOWING  FORMAT  IF  NOT  APPLICABLE  TO  THE  PROBLEM 
NOTHING  NEED  BE  ENTERED.  THE  CARDS  MUST  STILL  BE 

incllcec. 

1ST  CARD _ I5.3E10.5,2I5,2E10.5,A6 _ 

2ND  CARD  EIC.5,15 

FOR  THIS  PROBLEM  TFE  SUPPLEMENTARY  DATA  ENTRIES  ARE, 


NIL 


CATA  ENC 


_ 


— 


A  T  AO  TJ9/I 


WMflOl)  AT  AO  ]  H  A  D  JlflT/133 

jOS-00  13  )  S  0  0  !•  d  V I 

3JTIT 

l  2  JOIM  I  T  AID  )  M)0  J  C1323J0  3  IT  2U00J  TJOH 

(  .  1  •  A 1  \  ; 

j  .  ^  0*S  0.0  10.3  •  3  SO.D- 

j  ,  jc  SdY  t  J  .  1  <?  S  1  1  .  i  £  d  S  S  8  •  Q  ' .  t  .  .  >J  i  J  •  I 

(  !  T  'll  'i  U  *3  1  T  3  :J  *  T  0  i  V  -  1 

i  s  s  o  i  s 

« i)  j  ^  i'J  t  A\.  •■/) »  OY  3  t3MA  3  2  3  H  T  -MAO  T  UM/I I  OlrfA  20AJ3  V.OHTO 

.  H  i  )  J  »  4  J I  :  ,  ( 1  )  TM3  ,  T , T  I «  T '  \ /  » 303H* 3D3^ 1 

jHl  JT  0,  I  v  J  3  A  iv  iA  2  A  C1383T/13  J*A  23J‘3A1«AV  32  3  31 
,,-jj  j«M  3MT  )T  3J3A31JHMA  TO/  U  T  A  4P  33  J/IW0JJJ3 
^3  JJIT2  I2J3  23flAD  3HT  .03fl3T/l3  33  3333  OlrfMTD/1 

.  333  J  lOlAl 

dA,c,oi3Steis.e.ou£tei  3*ao  T2i 

-  .15*: 

,  .  |  .  •  A  I  ■  1  •  !  •  : 

31/ 


333  AT  A3 
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FI  NO.  CF  A-VECTCR  TERMS  N  *  17 


F 2  NO.  CF  M- VECTOR  TERMS  N  *  6 

F 3  NO.  GF  RLNS  TO  BE  MADE.  N  =  1 


FA  CNE  CF  FCUR  OPTIONS  SPECIFIED 

-IF  Z— T RANSFORM  TO  BE  COMPUTED  N*0 
-ThC  SAMPLER  SYSTEM.  N=N 

-ALL  TRANSFER  FUNCTIONS  IN  ONE  FORWARD  N  *  0 

LCOP  MTH  UNITY  FEEDBACK.  N^O 
-OTHERWISE _ N-N _ _ _ _ 


F 5  ONE  GF  THREE  OPTIONS  SPECIFIED. 

-FCR  RCOT  LOCUS  GF  CONTINUOUS  SYSTEM  OR 
FCR  A  SYSTEM  IN  Z-FORM.  N=0 
-FCR  RCCT  LCC I  PC  I  NTS  OF  A  ONEj-SAMPLER  N  *  0 

SYSTEM  BUT  NO  Z-FORM  N*t 
-FCR  Z-TRANSFORM  COMPUTED  OR  ROOT  LOCI 
PCINTS  FOR  TWO-SAMPLER  SYSTEM  N=- < LO+ . . . ) 

Ffc  HIGHEST  POWER  CF  S _  N  *  2 

F 7  HIGHEST  POWER  GF  Z  N  =  0 

H B  NUMBER  CF  VALUES  ASS  IGNEO _ TO  T _ N  =  0 


F 9  RCOT  LOCI.  -,4t0R  BOTH  FEEDBACK  OPTION  (l,3tO)  N  =  l 
F  1 0  SC  A  N  CONTROL  _ ( N  =  0  ,  l ,  -  l  )  V+H,Hj>V _ _ _  N  0 


F 1 1  MODIFIED  Z-FORM  OPTION  (N=0#l)  N  =  0 

F 1 2  REPORT  HEADING  OPTION  (N  =  +  f-2)  _  N  =  -2 

UNUSUAL  Z-FORM  OPTION  (N=L*2) 


■ 
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F  13 

LOCI  OPT  ION, USUALLY  N  =  0 

N  * 

0 

F  14 

TERMS  IN  SERIES  FOR  G*(S),N=0 

GIVES  19 

TERMS  N  - 

0 

F  16 

IF  Z-TR  0NSFORM  TO  BE  COMPUTED 

N.NE.O 

N=CEGREE  OF  RESULTING  Z-FORM  DENOMINATOR  N  = 0 


♦  N  RCGT  LOCUS  POINTS 


F  17 

B-MATRIX 

YES 

N  =  1 

F  19 

REAL  PART 

N=- 

10K 

M-VECTOR  DATA 


FEEDBACK  LCCP  NO.  I 

SPECIFICATION  OF  THE  FEEOBACK  LOOP  COMPONENTS 
AS  TO  CEGREE  OF  NUMERATOR  AND  DENOMINATOR 


UNITY  FEEDBACK 


FORWARD  LOOP  NO.  I 


NO.  CF  TERMS  IN  NUMERATOR  N  -  2 

DEGREES  OF  THESE  TERMS  -V E  INDICATES  Z-FORM 
1  0 


NO.  CF  TERMS  IN  DENOMINATOR  N  *  2 

DEGREES  CF  THESE  TERMS  -VE  INDICATES  Z-FORM 
2  1 


■ 


, 
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INPLT  DATA 


CONTROL  CARD  DATA  (FCRMAT  2413) 


3C  14  1 

TITLE 

C-1C  3  5 

2  1 

0  1-2  C  50  0  3 

TRANSIENT 

RESPONSE. 

ZERCE 

ORDER  HOLD. 

A-VECTCR  ENTRIES  (COEFFICIENT  DATA.  FORMAT  9E8.5) 


-0.05 

C  .  2 

1.0 

C.O 

2.0 

4.  C 

0.1 

0.05 

0.82263 

1.1125 

1.0 

0.0875 

1.0 

l.C 

10.0 

.072 

-1.0 

l.C 

C.O 

1.0 

0.82263 

1.  1125 

1.0 

C.0875 

1  .C 

C.O 

1.0 

C.O 

l.C 

M-VECTCR 

ENTR  IES 

(  FORMAT 

2413) 

2  2  1  0  4  1  C  -1-21  4  2  11-1 


OTHER  FLAGS  AND  INPUT  DATA.  THESE  ARE,  N YQ , CM EG A , DON  EG , 
CMEGF *NBOC ,NZRT » DT, T»  FMT( 1 ) , STEP, LOOP. 

THESE  VARIABLES  ARE  ENTERED  AS  ABOVE  ACCORDING  TO  THE 
F  C  L  L  C  k  I  N  G  FCRMAT  IF  NCT  APPLICABLE  TL  THE  PROBLEM 
NOTHING  NEED  BE  ENTERED.  THE  CARDS  MUST  STILL  BE 
INCLUDED. 

1ST  CARD  I5,3E10.5,2I 5,2E10.5,A6 

2NC  CARD  E1C.5,I5 

FOR  THIS  PROBLEM  THE  SUPPLEMENTARY  DATA  ENTRIES  ARE, 


D  T  =  0  •  1 ,  T  =  1 C  •  C  ♦  FM  T  (  1  )  =  M  I N.  »  STEP=2.0 
DATA  END 


. 

\  . 


JH1 


I  .0 

j .  I 


ATAO  TJ9I/I! 

^  <  I )  ATAO  i^AO  J  3  H  T  /10  0 

o  oe  o  s-  i  o  i  $  £  31-*3  1  M  s£ 

3  JT  I  T 

•  OiOH  *  092  III  TH3I2WAHT 


s . o  ao.o- 

eSii.l  S3SSB.0 

3.1  o.i- 

j.I  2V80.0 


1  v  ,,H  .  AT  A  Cl 

T  /  J  I  Jl  »  TifO  ) 

J  •  £ 

o.s 

0.  i 

•  A 

0.  I 

rJ  Y  rt  0  .  J 

0.1 

I  i  .  1 

tdSS  ..0 

. 

0.3 

J  .  I 

J .  J 

0.  I 

0  .  j 

(  £  1 

t/ A^n  \ ) 

l-  I  i 


1S-J-  0  i  ^  0  1  *  s 


\z  JHT  .ATAQ  TU9HI  OiAA  23  m  mjhTO 

’  .  .  ...  •  ■ 
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a&,«!.oi35«ais.e.oi3t  ,ei  T^i 

eitr.013  3«A0  3IAS 

,  dl 1*11/13  ATAO  Y«ATH3M3J^^0a  3“»T  M3J83HM  2IHT  «□=» 
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FI  NO.  OF  A-VECTOR  TERNS 


N  =  30 


F2  NO.  CF  M- VECTOR  TERNS 


N 


=  14 


F  3  NO.  CF  RLNS  TO  BE  NAOE. 


N  =  1 


F4  CNF.  CF  FCUR  OPTIONS  SPECIFIED 

-IF  Z- TRANSFORM  TO  BE  COMPUTED  N=0 
-TWO  SAMPLER  SYSTEM.  N=N 

-ALL  TRANSFER  FUNCTIONS  IN  ONE  FORWARD  N  =  0 

LCOP  WITH  UNITY  FEEDBACK.  N*0 
-OTHERWISE  N=N 


F  5  ONE  CF  THREE  CPTIONS  SPECIFIED. 

-FCR  RCOT  LOCUS  CF  CONTINUOUS  SYSTEM  OR 
FCR  A  SYSTEM  IN  Z-FORM.  N=0 
-FCR  RCOI  LOCI  POINTS  OF  A  ONE-SAMPLER  N  =-10 

SYSTEM  BUT  NO  Z-FORM  N=1 
-FCR  Z-TRANSFORM  COMPUTED  OR  ROOT  LOCI 
POINTS  FOR  TWO-SANPLFR  SYSTEM  N=-(10+...) 


F6  HIGHEST  POWER  CF  S 


N  =  3 


F 7  HIGHEST  POWER  OF  Z 


N  =  5 


F 8  NUMBER  OF  VALUES  ASSIGNED  TO  T 


F9  RCOT  LOCI.  -,+,0R  BOTH  FEEDBACK  OPTION  (1,3*0)  N  =  1 

F  10  SCAN  CONTROL  (N  =  0, 1,-1 )  V  +  H,H,  V _  N  * .  0 

H 1 1  MOO  I F I  EC  Z-FORN  OPTION  (N  =  0,l)  N  =  1 

FI  2  REPORT  HEADING  OPTION  (N=+,-2) 

UNUSUAL  Z-FORM  OPT ICN  (N*l ,2) 


N 


2 


.  • 


3 
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M3  LOCI  CP  I  ICN,  USUALLY  N  =  0  N  =  0 

FI4  TERNS  IN  SERIFS  NCR  G*(S),N=0  GIVES  19  TERMS  N  =  50 

\  16  If  /-TRANSFORM  Tfl  BE  COMPUTED  N.NE.O 

N=CEGREE  OF  RESULTING  Z-FORM  DENOMINATOR  N  =  3 

4N  ROOT  LOCUS  POINTS 


F  1  7 

R-MATR  I  X 

YES  N= 1 

N 

*  -0 

F  1  9 

REAL  PART 

N  =  -  l  OK 

N 

*  -0 

M-VECTOR  DATA 


NO FEEDBACK 


FORWARD  LOOP  NO.  1 


NO.  CF  TERMS  IN  NUMERATOR  N  =  A 

DEGREES  CF  THESE  TERMS  -VE  INDICATES  Z-FORM 
1  0  -1-21 


NO.  CF  TERMS  IN  DENOMINATOR  N  *  A 

DEGREES  OF  THESE  TERMS  -VE  INDICATES  Z-FORM 
2  11-1 


n  j 
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APPENDIX  A. 9. 6 

TRANSIENT  RESPONSE  OF  A  SAMPLED-DATA  SYSTEM 
DIGITAL  CONTROLLER  DESIGN 


Purpose : 

To  illustrate  the  procedure  followed  in  obtaining  the 
response  of  a  digitally  controlled  system  to  a  step  change  in 
set  point,  and,  in  the  design  of  a  digital  controller  using 
the  Z-plane. 

Transfer  Function: 

1  -  e~sT  0.5 

H  (s)G(s)  =  -  -  (1) 

°  p  s  (s+1) (5s+l) 

Outline:  First  Example: 

1.  Take  the  Z-transform  and  find  the  Z-plane  root  locus  for 
transfer  function  (1) .  To  do  this  a  sampling  period  must 
be  decided  on.  For  this  problem  a  sampling  period  of  one 
minute  was  chosen.  The  computed  Z-transform  is 

0.03428Z  +  0.023014 

(H  G  )  (Z)  =  - 

p  Z  -  1.1866Z  +  0.301194 

The  root  locus  was  also  computed  for  H  G  (Z)  in  the 

op 

Z-plane.  This  is  shown  in  Graph  (11) . 

2.  Using  this  root  locus  a  controller  can  be  specified  in 
the  Z-plane  which  will  control  the  system.  In  order  to 
obtain  zero  offset  at  steady-state,  a  pole  must  be 
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specified  at  Z  =  1.  Zeros  must  not  exceed  the  poles  in 
the  controller,  therefore,  one  zero  can  be  specified. 
Examining  the  root  locus  of  HQGp(Z),  Graph  (11),  it  can 
be  seen  that  system  poles  occur  at  Z  =  +0.819  and  Z  =  0.368 
and  a  system  zero  at  Z  =  -0.669.  If  a  controller  pole  is 
added  to  the  system  at  Z  =  1 ,  the  instability  of  the  sys¬ 
tem  is  greatly  increased.  To  offset  the  destabilizing  in¬ 
fluence  of  the  integration  the  controller  zero  must  be 
placed  between  Z  =  1  and  Z  =  0.819.  The  controller  root 
starts  at  Z  =  1  then  goes  to  the  controller  zero.  The 
two  other  roots  of  the  system  go  towards  each  other,  meet, 
split,  eventually  cross  the  unit  circle  boundary,  then 
circle  back  to  the  real  axis  where  one  goes  to  the  zero  at 
Z  =  -0.669  and  one  goes  to  minus  infinity.  Graph  (12). 

3.  Let  the  controller  zero  be  placed  at  Z  =  0.834.  The 
placing  of  this  zero  decides  the  integral  time  of  the  con¬ 
troller.  As  the  zero  is  placed  closer  and  closer  to  the 
pole,  the  less  integral  action  there  is  in  the  system. 

This  placing  of  the  zero  is  done  mainly  by  trial  and  error, 
no  criteria  has  been  followed. 

4 .  Now 

0  (Z)  KD ( Z ) G  ( Z ) 

I  (Z)  1  +  KD  (Z)  G  (Z) 

D ( Z )  has  been  designated  as 

Z  -  0.834 


D(Z) 


Z 


1.0 
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K  must  be  decided  upon. 

It  was  indicated  by  the  Root  Locus  for  the  controlled 
system,  (Graph  12) ,  that  for  K  =  1.0  the  system  would 
be  slightly  underdamped.  The  gain  was  printed  out  on 
paper  but  is  not  included  in  Graph  (12) .  To  test  the 
correspondence  between  the  root  locus  prediction  and  the 
transient  response  calculation,  the  value  K  =  1.0  was 

chosen  for  the  transient  calculation.  This  is  the  first 
trial,  the  input  data  used  for  this  example  is  listed 
on  pages  153  to  158. 

5.  Calculate  the  closed  loop  transfer  function  in  Z. 

This  operation  is  done  by  the  machine  if  there  is  unity 
feedback.  To  enter  this  calculation  branch,  the  flag 
LOOP  must  be  set  equal  to  one,  (Appendix  A. 3.1,  A. 5.1). 
This  transfer  function  is  not  printed  out. 

6.  The  transient  response  curve  is  shown  as  Graph  (13) . 

For  a  loop  gain  equal  to  one  the  roots  of  the  character¬ 
istic  equation  are  all  real,  and  there  is  no  overshoot 

in  the  response,  Graph  (13). 

\ 

The  disturbance  was  applied  to  the  set-point  variable  and 
was  a  step  of  1. 


Summary  of  Procedure: 

1.  Obtain  the  Z-transform  of  the  process  and  hold,  that  is 
obtain  HG^(Z)  using  the  C.S.A.  program. 
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.  S  ni  noldoru/i  i.->ic-nBid  qool  beeol o  jl Id  ©dsluo.  6 


.duo  bsdnxiq  don  it  nc  onul  ssianedd  airiT 


.  d 
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2.  Obtain  the  Z-transform  of  the  selected  digital  controller. 
This  is  formed  by  the  user  using  the  Z-plane  root  locus 

or  any  other  criteria  he  wishes. 

3.  Calculate  the  closed  loop  transfer  function  in  Z  and  then 
invert  the  transfer  function  to  get  the  transient  response. 
If  there  is  unity  feedback  and  the  system  is  error- 
sampled,  the  LOOP  option  may  be  used. 

Second  Example: 

For  this  example,  the  zero  of  the  controller  is  changed 
from  0.834  to  0.95.  The  root  locus  for  this  open  loop  trans¬ 
fer  function,  controller  included,  is  calculated.  Graph  (14) . 
From  this  root  locus  a  loop  gain  of  12.0  is  chosen.  For  this 
gain  the  roots  are  complex  but  inside  the  unit  circle,  so 
that  an  oscillatory  stable  response  should  result.  The  tran¬ 
sient  response  was  calculated,  as  previously,  and  the  result¬ 
ing  transient  is  shown  in  Graph  (15) . 

The  step  input  to  the  set-point  had  a  value  of  one  and  for 
this,  the  overshoot  is  approximately  30%.  The  integral  action 
of  the  controller  is  indicated  by  the  decreasing  offset.  In 
this  example  the  integral  action  is  very  slow  so  that  the 
controller  zero  should  be  moved  further  to  the  left  on  the 
real  axis  to  provide  more  integral  action.  When  this  is  done, 
the  corresponding  root  locus  should  be  calculated  in  case  the 
gain  of  12.0,  which  was  used  here,  be  too  high  for  the  modi¬ 
fied  system.  The  input  data  for  this  problem  is  listed  on 
pages  159  to  161 . 
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Step  Input  =  1.0 


V 


32!  0 


0. 


8 .0 


16.0  Z+.O 

Time  (min.) 


4-0.0 
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INPUT  DATA 


CONTROL  CARD  CATA  (FORMAT  2  A  1 3 ) 

24  12  1  C-10  141100-2  0100  0  2 

TITLE _ 

ROOT  LOCUS  OF  SYSTEM  WITH  SAMPLE-HOLD  INSERTED. 
A-VECTCR  ENTRIES  (COEFFICIENT  DATA.  FORMAT  9E8.5) 


— C  •  l 

0.2 

1.0 

o 

• 

o 

2.0 

4  .  C 

1.0 

0.2 

l.C 

1.0 

1.0 

0.1 

-1.0 

1.0 

0.2 

1.0 

1  .C 

1.0 

0.0 

l.C 

C.C 

1.0 

M- VECTOR  ENTRIES  (FORMAT  2413) 

2  1  1  0 _ 2  0  -l  4  111-1 

OTHER  FLAGS  AND  INPUT  DATA.  THESE  ARE,  NYU , OMEGA , DOMEG, 
CMEGF,NBCC,NZRT,DT,T,FMT(  1  )  , STEP, LOCH. 

THESE  VARIABLES  ARE  ENTERED  AS  ABOVE  ACCORDING  TO  THE 
FOLLOWING  FORMAT  IF  NOT  APPLICABLE  TO  THE  PROBLEM! 
NOTHING  NEED  BE  ENTERED.  THE  OARDS  MUST  STILL  BE 

_ INCLUDED. _ _ 

1ST  CARC  I  5, 3E10. 5,21 5,2E10.5,A6 

2NC  CARC  E10.5, 15 

FOR  THIS  PROBLEM  THE  SUPPLEMENTARY  CATA  ENTRIES  ARE, 

N  ZR  T  =  1 


DATA  ENC 


— 


AT  AO  TJM^I 


(  I  AS  r  A « 3  I  )  i  T  A  :i  .1  fl  A  0  J  ) « T  / 1 ') 

S  0  JjIj  s -  0  0  1  i  A  I  )!-)  I-  ^ 

3J1  IT 


,  |  i  -jjmma,'  mi^  H3T2Y2  y]  JU30J  T03H 

•  AT AO  TH3I3I  23I#Tirt3  - 


. 


0.  i 
L  .  I 
0.  i 


J  .  »* 
0.  I- 

J  .  J 


j.s; 
i  .0 
o.  i 


0.0 

0.  I 

0.0 


o.i 
o.  i 
o.i 


s  .0 

0.  1 
J.  I 


1.3- 
S.  0 

O.i 


(£Ii  ]  l  r-  I  *3  A0TD3V-H 


I-  I  i  i  A  I-  0  <>  0  11S 

» 0 3 M 00  »  AO 3 MO  »  0  Y /  «  J«A  383HT  .ATAO  TUM/I  3/A  23AJ3  H3HT0 

.ijjj,  j  T  <  ’  t (  i  )  T  M  3  »  T  » I Q  » TflV/t TOO/*  I0DM3 
-jHT  Cl  Or  i  I  Cl  >4  OJA  JvUdA  8  A  33H3T/3  3«A  8  3  J  M  A  I  WAV  3833T 

_HT  .if  3  J  J  A  0  1  J  3  3  A  TJ/  11  TAMAJ3  0/IW3JJ03 

.  0 3  >  i  T  /  j  3  3  33  3M  O/HI 

•  333 J JO/ 1 

dA,c*.oi  is.eis,e. 013^*81  o^ao  t 2 1 

£1,8.313  OflAO  3/S 

f  .  >; ,  -j  |  >  1  ■  j  A  I  A3  Y>.  AT  'A  AM  3  JS9  J8  3  3T  MJJM'1^3  8131  «  '3 


i  -  i  a  \  a 


3/3  ATAO 
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F  I  NO  .  CF  A-VECTCR  TFRMS 


N  *  24 


H2 

NC.  CF 

M-VECTCR  TERMS 

N  =  12 

F  3 

NO.  CF 

RUNS  TC  BE  MADL. 

N  =  l 

F 4  ONE  CF  FC UR  OPTIONS  SPECIFIED 

-IF  Z- TRANSFORM  TC  BE  COMPUTED  N=0 
—  T V> G  SAMPLER  SYSTEM.  N=N 

-ALL  TRANSFER  FUNCTIONS  IN  ONE  FORWARD  N  =  0 

LOOP  WITH  UNITY  FEEDBACK.  N  =  0 
-OTHERWISE  N  =  N 


F 5  ONE  CF  THREE  OPTIONS  SPECIFIED. 

- FC R  ROOT  LOCUS  OF  CONTINUOUS  SYSTEM  OR 
FC R  A  SYSTEM  IN  Z-FORM.  N  =  0 
-FCR  ROOT  LOCI  POINTS  OF  A  ONE-SAMPLER  N  =-10 

SYSTEM  BUT  NO  Z-FORM  N=1 
-FCR  Z-TRANSFORM  COMPUTED  OR  ROOT  LOCI 
PCINTS  FOR  TWO-SAMPLER  SYSTEM  N=-(10+...) 


Ft  HIGHEST  POWER  CF  S 


N  =  I 


F 7  H  IGHEST  FEWER  CF  Z 


N  =  4 


F  8  NUMBER  OF  VALUES  ASSIGNED  TO  T  N  =  1 

F 9  ROOT  LOCI.  -,+,CR  BCTH  FEEDBACK  OPTION  (lf3t0)  N  =  l 


F 10  SCAN  CONTROL  (N  =  0,l,-l)  V  +  H,H,V  N  =  0 

Fll  MOD  I F  I  EC  Z-FORM  OPTION  (N  =  0,i)  N  =  0 

F 12  REPORT  HEADING  OPTION  (N  =  +  ,-2)  N  =  -2 

UNUSUAL  Z-FORM  OPTION  (N=l,2> 


•  I 


/,  n  '■ .  r;  <■  ?  ;  - 

0*tf  .»DA8J33-I  YT1*U  HT  I  W  <T3*'J 


' 

IOOJ  TOD«  *0  q3Ta^MOO  M#P13WA*t-i  *H3~ 
(...♦Oil-  «  M3T2YE  fl3JMWA?-hlT  * 03  211^13^ 

•  1  i  13  »3W0=»  TZJHdn  <H 


*  U 


T  OT  ?3JJAV  33  33 

]  *  M  ( 0 « € « i  )  H31T93  )OA80333  rlT3Q  *0t^-  .1D0J  TOO*  ?3 


V,H,H*V  J3STV1D3  H02  01=1 
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F  1  3  L  0  C  I  OP  T  IC N ,  USUAL L  Y  N*0_ _ _ _  _N _* 0 

F  14  TERNS  IN  SERIES  FCR  G*(S  )  fN=0„ GIVES  19  TERMS  N  =100 

F 16  IF  Z-TRANSFORM  TO  BE  COMPUTED _ N.NE.O  _  _ 

N=CEGREE  OF  RESULTING  Z-FORM  DENOMINATOR  N  =  2 

4N  RCOT  LOCUS  POINTS 

F 1 7  B-M  A  TR I >  YES  N=1  N  =  -0 

F  19  REAL  PART  N=-  10K  N  =  -0 


M-VECTOR  DATA 


NOFJEDDACK 


FORWARD  LOOP  NO •  1 


NO.  CF  TERMS  IN  NUMERATOR  N  =  2 

DEGREES  OF  THESE  TERMS  -VE  INDICATES  Z-FORM 
C  -1 

NO._CF  TERMS  IN  DENOMINATOR  _N__r _ A  _ 

DEGREES  OF  THESE  TERMS  -VE  INDICATES  Z-FORM 
1  l  1  -1 


TPA9  Jaba 

' 


ATA3  SOT33V-M 


f  .OH  <1jJJ  OAAirfflOq 


HA03-S  2aTA010MI  3V- 


-  1  56- 


INPLT  DATA 


CONTROL  CARC  CATA  ( FORMAT  2413) 

24  12  1  C-1C  141100-2  01C0  0  2 

TITLE 

TRANSIENT  RESPONSE.  DIGITAL  CONTROLLER  INCLUDED. 


A-VECTCR  ENTRIES  (COEFFICIENT  DATA.  FORMAT  9E8.5) 


-C.  1 

C  .  2 

l.C 

0.0 

2.0 

4.  C 

l  .0 

0.2 

l.C 

1.0 

1.0 

0.  1 

-1.0 

l  .0 

0.2 

l.C 

l.C 

1.0 

c.o 

1.0 

C.C 

1.0 

M- V  EC  TO  R  ENTRIES  (FORMAT  2413) 

21  1C2C  -14111-1 

OTHER  FLAGS  AND  INPUT  DATA.  THESE  ARE,  N VQ , OMEGA , DOMEG t 
CMEGF,NBOO,NZRT,DT,T,FMT( 1  )  , STEP, LOOP. 

1HESE  VARIABLES  ARE  ENTERED  AS  ABOVE  ACCORDING  TO  THE 
FOLLOWING  FORMAT  IF  NOT  APPLICABLE  TO  THE  PROBLEM 
NOTHING  NEED  BE  ENTERED.  THE  CARDS  MUST  STILL  BE 
INCLLCEC. 

IS T  C  A R C _ 15 ^ 3 E 1 0 . 5 , 2  I  5 , 2 E 1 0 . 5 , A 6 

2NC  CARC  E1C.5,  15 

FOR  THIS  PROBLEM  THE  SUPPLEMENTARY  CATA  ENTRIES  ARE, 

NZR  T= 1 ,  C  T  =  l  .  C ,  T  =  60 . 0 ,  FMT(1)=MIN.  STEP=1.C,  LCOP=l 

Ml  VECTOR  (FORMAT  1013) 


1111 

A 1  VECTOR  (FORMAT  9E8.5) 


-0.834  l.C  - 1 . 0  1.0 


FOR  INFORMATION  ON  LOOP  =  1  DATA  INPUT 
REFER  TO  SUBROUTINE  AMALG ,  APPENDIX 

A- 


CATA  ENC 


AT  AC)  7JS/1 


(£lt*S  I  A  ■  t,  |)  T  A  )  HAD  1 3  5A  7  /T 11 ; : 

S  0  0 J 10  s-  0  3  I  I  A  i  Oi-3  I  S!  AS 

3J717 

•  CJOUJD*  v.jJ  ) >  7  /)  j 0  J  \  7  I  - 1  •  .32  3023*  7  /! 3  I  ^  ' A  :  T 


l 

♦ 

b  3 1  7A/M  H 

.  A  T  A  C I 

7/31  3133  jH3  ) 

2  3  1  9  T  /I  3 

'•  JT3JV-A 

u  • 

i 

J  .  A 

o.s 

0.0  0.1 

s  • 

0  1.3- 

. 

j 

0.  i- 

1  .0 

0.1  0.1 

3. 

S.O 

o . 

i 

3  .  J 

J .  i 

0.0  o. 

3. 

1  0.1 

(  £  )  AS  7  a  /  >  n 

231*7/3 

- 

1-  1 

1  1  A  i-  0 

'  S  3 

I  I  s 

,  v)  ,a  j  ?  J  /  /  i'‘  \  7  3h7  .A!  A]  C  JM/I  3  A  2DAJ1  «JHT3 

.  U  .UtO  12  ,  ( 1  )TM3f 7, 70, 7  '♦  T.7  i/,  3  3  3/3 

jH7  07  «  1C  0  3.  JVOdA  2A  Q3M37/3  IMA  23J3AIMAV  J2HI 

M  -)  J  Ti  M  ^  -ill  JT  dJIAOIJMSA  10H  31  7AMMH  0/IWJJJ37 
i  j  1  1  ].  .■■■  ..  .  •'  1  18  3  j  i  \A  0/1  3 1  J  / 

•  3  j  3  j j on i 

aAfe.  >ids,c  !Ste.oiif  ,ei  3«a3  t 2 r 

, c*  .>1  3  1*A3  3/S 

,j  A  2  J  1  >  I  /  J  *■>  7  A  J  YMAT/3M3  J<19J2  337  M  3  J  8  3  M  9  2IHT  M33 

. .  I  •  -  (  I  )  T  M  7  »  .  «  0  •  I  *  1 3  *  1  =  T  M  A  / 

(U  01  7A/M03 )  M3703V  1/ 

I  1  1  1 

(  .339  TAMM03)  M3T33V  IA 

0.1  0.1-  3.1  A£8.0- 

[  u4iy  I  A  T  AG  1  =  900J  /O  /□I7AMMD3/1  M03 

XI  i  J  0  0  A  t  0  J  A/A  3UTUOMOJ2  07  M333H 

-A 


3/3  A  7  A 3 
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FI  NO.  CF  A-VECTCR  TERMS  N  -  24 

F  2  NC.  CF  M-VECTCR  TERMS  N  =  12 

F  3  NC.  CF  RINS  TO  BE  MADE.  N  =  l 


FA  CNE  CF  FCUR  OPTIONS  SPECIFIED 

-IF  Z-TRANSFCRM  TC  BE  COMPUTED  N=0 
—  T  V\  C  SAMPLER  SYSTEM.  N=N 

-ALL  TRANSFER  FUNCTIONS  IN  ONE  FORWARD  N  -  0 

LCCP  WITH  UNITY  FEEDBACK.  N=0 
-OTHERWISE  N  =  N 


H5  CNE  CF  THREE  OPTIONS  SPECIFIED. 

-FCR  RCCT  LOCUS  CF  CONTINUOUS  SYSTEM  OR 
FCR  A  SYSTEM  IN  Z-FORM.  N  =  0 
-FCR  RCCT  LOCI  POINTS  OF  A  ONE-SAMPLER  N  =-10 

SYSTEM  BUT  NO  Z-FORM  N=1 
-FCR  Z-TRANSFCRM  COMPUTED  OR  ROOT  LOCI 
POINTS  FOR  TWO-SAMPLER  SYSTEM  N=-(10+...J 


H 6  HIGHEST  FEWER  CF  S  N  =  1 

F 7  HIGHEST  POWER  CF  Z  N  =  4 

\Q  NUMBER  OF  VALUES  ASSIGNED  TC  T__  N  =  1 


H  S  ROOT  LCC  I .  -,+,CR  BCTH  FEEDBACK  OPTION  11*3, 0)  N  =  l 


F 1 0  SCAN  CONTROL  (N=0*l,-l)  V+H,H,V 


Eli  iVGCIFIEC  Z-FORM  OPTION  (N=0,1) 


N  =  0 


F l 2  REPORT  HEADING  OPTION  (N=+,-2) 
UNUSUAL  Z-FCRM  OPTION  (N=-l,2) 


N  =  -2 


1 


' 


A  -0  ?)"  >*>  IM.J  131*  '  ' 
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M3  LOCI  CPI  ICN  t  USUAL  L  Y  N=0  N  =  __  0 


F  14 

TERNS  IN 

SERIES  FCR  G*(S)tN=0 

GIVES 

19  TERMS 

N  =100 

F  16 

IF  Z- TR  £NSFORN  TO  BE  COMPUTED 

N.NE  • 

0 

N=CEGR  EE 

OF  RESULTING  Z-FORM  DENOMINATOR 
♦ N  ROOT  LOCUS  POINTS 

N  =  2 

F  17 

B-Mfl  IR  I  > 

YES  N  =  1 

N  =  -0 

F  19 

REAL  PART 

N  =  -  10K 

N  =  -0 

M- VECTOR  DATA 


\'P  FEEDBACK 


FORWARD  LCCP  NC.  1 


NC.  CF  TERNS  IN  NUMERATOR  N  =  2 

DEGREES  CF  THESE  TERMS  -VE  INDICATES  Z-FORM 
C  -1 


NC.  CF  TERMS  IN  DENOMINATOR  N  =  4 

DEGREES  CF  THESE  TERMS  -VE  INDICATES  Z-FORM 
1  I  1  -l 


Q3TU3M30  3a  OT 

I  *  A  23Y 


< 


*0*80331 


INPUT  DATA 


CONTROL  CARC  DATA  (FORMAT  2-413) 

2 4  12  1  C-iC  1  4  1  l  C  0  -2  01C0  0  2 

TITLE 

DIGITAL  CCNTRCLL ER  CESI6N 
A-  V  EC  TC  R  ENTRIES  (COEFFICIENT  DATA.  FCRMAI  9E8.5) 


-C.  1 

C.2 

1.0  0.0 

2.0 

4.  C 

1.0 

0.2 

l.C 

1.0  1.0 

0.  1 

-1.0 

1.0 

l.C 

l.C 

1.0  0.0 

1.0 

C.C 

1.0 

M- VECTC  R 

ENTRIES 

(FORMAT  2413) 

2  1 

1  C  2 

C  -1  4  1  1 

1  -1 

OTHER  FLAGS  AND  INPUT  DATA,  THE^E  AREt  N YU t OMEGA , DQMEGt 
CMEGF ,NBOC ,NZR T ,DT, T»  FMT( 1 )» STEP* LOOP, 

THESE  VARIABLES  ARE  ENTERED  AS  ABOVE  ACCORDING  TO  THE 
FOLLOWING  FORMAT  IF  NOT  APPLICABLE  TC  THE  PROBLEM 
NOTHING  NEEC  BE  ENTERED.  THE  CARDS  FUST  STILL  BE 
INCLLCEC. 

1ST  CARC _ I5,3E10.5.2I5.2£10.5,A6 

2ND  CARD  E1C.5,  15 

FOR  THIS  PROBLEM  THE  SUPPLEMENTARY  DATA  ENTRIES  ARE ♦ 
NZRT=1,  C  T= 1 . 0  »  T  =  6  C • 0  »  FMT(l)=MIN.,  STEP^l.O,  LCUP=1 

Ml  VECTOR  (FORMAT  1013) 


2  1  C  1  1 

A 1  VECTOR  (FORMAT  9E8.5) 

-C.55  l.C  12.0  -1.0  1.0 


FOR  INFORMATION  ON  LOOP  =  1  DATA  INPUT 
REFER  TC  SUBROUTINE  AM ALG »  APPENDIX 

A— 4 


CATA  END 


!  - 


A  TAG  Tjq/l 


(MAS  7  A  1 51  i  T  )  A  T  A  1  1HA3  JJHT/33 

S  0  k;  I  J  S  -  )  0  I  I  >  I  oi-::  i  r>  1  •>  S 

3JTI  T 

i  a  j  i  h 3  j  j  ) h  t  j  3  j  t  i  o  i  n 


i  .  b-ii  1  A  v*  i  1  .  TAJ 

T  'A  3 1  3  1  1 

1  i  ) 

;  ;  1  1  8  H3TD3V-A 

w.  !  J  .  A 

O.S 

0.0 

O.i 

S.O 

I  .0- 

S . t  .1  0.1- 

i  .0 

0.  1 

O.I 

3.1 

s.o 

0 .1  J .  3 

O.I 

0.0 

0 .  I 

3.1 

3.1 

(MAS  T  A  / '  0  T ) 

231  AT /3  HDTD3V-# 

i-  J  1 

I  A 

1-  3 

S  3  1  I 

s 

.  ■  i  t ,  < 

, J  HA  3hT 

.  A  T  A  0 

Tl>3/I 

)  /  A  aDAJI 

H33T3 

.  i  JJJtHJTa*  {  i  )  1  13,1  *  T 

, i h :  a* 

3  i  M,  333MO 

jHT  1  « i i<’  1  -  .  J3A  \  4 

iA  aA  33H3T/3  3Ha 

a  3  j  ‘j  a  r  -  a  v  j  a  j  h  t 

-n  •  iH'i  3i-U  JT 

;  J  8  A  3 I JHHA 

T  )M  31 

T  A  M  H  o  3 

0/ 1 W  3 J  J  33 

IJ11 

03H3TH3 

Jj  33  3/  3/MT3/1 

.  i  i  1  J  J D/l  I 

a  a 

,  •  . :  i  i (  » «  i  s 

t  <*'  .  0  l  .  . 

,  ?T 

J  *  A  3  T  a  1 

M  ,c  . 

D  n 

3HA3  3/S 

,  |  ,  A T A  i « AT l  i '  • 

3  3i  m 3 j •  ■ h q  am 

H33 

1=1  ,.)j  t  J.  1  ^riTa  ,  .HIM*  (  i  >TM3  ,0 

. 

,  0  .  I  =  T  3  1 1  » 

T  HS  / 

(El  J I  1  AMH03  )  fl  U  03  V  M 


i  I  0  I  s 
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F  1  NO*  CF  A-VECTOK  TERMS  N  *  24 


F  2 

NO.  f  F  M-VECTCR  TERMS 

N  = 

12 

F  3 

NO.  LF  RUNS  TO  BE  MADE. 

N  s 

l 

FA 

ONE  CF  FCUR  OPTIONS  SPECIFIED 

-IF  Z-TPANSFORM  TC  BE  COMPUTED  N=0 
-TWO  SAMPLER  SYSTEM.  N=N 

i 

-ALL  TRANS!  F  R  FUNCTIONS  IN  ONE  FORWARD 
LCOP  WITH  UNITY  FEEDBACK.  N=0 
-OTHERWISE  N=N 

N  = 

0 

F  5  ONE  CF  THREE  OPTIONS  SPECIFIED. 

~FCR  RCOT  LOCOS  OF  CONTINUOUS  SYSTEM  OR 
FCR  A  SYSTEM  IN  Z-FORM.  N=0 

_ _ HFC R  RCOT  LOCI  POINTS  OF  A  ONE-SAMPLER  _ N  =-LO 

SYSTEM  BUT  NO  Z-FURM  N=  L 
-FCR  Z-TRANSFORM  COMPUTED  OR  ROOT  LOCI 
PCINTS  FOR  TWO- SAMPLER  SYSTEM  N=-(10+...) 


_ F.6 _ HIGHEST  P.O WER  C F_S _ _ ..... _ N  =.  1 


F 7  HIGHEST  POWER  OF  Z  .  N  =  4 


F8  NUMBER  OF  VALUES  ASSIGNED  TO  T  _ N  =  1 


F  9 

ROOT  LOCI.  -t+tOR  BOTH 

FEEDBACK  OPTION 

(  1*3,0)  N 

a  1 

F  10 

SCAN  CONTROL  (N 

=0,1,- 

l)  V+H»  H » V 

'  N 

=  0 

FIl 

MODIFIED  Z-FORM 

OPTION 

( N  =  0,  l  ) 

'  N 

a  0 

F  12 

REPORT  HEADING 

OPTION 

(Na+,-2) 

N 

=  “2 

UNUSUAL  Z-FORM  OPTION  (N=i,2) 


I 


. J04M  38 


.  -I  '  \M  <  1  .  '  1 


(  |<  o:  30  -I  "'ll  *  ’I*  , 


:  -  Vl  ,H 
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F13  LOCI  OPT  ION i USUALLY  l|*0 _ _ _  N  *  0 

F 1-4  TERNS  IN  SERIES  FOR  G*(S),N=G  GIVES  19  TERMS  N  =100 

F  16  IF  Z-TRANSFORM  TO  13 E  COMPUTED  N.NE.O 

N=CEGREE  OF  RESULTING  Z-FORM  DENOMINATOR  N  =  2 

4N  RCCT  LOCUS  POINTS 

F  1 7  D-MATR l *  YES  N=1  N  =  -0 

F19  REAL  PART  N=-10K  N  =  -0 


M-VECTOR  DATA 


NO  FEEDBACK 


FCRWARC  L C CP  NO.  1 


NO.  CF  TERMS  IN  NUMERATOR  N  =  2 

DEGREES  OF-  THESE  TERMS  -VE  INDICATES  Z-FORM 
C  -1 

NO.  CF  TERMS  IN  DENOMINATOR  N  =  4 

DECREES  CF  THESE  TERMS  -VE  INDICATES  Z-FGRM 
1  I  1  -l 


iM«31  PI  23VI 0  0*M.(2)*D  m3  iil*32  2'«1 
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APPENDIX  A. 9. 7 

DIRECT  REPLACEMENT  OF  A  CONTINUOUS  CONTROLLER  BY 

A  DIGITAL  CONTROLLER 


Purpose : 

To  illustrate  the  method  outlined  in  Section  (7.4.2) 
for  the  design  of  digital  controllers. 

Transfer  Functions: 

Continuous  System  with  Controller 

Gc(s)Gp(s) 

Process  Transfer 

G  (s) 

P 

Process  Transfer 
Ho(s)Gp(s)  = 

Outline : 

The  first  step  in  the  design  is  to  derive  the  Z-transform 
of  the  continuous  system  with  the  continuous  controller  in¬ 
cluded.  The  root  locus  for  this  system  (Graph  16)  ,  is  also 
obtained  so  that  a  loop  gain  may  be  decided  upon  for  the 
transient  response  calculation.  The  associated  gains  are 
not  included  on  the  graph.  This  Z-transform 


s+0 . 1 

s  (s+1. 0)  (s+0 . 2) 
Function 

0.1 

(s+1 . 0)  (s+0 . 2) 

Function  with  Hold 

1  -  e“sT  0.1 

( - )  ( - ) 

s  (s+1 . 0) (s+0 . 2) 


. 


(S.fc.V)  nc  i.rj#;  i  i  ^9£±I3uo  bortt  n  uJ  s.^siieuIXi  oT 


.  >:  .  .  •:  '  ■ 


inlloi^aop  l  i  '  ^  oum  *n  O 


(S.O+e) (0. I+e)a 
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Graph  16 

s-Plane  Root  Locus  for  the  Z-Transform 
of  the  Continuous  System 


-• 


]6  4 


0. 59555678Z‘*  -  0.5391119Z  -  0.000269284 
(Z  -  0.81873076)  (Z  -  0 . 36787945)  (Z-l) 

was  derived  by  the  C.S.A.  program.  The  derivation  consti¬ 
tutes  one  run  of  the  program. 

The  second  run  is  needed  to  calculate  the  Z-transform 
of  the  process  transfer  function  and  hold.  This  is  the  term 

HG  (Z)  in  equation  (43) .  For  this  problem 
P 

0 . 034278Z  +  0.023014 

HG  (Z)  =  - 

P  (Z  -  0.81873)  (Z  -  0.36788) 

Using  these  two  transfer  functions,  the  Z-transform  of 
the  digital  controller  may  be  calculated  (Equation  45).  Thus, 

0.59557Z2  -  0 . 539112Z  -  0.000269 

D  (Z)  =  - - - 

(Z-l)  (0.034278Z  +  0.023014) 

This  is  the  Z-transform  expression  for  the  digital  con¬ 
troller.  It  is  not  necessary  to  carry  this  solution  through 
to  the  transient  response,  however,  for  illustration  this  is 
included. 

From  a  previous  problem,  not  included  in  this  appendix, 
the  transient  response  of  the  system  represented  as 

1.44  (s+0.1) 

G  (s)G  (s)  =  - - 

p  (s)  (s+1.0)  (s+0.2) 

is  available.  The  only  difference  between  this  system  and 
the  one  currently  being  treated  is  the  difference  in  gain. 
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Changing  the  process  transfer  function  numerator  from  0.1 
to  0.144  will  effectively  match  the  gains  since  a  gain  of  10 
is  already  included  in  the  digital  controller  transfer  func¬ 
tion.  With  the  gain  matching  completed  one  more  run  is  neces¬ 
sary  to  arrive  at  the  transient  response. 

This  particular  value  of  gain  was  picked  since  a  com¬ 
parison  was  available , however ,  different  values  can  be  used 
without  redesigning  the  controller.  It  is  only  necessary  to 
multiply  either  the  numerator  of  the  controller  or  that  of 
the  process  by  the  suitable  constant  to  give  the  gain  chosen 
from  the  Root  hocus  Plot,  Graph  (16) . 

The  input  data  for  each  of  the  three  runs  is  shown 
with  the  appropriate  titles.  The  root  locus  for  the  process 
under  continuous  control  (Graph  16)  is  shown.  The  root  locus 
for  the  system  under  the  equivalent  digital  control  was  cal¬ 
culated  as  a  check  for  error  but  was  the  same  as  the  con¬ 
tinuous  case  and  consequently  is  not  shown. 

Graph  (17)  gives  the  results  of  the  transient  response 
calculation.  This  response  was  identical  to  that  found  for 
the  continuously  controlled  system. 

In  Graph  (16)  a  horizontal  line  is  shown  at  Re(s) 

=  3.139.  This  line  marks  w  /2  which  is  the  sampling  fre- 

s 

quency  divided  by  two.  When  mapping  the  s-plane  into  the 
Z-plane  the  mapping  is  initially  carried  out  over  a  primary 
strip  in  the  left  half  of  the  s-plane.  The  rest  of  the 
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Graph  17 

Transient  Response  Continuous 
Controller  Replacement 
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s-plane  is  divided  into  an  infinite  number  of  periodic  strips 
of  width  w  which  also  map  into  the  unit  circle.  This  is 
shown  in  Kuo (12).  The  horizontal  line  in  this  plot  simply 
indicates  that  the  scan  has  been  carried  over  the  border  of 
the  primary  strip.  This  line  would  not  be  obtained  if  the 
root  locus  in  s  was  found  usinq  the  continuous  transfer  func¬ 
tion,  rather  than  the  Z-transform  equivalent. 

The  input  data  for  the  steps  in  this  problem  are  listed 
on  pages  168  to  173. 
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INPUT  DATA 


CONTROL  CARD  DATA  (FORMAT  2 A  I  3 ) 

24  12  I  C-1C  13  1000  -2  0  50  03 

TITLE _ 


Z-TRANSFCRM  CF  CONTINUOUS  SYSTEM 


A- VECTOR 

ENTRIES 

(COEFFICIENT 

DATA. 

FCRMAT 

9E8.5  ) 

-C  .  1 

0.2 

l.C 

0.0 

o 

• 

* 

4  •  C 

1. 

0 

0.2 

l.C 

1.0 

1.0 

-1.0 

l.C 

5. 

0 

0.72 

C.2 

1.0 

1.0 

1.0 

C.C 

!  . 

0 

M-VECTCR 

ENTRIES 

( FORMAT 

24  13) 

3-1  L=l  Q  2  I  Q  3  1  1  I 

OTHER  FLAGS  ANC  INPUT  DATA.  THESE  ARE,  NYQ , OMEGA , DOMEG, 

CMEGF,N80C,NZRT,DT,T,FMT(1) , STEP, LOOP. 

THESE  VARIABLES  ARE  ENTERED  AS  ABOVE  ACCORDING  TO  THE 
FOLLOWING  FORMAT  IF  NOT  APPLICABLE  TO  THE  PROBLEM 

_ NOTHING  NEED  BE  ENTERED.  THE  CARDS  MUST  STILL  BE _ 

INCLLCEC. 

1ST  CARC  15, 3E1C. 5,215, 2E10.5»A6 

2ND  CARD  E1C.5,I5 

FOR  THIS  PROBLEM  THE  SUPPLEMENTARY  DATA  ENTRIES  ARE, 

_ NLLL _ 


CATA  ENC 


ATAQ  TJM’/I 


(£]AS  TAMM  DU  i  TAJ  3  WAD  J3HT/13D 

i  k  )  s-  o  o  o  i  e  i  oi-o  i.  si  *s 

3JT  1  I 


:  . 

0.  I 

o.e 

o.  i 


I  A  •  VJ  3 

J  .  A 
D.  i 

D  .  J 


■  ' 

i  1303)  2  J  U  T  /1 3  W3TD3V-A 

J  .  A  0.0  '•i  *1  S  •  j  I  •  D~ 

o.l-  0.1  0 . 1  J.i  s .  0 

0.1  0.1  0.1  S.D  SY.O 

(  )  Z 

I  I  i  6  0  1  s  0  i-  I  I  c 


,  ,  r  •  AT  AO  T  0 q  ■"  1  1 1/1 A  2  D  A  J  T  '  1  ( 1 

.3  JJJUU7  2t  (I)TM3,7  ,  7  CUT  fl\/-.30S/.  3  33  0 
07  ;  /I  M  3.  A  JVOttA  2A  33fl3T/1J  3«A  2  3  J  3  A  1  WAV  32331 
M3JflO*H  3rt 7  37  3jyA31 JM^A  TOl*  31  TAMM33  OHIw3JJ33 
jjH  2  7  2  U v  2  0  M  A  0  3HT  .03M3T/3  3fl  333lA  D/U3TD/1 

.333JJ3/I 

oa.c.oi  iS,eis,e.oi3€,ei  3«ad  i 2 1 

21*2.313  o w a 3  3 /is 

231  AT  J  AT  A  3  YMAT^jM3JS3J2  33T  M3J6  3V9  21  3  T  03 


J  I  I/: 


J/13  A T  A 3 
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FI  NO.  CF  A- VEC TCR  TERRS 


N  *  24 


F2  NO.  CF  M-VECTGR  TERRS  N  *  12 


F 3  NO.  CF  RINS  TO  BE  RADE. 


F4  ONE  CF  FOUR  OPTIONS  SPECIFIED 

-IF  Z-TRANSFCRM  TC  BE  COMPUTED  N  =  0 
-TWO  SAMPLER  SYSTEM.  N=N 

-ALL  TRANSFER  FUNCTIONS  IN  ONE  FORWARD  N  -  0 

LOOP  WITH  UNITY  FEEDBACK.  N*0 
_-  C  T  H  E  R  W  I S  E  _  N  =  N  _ _ 


F 5  ONE  CF  TFREE  OPTIONS  SPECIFIED. 

-FCR  RCCT  LOCUS  CF  CONTINUOUS  SYSTEM  OR 
FCR  A  SYSTEM  IN  Z-FORM.  N=0 
-FCR  RCOT  LOCI  PCINTS  OF  A  ONE-SAMPLER  N  *-10 

SYSTEM  BUT  NO  Z-FORM  N=l 
-FCR  Z-TRANSFORM  COMPUTED  OR  ROOT  LOCI 
PCINTS  FOR  TWO-SAMPLER  SYSTEM  N=-(lO+...) 


F6  HIGHEST  POWER  OF  S _  N  *  l 

F 7  HIGHEST  POWER  CF  Z  N  =  4 


FR  NUMBER  OF  VALUES  ASSIGNED  TO  T _  N  =  1 

F 9  ROOT  LOCI.  -,+,OR  BOTH  FEEDBACK  OPTION  (1,3,0)  N  =  I 

F 10  SCAN  CONTROL  (N=0, 1,-1)  V  +  H,H,V  _ ;  N  =  0 

F 1 1  MODIFIED  Z-FORM  OPTION  (N  =  0,i)  N  *  0 

F 12  REPORT  HEADING  OPTION  (N  =  +,-2) _  __  N  *  -2 

UNUSUAL  Z-FORM  OPTION  (N=l,2) 


: 


wo  I  ’"u  Ov  10  H  1  '  t:*... 

tA3 1  1  SC  *  JS-  S  JAU2  ’  'U 
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I  13  LOCI  OPT  ION, USUALLY  N  =  0  .  N  =  0 

F  l  A  TERMS  IN  SERIES  FOR  G*(S),N=0  GIVES  19  TERMS  N  =  50 

F 16  IF  Z-TRANSFORM  TO  BE  COMPUTED  N.NE.O _ 

N=CEGREE  OF  RESULTING  Z-FORM  DENOMINATOR  N  *  3 

♦N  RCOT  LOCUS  POINTS 

F 1 7  B-MATRIX  YES  N  =  1  N  =  -0 

F 1 9  REAL  PART  N  =  -  IOK  N  =  -0 


M-VECTOR  DATA 


NO  FEEDBACK 


FORWARD  LOOP  NO.  I 


NO.  OF  TERMS  IN  NUMERATOR  N  =  2 

DEGREES  OF  THESE  TERMS  -VE  INDICATES  Z-FORM 
l  0 

NO.  CF  TERMS  IN  DENOMINATOR  N  =  3 

DEGREES  OF  THESE  TERMS  -VE  INDICATES Z-FORM 
1  1  1 


VJ' 


2m«i  pi  eavta  »03  Z3i«2  vu  mmt  pm 

■ 

I*M  23V  X I  UTAH-8 
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INPUF  DATA 


CONTROL  CARD  DATA  (  FORMAT  2 A  I  3 ) 

2 A  12  l  C - id  1  4  1  1  C  0-2  0100  0  2 


TITLE _ 

ANALOGUE  TO  DIGITAL  CONTROLLER 
A-VECTCR  ENTRIES  (COEFFICIENT  DATA.  FORMAT  9E8.5) 


-0.1 

C.2 

1.0  0.0 

2.0 

4  .  C 

1.00 

0.2 

l.C 

1.0  1.0 

0.  144 

-1.0 

1.0 

0.2 

1.0 

l.C 

1.0  0.0 

1.0 

C.C 

1.0 

M-VECTCR 

2  1 

ENTRIES 

1  0  2 

(FORMAT  2413) 

C  -1  4  1  1 

1  -1 

OTHER  FLAGS  AND  INPUT  DATA.  THESE  ARE,  N YQ , OMEGA , DGMEG , 
_____  CMEGF,NBOD,NZRT,DT,T,FMT(  1  ),  STEP,  LOOP. 

THESE  VARIABLES  ARE  ENTERED  AS  ABOVE  ACCORDING  TO  THE 
FOLLOWING  FORMAT  IF  NOT  APPLICABLE  TC  THE  PROBLEM 
NOTHING  NEED  BE  ENTERED.  THE  CARDS  MUST  STILL  BE 
INCLUDED. 

1ST  CARD  I5,3E10.5,2I5,2E10.5,A6 

2NC  CARC  E1C.5,  15 

FOR  THIS  PROBLEM  THE  SUPPLEMENTARY  DATA  ENTRIES  ARE, 

DT= 1 • 0  T  =  60 • C  *  FM  T (  1  ) =M I N • ,  STEP=1.0,  LC0P=1 


Ml  VECTOR  (FORMAT  1013) 
12  2  11 


A 1  VECTOR  (FORMAT  9E8.5) 

- •  C C  0  2  6  9-  • 5  39 112.59557  .023014  .034278  -1.0  1.0 

FOR  INFORMATION  ON  LOOP  =  1  DATA  INPUT 
REFER  TO  SUBROUTINE  AM ALG ,  APPENDIX 

A— 4 


CATA  ENC 


AT  AQ  IJ9/I 


(MAS  |  A  si  I  )  r  A  O  jT  1  10 

001 J  S-0  0  I  I  A  i  )  I  -  )  1.  Si  AS 


j  r  i  t 

P3JJOP7ZJO  JATI010  OT  3J03JAZA 
(  .  3<e  I  .w>  l  -  A  T  A  0  T/310M33  10)  20  I PT/3  POTOJm' 


00.  1 

Dm 

O.S 

0.0 

0.  1 

s.o 

I  .0- 

,  .  0 

0. 1 

.  t  - 

♦>AI  .0 

0.  1 

0.1 

0.1 

s.o 

0. 1 

0.3 

O.i 

0.0 

0.1 

0.1 

0.1 

<  €  I  AS  TAPPD3)  23 1  AT  43  A3T33V-* 


1-  1  i  i  A  I-  0  s  0  i  i  s 

i  MY/  M«A  323MT  .  AT  A  M  IOS/I  QUA  2DAJ3  P  3  H  T  3 

.  0  .  J  ,  S3  I  2  ,  (  I  )  T  H  ,  T  t  TO  »  T  PS  ' ,  0:  8/  •  00  iM  3 

0(7  u7  O/UPDOOA  3  V  0  a  A  2A  33P3T/I3  3PA  3J3AIPAV  3230T 

MJQlPS  (T  JT  3 J 0 A 0 1  J  3  9  A  T 314  31  TA’*  00  0/M3JJ03 
ji  J  j  1  T  i.  Tc  04  20  A  )  jHT  .  C.  3  503  MM  (i  133/  0/  I  TO/ 

.333JJ0/I 

dAfe.0i3St5ISte.0I J(  ,ei  3PA0  T  2  I 

<M  ,r  .01  1  3PA0  J/S 

tOP  2JJPT/3  ATAJ  YPAT/3M3 JH9U2  30T  P3J30PS  2IOT  P03 

,  .  |  1  .  .  i  -(ill  •  « 3 . 0  d  *  1  . 1  *  T  0 

(  £  1  0 J  T  A M 5303  )  P0T03V  1  M 

i  1  S  S  1 

(M9M  f  A MPG3 )  P3T03V  IA 
.1-  V  S  A  0  •  -  'IOtS3.  VgM  .SliM  1  .  -  /  d  S  0  0  0  .  - 

TUS/1  A  1  AO  I  =  SOOJ  140  ZG  I  T  A ‘OP  03  /  I  PG3 
XI  /MSA  » <  JAMA  OHITUOPU2  OT  P33J- 

A-A 


3/3  A  T  A3 
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FI  MO*  OF  A-VECTOR  TERMS  N  *  24 


F 2  NO.  OF  M-VECTGR  TERMS  N  *  12 

F 3  NO.  OF  RUNS  TO  BE  MADE.  N  =  1 


F 4  ONE  OF  FOUR  OPTIONS  SPECIFIED 

-IF  Z-TRANSFORM  TO  BE  COMPUTED  N=0 
-TWO  SAMPLER  SYSTEM.  N=N 

-ALL  TRANSFER  FUNCTIONS  IN  ONE  FORWARD  N  =  0 

LOOP  WITH  UNITY  FEEDBACK.  N=0 
-OTHERWISE  N  =  N _ 


F5  ONE  OF  THREE  OPTIONS  SPECIFIED. 

-FOR  ROOT  LOCUS  OF  CONTINUOUS  SYSTEM  OR 
FOR  A  SYSTEM  IN  Z-FORM.  N=0 

-FOR  ROOT  LOC I  POJNTS  OF  A  ONE- SAMPLER _  N  =-10 

SYSTEM  BUT  NO  Z-FORM  N=l 
-FOR  Z-TRANSFORM  COMPUTED  OR  ROOT  LOCI 
POINTS  FOR  TWO-SAMPLER  SYSTEM  N=-(10+...) 

F 6  HIGHEST  POWER  OF  S _ _ N  • . 1 

F 7  HIGHEST  POWER  OF  Z  N  *  3 

F 8  NUMBER  OF  VALUES  ASSIGNED  TO  T  _  N  =1 


F9  ROOT  LOCI.  -,+,OR  BOTH  FEEDBACK  OPTION  (1,3,0)  N  =  0 

F 10  SC  AN  C  ON  JR  0  L (N«0»  1  t-l )  '  V+Ht.H»  V .  . N  =  0 

Fll  MODIFIED  Z-FORM  OPTION  (N=0,1)  N  *  0 

F 12  REPORT  HEADING  OPTION  _(N«+j-2)  N  «  -2 

UNUSUAL  Z-FORM  OPTION  (N=l ,2) 


•  3 Q A M  01  • 


J  t  t  '  J  j  V  ....  .. 


■.  32 1  »<»3nT0- 


-  17  3  - 

F  1  3  LCC  I  CPTICN, USUALLY  N  =  0  N  =  0 

F-H  TERNS  IN  SERIES  FC'R  G*(S),N=0  GIVES  IV  TERMS  N  =100 

E  16  IF  Z-TRANSFORN  TO  BE  COMPUTED  N.NE.G 

N  =  C  EGR  E  E  CF  RESULTING  Z-FORM  DE  IMEjN  I  N  A  TCR.  N  =  2 

+  N  RCCI  LOCUS  POINTS 

F 1 7  B-NATRIX  YES  N=1  N  =  -0 

F  1 9  REAL  PART  N  =  -  ICK  N  =  -0 

M-VFCTGR  DATA 


NO  FEEDBACK 


FORWARD  LCCP  NC.  l 

NO.  CF  TERNS  IN  NUMERATOR  N  =  2 

DECREES  OF  TEESE  TERMS  -VE  INDICATES  Z-FORM 
C  -I 

NC.  CF  TERMS  IN  DENOMINATOR  N  =  4 

DEGREES  CF  TEESE  TERMS  -VL  INDICATES  Z-FORM 

1  I  1  -1 


' 
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